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Abstract 

This review is focused on neutrino mixing and neutrino oscillations in the light 
of the recent experimental developments. After discussing possible types of 
neutrino mixing for Dirac and Majorana neutrinos and considering in detail the 
phenomenology of neutrino oscillations in vacuum and matter, we review all 
existing evidence and indications in favour of neutrino oscillations that have 
been obtained in the atmospheric, solar and LSND experiments. We present 
the results of the analyses of the neutrino oscillation data in the framework of 
mixing of three and four massive neutrinos and investigate possibilities to test 
the different neutrino mass and mixing schemes obtained in this way. We also 
discuss briefly future neutrino oscillation experiments. 
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1 Introduction 



The strong evidence in favour of oscillations of atmospheric neutrinos found by the Super- 
Kamiokande Collaboration ||T],|2|] opened a new era in particle physics. There is no doubt 
that new experiments are necessary to understand the nature of neutrino masses and 
mixing which are intimately connected with neutrino oscillations, but the first decisive 
step has been done: massive and mixed neutrinos can now be considered as real physical 
objects. 

The problem of neutrino mass has a long history. Originally, Pauli considered the 
neutrino as a particle with a small but non-zero mass (smaller than the electron mass) 
and the method for the measurement of the neutrino mass through the investigation of 
the /9-spectrum near the end point was proposed in the first theoretical papers on /3-decay 
of Fermi [|] and Perrin || . 

The first experiments on the measurement of the neutrino mass, based on the Fermi- 
Perrin method, yielded the upper bound m v < 500 eV |7j which was improved in the fifties 
to m v < 250 eV ||. Therefore, it became evident that the neutrino mass (if non-zero at 
all) is much smaller than the electron mass. This was the main reason that in 1957, after 
the discovery of parity violation in /3-decay, the authors of the two-component theory of 
the neutrino (Landau ||, Lee and Yang [HJ, Salam [ll]) assumed that the neutrino is a 
massless particle, the field of which is either a left-handed field vl or a right-handed field 

In 1958, Goldhaber et al. JT2J measured the helicity of the neutrino. The result of 
this experiment was in agreement with the two-component neutrino theory and it was 
established that the neutrino field is The results of the experiment of Goldhaber et 
al. could not exclude, however, the possibility of a small neutrino mass. In the V— A 
theory (Feynman and Gell-Mann [lq|, Sudarshan and Marshak [|TI|]) the Hamiltonian of 
weak interactions contains the left-handed component of the neutrino field vl, and also 
the left-handed components of all massive fields. Therefore the possibility for the neutrino 
to be nevertheless a massive particle became more natural [|HJ after the confirmation of 
the V— A theory. 

In 1957, B. Pontecorvo JTB,TP] proposed the idea that the state of neutrinos produced 



in weak interaction processes is a superposition of states of two Majorana neutrinos f20|| 
with definite masses (analogous to the states \K°) and \K°) which are the superposition 
of \Ki) and \K 2 ), the states of particles with definite masses and widths). In this way, 
B. Pontecorvo arrived at the hypothesis of neutrino oscillations (analogous to K° ^ K° 
oscillations). At that time only one type of neutrino was known. The possibility of mixing 
of the two species of neutrinos v e and was considered in Ref. pip . All possible types 
of neutrino oscillations for this case were investigated by Pontecorvo in 1967 p^j . 

Gribov and Pontecorvo proposed in 1969 JSJ the first phenomenological theory of 



neutrino mixing and oscillations. In this theory, the two left-handed neutrino fields v e ^ 



Notice that in the Goldhaber et al. experiment the helicity of the electron neutrino was measured. 
The helicity of the muon neutrino was measured in several experiments (for the references see the review 
of V.L. Tclcgdi fll3|). The best accuracy in the measurement of the muon neutrino was achieved in the 
experiment by Grenacs et al. Hj. 
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and i/pL are linear combinations of the left-handed components of the fields of Majorana 
neutrinos with definite masses and the neutrino mass term contains only the left-handed 
fields v e L and v^. 

In 1976, neutrino oscillations were considered in the scheme of mixing of two Dirac 



neutrinos based on the analogy between quarks and leptons pi,p5|| and in the same year 
in the general Dirac-Majorana scheme |26j (for later works see Refs. [f27j-|30"|1). 



The theoretical arguments in favour of non-zero neutrino masses and mixing are based 
on the models beyond the Standard Model (see, for example, Ref. ||31|| ). In such models 
the fields of quarks, charged leptons and neutrinos are grouped in the same multiplets and 
the generation of the masses of quarks and charged leptons with the Higgs mechanism as 
a rule provides also non-zero neutrino masses. 

In 1979, the see-saw mechanism for the generation of neutrino masses was proposed 
1 32 -[34 . This mechanism connects the smallness of neutrino masses with the possible 
violation of lepton number conservation at a very large energy scale. 

At present the effects of neutrino masses and mixing are investigated in many different 
experiments. There are three types of experiments in which the effects of small neutrino 
masses (say, of the order of 1 eV or smaller) and mixing can be revealed (for a review 
and references see Ref. 



1. Neutrino oscillations experiments. 

2. Experiments on the search for neutrinoless double /3-decay. 

3. Experiments on the measurement of the electron neutrino mass with the precise in- 
vestigation of the high energy part of the /3-spectrum of 3 H. 

Three indications in favour of neutrino masses and mixing have been found so far. 
These indications were obtained in the following experiments: 

1. Solar neutrino experiments (Homestake [j36H38|1, Kamiokande [^9|-^l|, GALLEX |42 
|43|, SAGE g§@, Super-Kamiokande 



2. Atmospheric neutrino experiments (Super-Kamiokande []1],[§, Kamiokande [|9|, 1MB 



H, Soudan fTLfl, MACRO ||]). 
3. The accelerator LSND experiment p3j|54 



Many other neutrino oscillation experiments with neutrinos from reactor and accel- 
erators did not find any evidence for neutrino oscillations. In the experiments on the 
search for neutrinoless double /3-decay no indications for non-zero neutrino masses were 
found (see Section 6.2). The present upper bound for the electron neutrino mass ob- 
tained in the Troitsk experiment [55] is 2.7 eV (see also the Mainz experiment f56|). The 
upper limits on the masses of and v T are 170 keV (90% CL) and 18.2 MeV (95% CL), 
respectively |j35 |. Neutrinos play an important role in cosmology and astrophysics and 
many bounds on neutrino properties can be derived in this context. For reviews see, e.g., 
Refs. |31|,|57 . 
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In this review we discuss the phenomenological theory of neutrino mixing (Section 
2), neutrino oscillations in vacuum (Section 3), neutrino oscillations and transitions in 
matter (Section 4) and experimental data and results of analyses of the data (Sections 5 
and 6). We also consider the implications of the existing experimental results on neutrino 
oscillations for experiments in preparation. After the conclusions (Section 7) we discuss 
some properties of Majorana neutrinos and fields in Appendix A. 

We hope that this review will be useful not only for the physicists that are working 
in the field but also for those who are interested in this exciting field of physics. In many 
cases we present not only results but also derivations of the results. For those who start 



to study the subject we refer to the books Refs. [31, 58 - SI and the reviews Refs. [62-73 



2 Neutrino Mixing 

All the numerous data on weak interaction processes are perfectly well described by the 



Standard Model |74H76|. The standard weak interactions are due to the coupling of quarks 
and leptons with the gauge W and Z vector bosons, described by the charged- current 
(CC) and neutral-current (NC) interaction Lagrangians 

C™ = --±=jf°W» + h.c, (2.1) 

£ NC = __ff.NC ZP< (2 _ 2) 



2cos# 



ir 



Here g is the SU(2)l gauge coupling constant, 9w is the weak angle and the charged and 
neutral currents and are given by the expressions 



•cc 



■NC 
Jp 



2 ^7p4 + -- - , (2.3) 

£=e,fj.,T 

^2 ^^pVeL + • • • , (2.4) 



-e,p,T 



where the I are the physical charged lepton fields with masses mg and we have written 
explicitly only the terms containing the neutrino fields. The flavour neutrinos z/ e , v^,... 
are determined by CC weak interactions: for example, the is the particle produced 
in the decay n + — > fi + + and so on. The number of light flavour neutrinos is given 
by the invisible width of the Z boson J77| in the Standard Model. It was shown in the 
famous LEP experiments on the measurement of the invisible width of the Z boson that 
the number of the light flavour neutrinos is equal to 3. The most recent experimental 
value of the number of neutrino flavours is 2.994 ± 0.012 |[35| , showing that are no other 
neutrino flavours than the well-known u e , u^, u T . 

The CC and NC interactions conserve the electron L e , muon and tau L T lepton 



numbers, which are assigned as shown in Table p7T . 

There are no indications in favour of violation of the law of conservation of lepton 
numbers in weak processes and very strong bounds on the probabilities of the lepton 



5 















e") 


+ 1 















+1 





{v T 










+ 1 



Table Assignment of lepton numbers. The corresponding antiparticles have opposite lepton numbers. 



number violating processes have been obtained from the experimental data. The most 
stringent limits (90% CL) are (see Ref. JHJ): 

T(/i -> e 7 )/r(/i -> all) < 4.9 x 10 -11 , (2.5) 

r(/i -> 3 e)/r(/i -> all) < 1.0 x 10~ 12 , (2.6) 

a{n~ Ti -> e~ Ti)/cr(/x~ Ti -> capture) < 4.3 x 10~ 12 , (2.7) 

r(fsT L -> e n)/T(K L -> all) < 3.3 x 10" 11 , (2.8) 

r(K+ -> 7r + e~/i + )/r( J RT + -> all) < 2.1 x 10~ 10 . (2.9) 

According to the neutrino mixing hypothesis fl8],|l9", 21], the conservation of the lepton 



numbers is only approximate. It is violated because of non-zero neutrino masses and 
neutrino mixing. In the case of neutrino mixing, the left-handed flavour neutrino fields 
v iL are superpositions of the left-handed components v kL of the fields of neutrinos with 
definite masses m^: 



k=l 



Ik VkL 



{£ = e,/i,r) , 



(2.10) 



where U is a unitary mixing matrix. The number of massive neutrino fields n could be 
equal or more than 3. If n is larger than 3 there are sterile neutrinos that do not take 
part in the standard weak interactions ( |2.3| ) and ( |2.4| ). 

It is well-established that quarks take part in CC weak interactions in mixed form 
with the V— A current 

q'=u,c,t q=d,s,b 



where V is the Cabibbo-Kobayashi-Maskawa mixing matrix |f7S, 75[]. The relation ( |2.10D 
is analogous to the mixing in Eq. (|2.11| ). However, between the mixing of neutrinos and 
quarks there can be a fundamental difference. In fact, quarks are four-component Dirac 
particles: quarks and antiquarks have opposite electric and baryonic charges. Instead, 
neutrinos are electrically neutral particles. If the total lepton charge 

L = L e + L^ + L T (2.12) 

is conserved, neutrinos with definite masses are four-component Dirac particles like quarks 
(in this case a neutrino differs from an antineutrino by the opposite value of L). If the 
total lepton number (|2.12| ) is not conserved, massive neutrinos are truly neutral two- 
component Majorana particles. These possibilities are realized in different models and 
correspond to different neutrino mass terms. 
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2.1 Dirac mass term 



A Dirac neutrino mass term can be generated by the Higgs mechanism with the standard 
Higgs doublet which is responsible for the generation of the masses of quarks and charged 
leptons.0 In this case the neutrino mass term is given byfj 

£ D = -^%M° v eL +h.c. (£ = e,/i,r), (2.13) 

£,£' 

and M D is a complex 3x3 matrix. 

The Dirac mass term ( |2.13| ) can be diagonalized taking into account that the complex 
matrix M D can be written as 

M D = VmU\ (2.14) 
where V and U are unitary matrices and m is a positive definite diagonal matrix, m^j = 



rrikSkj, with > 0. Hence, the Dirac mass term ( 2.13 ) can be written in the diagonal 
form 



3 



£ D = - ^ m k v k Vk , (2.15) 



k=l 



with 



3 

J2 U ^^kL (£ = e,/i,r). (2.16) 

k=l 

Therefore, in the case of the Dirac mass term ( |2.13j ) the three flavour fields (i = 
e, n, t) are linear unitary combinations of the left-handed components Vki of three fields 
of neutrinos with masses (k = 1, 2, 3). On the other hand, we also have 

3 

VtR = ^2 Vek VkR ^ = e ' /i ' T ) ' ( 2 - 17 ) 
k=l 

but these right-handed fields do not occur in the standard weak interaction Lagrangian. 
Therefore, right-handed singlets are sterile and are not mixed in this scheme with the 
active neutrinos.^ 

The field is a Dirac field if not only the mass term ( |2.13| ) but also the total La- 
grangian is invariant under the global U(l) transformation 

ui^e^ut, £^e iv £ (£ = e,/i,r), (2.18) 



2 One needs to assume the existence of right-handed SU(2)l singlet fields vgR. If these fields do not 
exist, it is not possible to construct a Dirac mass term for the neutrinos, i.e., neutrinos are massless 



particles. This model is sometimes called Minimal Standard Model [|74|~|76|. 

3 Note that in this review we use four-component spinors. For the equivalent alternative of two- 
component spinors see, e.g., Refs. p0|,pt|. 

4 In the Dirac case it is possible to conceive also the existence of left-handed sterile singlet fields mixed 
with the active flavour neutrino fields, in contrast to the right-handed singlets which cannot be mixed 
with the active neutrinos because of lepton number conservation. 
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where the phase <p is the same for all neutrino and charged lepton fields. Then, using 
Noether's theorem, one can see that the invariance of the Lagrangian under the trans- 
formation (|2.18|) implies that the total lepton charge L is conserved. Therefore, L is the 
quantum number that distinguishes a neutrino from an antineutrino. 

The Dirac mass term ( |2.13| ) allows processes like /i — ► e + 7, fi~ — > e~ + e + + e~ . 



However, the contribution of neutrino mixing to the probabilities of such processes is 
negligibly small H81-83T1. 



The unitary 3x3 mixing matrix U can be written in terms of 3 mixing angles and 6 



phases. However, only one phase is measurable [ff9| . This is due to the fact that in the 
Standard Model the only term in the Lagrangian where the mixing matrix U enters is 
the CC interaction Lagrangian fl2.1|) . With neutrino mixing the lepton charged- current 
is given by 

3 

J cct = 2 J2 h lp u eL = 2 ^hlpUikU^. (2.19) 

£=e,fj.,T £=e,fj,,r fc=l 

Because the phases of Dirac fields are arbitrary one can eliminate from the mixing matrix 
U five phases by a redefinition of the phases of the charged lepton and neutrino fields, 
with only one physical phase remaining in U. The presence of this phase causes the 
violation of CP invariance in the lepton sector. A convenient parameterization of the 



mixing matrix U is the one proposed in Ref. | 8"3}| 



C12C13 S12C13 s 13 e lSls 

U = ( -S12C23 - ci2S 2 3Si 3 e l<513 C12C23 - s 12 s 2 3S 13 e lSl3 S23C13 I , (2.20) 

S12S23 - Ci2C23Si 3 e J<513 -C12S23 - Si 2 C23Si 3 e l<513 C23C13 

where Cy = cos i?^- and = sin^y and 8i 3 is the CP- violating phase. 

Since in the parameterization ( p.20| ) of the mixing matrix the CP-violating phase ^13 



is associated with s 13 , it is clear that CP violation is negligible in the lepton sector if 
the mixing angle $13 is small. More generally, it is possible to show that if any of the 
elements of the mixing matrix is zero, the CP-violating phase can be rotated away by a 
suitable rephasing of the charged lepton and neutrino fields. 5 



2.2 Dirac— Majorana mass term 

If none of the lepton numbers is conserved and both, left-handed flavour fields (£ = 
e, [A, t) and sterile right-handed gauge singlet fields u s r (s = S\,S2,...) enter into the 
mass term we have the so-called Dirac-Majorana mass term 

£ D+M = C M + £ D + C M ^ (2 21) 

with 

CP = -^7^M s >, L + h.c., (2.22) 

s,£ 

5 This can also be seen by noticing that the Jarlskog rephasing-invariant parameter ]85|-|87[| vanishes 
if one of the elements of the mixing matrix is zero. 
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£ 



£ 



i^l^M* (^) c + h.c. (2.24) 



Here M D , M L and M R are complex matrices and the indices s, s' run over n R values. 
Let us stress that the number of right-handed singlet fields could be different from the 



number of neutrino flavours (see [30,8 



The charge-conjugate fields are defined by 

{v lL ) c = CVE T , (u sR ) c = Ci7^ T , (2.25) 

where C is the charge conjugation matrix (see Appendix A). Notice that {yih) c and {v s r) c 
are right-handed and left-handed components, respectively. Indeed, with C~ lr y 5 C = 7 J 
one finds 

^^L) c = c(^4p) =Cve T = {^l) c (2.26) 
and the analogous reasoning holds for {v s r) c - Furthermore, we have 



-uJ L C-\ (v sR y = -vJrC' 1 . (2.27) 



The matrices M L and M R are symmetric. This can be shown with the help of the 
relation vJ l C~ 1 i>£il = vJ, l C~ 1 V£l, which follows from the fact that C is an antisymmetric 
matrix and from the anticommutation property of fermion fields. An analogous relation 
holds for the right-handed fields. Then, using Eq.( [2.27|) , one obtains 

("il) c Mf e v tL = ~J2 v t i>l K< nL = E fa) M ^ Ui ' L ■ ( 2 - 28 ) 

£,£' £,£' £,£' 

Let us introduce the left-handed column vector 

n L = ( , V \ c ) with v L = I v^ L I , (2.29) 




where u R is the column vector of right-handed fields. With the relation 

Yf* m * v * = - E "a M ° ^ T = E ( MD )£ (^) C ' ( 2 - 3 °) 

s,£ s,£ s,£ 

the Dirac-Majorana mass term £ D+M can be written in the form 

1 



£ D+M = (nL)c M D+M ^ + h c ^ (231) 



with the symmetric (3 + n R ) x (3 + n R ) matrix 



S ( £ %I ) • (2.32) 
9 



The complex symmetric matrix M D+M can be diagonalized with the help of a unitary 
matrix U H,l0|: 

M D+M = ([/t) T mf/t, (2.33) 

with rhkj = m^Skj and > 0. Using this relation the mass term ( |2.31| ) can be written 

as 

C D+u = --NmN=-- J2 m ^ u k, (2.34) 



2 2 

k=l 



with 

N 



\ [ ) 

The fields satisfy the Majorana condition 



17+ n L + (U ] n L ) c . (2.35) 



(z/ fc ) c = z/ fc (A; = l,...,3 + n fl ). (2.36) 

Therefore, in the general case of a Dirac-Majorana mass term the fields of particles with 
definite masses are Majorana fields. Majorana particles are particles with spin 1/2 and 
all charges equal to zero (particle = antiparticle). Some properties of Majorana fields are 
discussed in Appendix A. 

We want to emphasize that it is natural that the diagonalization of the mass term 
( |2.31| ) leads to fields of Majorana particles with definite masses: the mass term ( |2.31| ) for 



the case of a general matrix M D+M is not invariant under any global phase transformation. 
In other words, in the general case of the Dirac and Majorana mass term there are no 
conserved quantum numbers that allow to distinguish a particle from its antiparticle. 
From Eq.(|2.35|), for the left-handed components of the neutrino fields we have the 



mixing relations 

3+n R 

vtL = ^2 Uik VkL ^ = e ' ^' r ) ' (2.37) 

k=l 
3+n R 

(vsrY = Usk VkL ■ ( 2 - 38 ) 

k=l 

Thus, in the case of a Dirac-Majorana mass term the flavour fields un are linear uni- 
tary combinations of the left-handed components of the Majorana fields of neutrinos 
with definite masses. These components are also connected with the sterile fields (v s r) c 
through the relation ( |2.3S ). If the masses rrik are small, the mixing relations (|2.37|) and 



( |2.38| ) imply that flavour neutrinos u e , and v T can oscillate into sterile states that are 
quanta of the right-handed fields v sR . We will discuss such transitions in Section 6. 

Let us stress that in order to have all three terms in the Dirac-Majorana mass term 
( p.31| ) not only a Higgs doublet but also a Higgs triplet [|9lH 93H and a Higgs singlet 
are necessary. Thus, it can be generated only in the framework of models beyond the 
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Standard Model. A typical example is the SO(10) model (see, for example, For a 

discussion of radiative corrections to the Dirac-Majorana mass term see Ref. |M . 



Up to now we did not make any special assumption about the Dirac-Majorana mass 
term (|2.31|) . Let us now consider the possibility that CP invariance in the lepton sector 
holds. In this case we have 

U CP C D+M (x) = C D+M (x P ) . (2.39) 

where Ucp is the operator of CP conjugation, x = (x°,x) and xp = (x°, —x). For the 
fields ni(x) we take the usual CP transformation 

U CP n L {x) Ucl = r 11 °CnZ T {xp) , (2.40) 

where rj is a diagonal matrix of phase factors. The mass term £ D+M can be written in 
the form 

C D+M (x) = ]pl{x)C- l M* +u n L {x) - ^nZ(x) (M D+M )t Cn?(x) . (2.41) 

Using Eq. ( [2.40| ) we obtain 

Ucp n T L (x) C~ l M D+M n L (x) = nZ(x P ) T 1 M T ' +M T 1 CnZ T {x F ) . (2.42) 

Hence, the CP invariance condition ( p.39| ) is satisfied if 

r/ M D+M r ? = -(M D+M ) t . (2.43) 

Up to now 7] was an arbitrary diagonal matrix of phase factors. If we choose rj = i, in 
this case CP invariance in the lepton sector implies that 

M D+M = ( M D+M )t . (2.44) 

Taking into account that M D+M is a symmetric matrix, the condition ( |2.44j ) is equivalent 
to 

M D+M = (M D +M) * (245) 

The real symmetric matrix M D+M can be diagonalized with the transformation 

M D+M = Om'0 T , (2.46) 

where O is an orthogonal matrix [O t = 0~ r ) and m! is a diagonal matrix, m' k - = m' k 5kj- 
The eigenvalues m' k can be positive or negative and the neutrino masses are thus given 
by ink = \ m 'k\- Therefore, we write the eigenvalues as 

m k = m k p k , (2.47) 

where p k = ±1 is the sign of the k th eigenvalue of the matrix M D+M . Then the relation 
( p.46| ) can be rewritten in the form 

M°+ u = (rf) T mrf (2.48) 
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with 

C/t = ^pQ T (2.49) 

and rhkj = m k bkj- Here p is the diagonal matrix with elements pkj = Pk$kj- Therefore, if 
CP is conserved in the lepton sector, the neutrino mixing matrix has the simple structure 
shown in Eq.( |2.49| ), which implies that 

U* = Up. (2.50) 

The CP transformation of the Majorana field v k is given by (see Appendix A) 

U C p v k {x) U C p = r§ p 70 v k (x P ) , (2.51) 

where r]^ p is the CP parity of the Majorana field v k which is 7/{? p = ±i (see Ref. [|35|j). 
From Eq. ( |2.51| ) we have 

U C p v kL (x) U C p = r]^ p 70 v kR {x P ) . (2.52) 

This relation and Eqs.( |2.35| ), (|2.40|) and ( |2.50| ) imply that the CP parity of the Majorana 
field v k is 

Vk P = *Pk- (2-53) 

Indeed, we have 

Ucp N L (x) U c l = UU l0 C nZ T (x P ) =U^U*t lo C \x P ) =ip lo N R (x P ) . (2.54) 
2.3 Majorana mass term 

If only left-handed neutrino flavour fields [I = e, p, r) enter into the Lagrangian we 



can write down the mass term [23. 98 



CM = -\H M C M «' v ^ > ( 2 - 55 ) 

l,t' 

where M L is a symmetric complex matrix (see previous section). Then we have for the 
mixing 

3 

vtL = ^2 Uek VkL ' 

k=l 

where v k is the field of Majorana neutrinos with mass rrik- In this case the number of 
massive Majorana neutrinos is equal to the number of lepton flavours (three). Note that 
the generation of the Majorana mass term ( |2.55| ) requires an enlargement of the scalar 
sector of the Minimal Standard Model p9f : with a Higgs triplet like in the Majoron 
models Ref. ]9~1"|-|9"3]] Majorana masses are obtained at the tree level, whereas radiative 
generation is possible at the 1-loop level with a singly charged Higgs singlet (plus an 
additional Higgs doublet) |3"T|, |1UU| - [1U4| |) or at the 2- loop level with a doubly charged 
Higgs singlet (plus an additional singly charged scalar) ||105| . 
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Since the Majorana condition ( |2.36| ) does not allow the rephasing of the neutrino 
fields, only three of the six phases in the 3x3 mixing matrix U can be absorbed into 
the charged lepton fields in the charged current ( |2.19| ). Therefore, in the Majorana 
case the mixing matrix contains three CP- violating phases [28, [29|, |106|| in contrast to 
the single CP-violating phase of the Dirac case discussed in Section |2.1|. However, the 



additional CP-violating phases in the Majorana case have no effect on neutrino oscillations 
in vacuum p8|,p9|,|T06| (see Section 0) as well as in matter ||107|| . 



2.4 The one-generation case 



Let us consider the Dirac-Majorana mass term in the simplest case of one generation. 
We have 



C 



D+M 



1 



1 



m L {v L ) c v L - m D v R v L - - m R v R [v R ) c + h.c. 



~ m R (n L ) c Mn L + h.c. 



with 



M = 



m L m D 
m D m R 



(2.57) 



(2.58) 



For simplicity we assume CP invariance in the lepton sector (see Subsection 2.2). In 
this case mi, and m R are real parameters. In order to diagonalize the matrix M, let 
us write it in the form 



M = - TrM + M° = - (m L + m R ) + M° 
2 2 



with 



^ ~\ ( m R - m L) 



m D 



(2.59) 



^2.60) 



The eigenvalues of the matrix M° are 



m i,2 = =F^ \l (rn R - m L ) 2 + 4m^ . 



Furthermore, we have 

M° = Om O T , 
where m° = diag(m5,m2) and O is the orthogonal matrix 



O 

with the mixing angle $ given by 

m R - m L 



cos 2$ 



^(m R -m L ) 2 + 4ml 



cos $ sin $ 
■ sin $ cos $ 



sin 2d 



2m 



D 



y/(m R - m L ) 2 + 4m 2 j 



(2.61) 
(2.62) 

(2.63) 
(2.64) 



13 



Hence, the matrix M can be written as 



M = Orri O t 



(2.65) 



where ml = diag(m' 1 , m' 2 ) and 



m 



1,2 



(m R + m L ) =F J {m R - m L ) 2 + 4m^ 



(2.66) 



The eigenvalues of the matrix M are real but can have positive or negative sign. Let us 
write them as 



m k = m k p k , 



(2.67) 



where 



m. 



and pk is the sign of the k eigenvalue of the matrix M. As shown 



in Eq.( [2.53| ), the CP parity of the Majorana field i/& with definite mass uik is vfc 
The relation ( |2.65| ) can be written in the form 

M = (U ] ) T m 

with 77i = diag(mi, m 2 ) and 



CP 



«Pfc- 



(2.68) 



(2.69) 



with p = diag(pi,p 2 ). Using now the general formulas obtained in Section |2~2| , we have 
the mixing relation 



U 



"2L 



(2.70) 



with 



U 



Pi) 



(2.71) 



cost? (\fp~2)* suit? 
sin 7? (v 7 /^)* cost? 

Therefore, the three parameters mi, 77?d, 77?^ are related with the mixing angle d 
and the neutrino masses mk by the relations ( |2.64j ) and fl2.66| ), ( |2.67| ). The signs of the 
eigenvalues of M determine the CP parities of the massive Majorana fields Uk- 

In the framework of CP conservation, the relations obtained so far are general. In 
the following part of this section we consider some particular cases with special physical 
significance. 



2.5 The see-saw mechanism 



Let us consider the Dirac-Majorana mass term (|2.57|) and assume |]32|-p4| that = 0, 



TTiD — 7??/, where 77?/ is the mass of a quark or a charged lepton of the same generation, 



and m R ~ M. ^> m' . In this case, from the relations ( |2.64j ), ( |2.66| ) and ( j2.67|) we have 

(2.72) 

< m f , pi = -1 , 



n m 

7? ~ — < 1 

M 



777 



7771 ~ 

7772 — 



M 

M, 



P2 



Pi 
1 . 



(2.73) 
(2.74) 
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These relations imply that approximately 



\VR) 



V-IL 



(2.75) 



and the Majorana fields v\, v 2 are connected with the fields vl an d vr by 



v x =i[u L - {v L ) c \ , u 2 = u R + (u R ) e . (2.76) 

The mechanism which we consider here is called see-saw mechanism []32"l-j34j1 . It is 
based on the assumption that the conservation of the total lepton number L is violated 
by the right-handed Majorana mass term at the scale Ai that is much larger than the 
scale of the electroweak symmetry breaking. Several models which implement the see- 
saw mechanism are possible (see, for example, PI, 66], |68]| and references therein) and 



the scale A4 depends on the model. This scale could be as low as the TeV scale (for 
example, in left-right symmetric models [31 , 108 , 109|| ) or an intermediate scale, or as high 
as the grand unification scale ~ 10 15 GeV or even the Plank scale ~ 10 19 GeV. The great 
attractiveness of the see-saw model lies in the fact that, through the relation ( |2.73|) , it 
gives an explanation of the smallness of neutrino masses with respect to the masses of 
other fundamental fermions.f] 

In the case of three generations the see-saw mechanism leads to a spectrum of masses 
of Majorana particles with three light neutrino masses and three very heavy masses 
Mk (k = 1,2,3) of the order of the scale of violation of the lepton numbers. This is 
realized if the mass matrix (|2T32] ) has M L = 0, 

(M D ) T 
M R 



M 



D+M 







(2.77) 



and M R is such that all its eigenvalues are much bigger than the elements of M D . In 
this case the mass matrix is block-diagonalized (up to corrections of order (M i? )~ 1 M D ) 
by the unitary transformation 



W T M D+M W 



Mi ight 






-^heavy 



(2.78) 



with 



W 



( 



\ 



i(M D ) t (M fl (M 



R\U-i 



M 



-(M 



R\-l 



M 



1 



(M D ) f (M i? ) t 
flw M D (M D )t(M^)t 



-i 



/ 



The matrices for the light and heavy masses are given by | 112 . 1X3 
M Ught ~ -(M D ) T (M R )~ X M D , M heavy ~ M R 



(2.79) 



(2.80) 



The mass eigenvalues of the light neutrinos are determined by the specific form of M D 
and M R . Note that in left-right symmetric models and in SO(10) models the matrix M L 
in the big matrix M D+M ( |2.32 ) can be important (see, e.g., Ref. ||114 |). 

Two simple possibilities are discussed in the literature (see Refs. [|115| , |116| ): 



'Note that the see-saw mechanism can also be realized for Dirac neutrinos 1 1 1 , 111| 
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1. If M R = A4 I, where I is the identity matrix, one obtains the quadratic see-saw, 



(M D ) T M D , . 

M light ^-^-^ , (2.81) 



and the light neutrino masses are given by 



f\2 



where m{ is the mass of a quark or a charged lepton of the k th generation. In this case 

J\2 . f^f\2 . f^f\2 



the neutrino masses m k scale as the squares of the masses m{: 



mi '■ m 2 '■ m 3 — ( m {) '■ ( m i) '■ ( m D • (2.83) 

re A4b chara 
see-saw (see, for example, Ref. \\L17\\ ) 



2. If M R = -0- Md, where Air> characterizes the scale of Md, one obtains the linear 



M 

and the light neutrino masses are given by 

-Md 



(2-84) 



m{ (k = 1, 2, 3) . (2.85) 



M 

In this case the neutrino masses scale as the masses m{: 

TOi : m 2 : = m{ : vrJ 2 : . (2.86) 

Let us stress that in any case the see-saw mechanism implies the hierarchical relation 

rrii rn 2 m 3 (2.87) 
for the three light Majorana neutrino masses. 

2.6 Effective Lagrangians 

In the Standard Model without right-handed singlet neutrino fields there are no renor- 
malizable interactions that give masses to the neutrinos after the spontaneous breaking 
of the SU(2) L x U(l) y symmetry with the Higgs doublet mechanism. However, there is 
a general belief that the Standard Model is the low-energy manifestation of a more com- 
plete theory ||118| , |119|| (for reviews see Refs. |pTj|120|| ). The effect of this new theory is to 



induce in the Lagrangian of the Standard Model non-renormalizable interactions which 
preserve the SU(2) L x U(l) y symmetry above the electroweak symmetry breaking scale, 
but violate the conservation of lepton and baryon numbers (see Ref. ||121[| and references 



therein). These non-renormalizable interactions are operators of dimension d > 4 and 
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must be multiplied by coupling constants that have dimension Ai 4 ~ d , where M. is a mass 
scale characteristic of the new theory. It is clear that the dominant effects at low energies 
are produced by the operators with lowest dimension. 

In the Standard Model the lepton number non-conserving operator with minimum 
dimension that can generate a neutrino mass is£] ||122[ - |127|| 



77 ^ ( L l C l(J 2vL l ,) (<f) T a 2 acj))+ h.c. 



M 2 



(2.88) 



where g is a 3 x 3 matrix of coupling constants, a are the Pauli matrices, Li are the 
standard leptonic doublets and is the standard Higgs doublet: 



L, 



h 



(£ = e,fi,r) 



'■P 



(2.89) 



The effective operator in Eq. fl2.88|) has dimension five and the coupling constant is pro- 
portional to M.^ 1 . 

When SU(2) L x U(l) y is broken by the vacuum expectation value v/y2 of ip°, the 
effective interaction (|2.88j) generates the Majorana mass term 



(2.90) 



£M = \ ~M ^ 9u ' U * L C lvi ' L + h ' C ' = ~\~~M^ 9u ' V(JL + h ' C ' 



2 M 



£,£' £,£' 

where v ~ 246 GeV. In this mass term there is a suppression factor v/M. which is 
responsible for the smallness of the neutrino masses [122,123]. 



2.7 Maximal mixing 

The expression ( |2.64 ) for the mixing angle d implies that the mixing is maximal, i.e., 
d = 7r/4, if rriR = m^. In this case, assuming that > [mr>| , the Majorana neutrino 
masses are given by 

m \,2 = m L T m D • (2-91) 
The fields u L and {vr) c are connected with the fields V\ and v% by the relations 

VL = {VlL + V2L) , = {~V\L + V2L) ■ (2.92) 

The massive Majorana fields are given by 

v lt2 = i= {[u L + [y L y\ =F [u R + (v R ) c }} . (2.93) 



7 This operator can be written also as 

9U> {Lj o 2 <£) C- 1 (<f> T a 2 L v ) + h.c. 
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If ttir = mi = and m-Q > 0, the mass term ( 2.57|) is simply a Dirac mass term. 
Applying Eqs. (|2.64j) and (|2.66|) , we obtain 



7T 

4 



m l 2 = =Fm D • 



(2.94) 



Using Eq.([2~7TD for the mixing matrix U, one can see that the fields vl and (vr) c are 
connected to the Majorana fields v\ and z/ 2 by the relations 



u L = -j= (-iviL + v 2 l) , 



V2 



(ivi L + v 2 l) , 



(2.95) 



where v\$, are Majorana fields with the same mass m\ = m-i = mo and with CP parities 
rji F = —i and r]^ = +i (see Eq.( |2.53|) ). Eq. (|2.93| ) implies for a Dirac field 

1 



(2.96) 



Thus we arrive at the well-known result that a Dirac field can always be represented as 
an equal mixture of two Majorana fields with the same mass and opposite CP parities. 
One can see this result also directly: 



v — — - 



V — V 



V2\ \/2 



+ 



1 z/ + z/ c 



y/2 V2 



where 



v — v 



V + v c 



(2.97) 



(2.98) 



V2 ' Z V2 
are Majorana fields. 

Finally, there is the possibility that \mn\ <C m-£> but at least one of the param- 
eters w,L t R is non-zero. In this case Eqs. fl2.95 ) and (|2.96|) are approximately valid and 
v>i } 2 are two Majorana neutrinos with opposite CP parities and almost degenerate masses 
given by 

mi 2 ~ m D =f - (mi + m R ) . (2.99) 



The field v is called in this case a pseudo-Dirac neutrino field p2|,|65|,p^8Hl30||. 



3 Neutrino oscillations in vacuum 
3.1 The general formalism 

If there is neutrino mixing, the left-handed components of the neutrino fields v a ^ (a = 
e, fi, t, s\, S2, ■ ■ •) are unitary linear combinations of the left-handed components of the n 
(Dirac or Majorana) neutrino fields (k = 1, . . . , n) with masses m^: 

n 

VgL = ^ Ugh VhL ■ (3.1) 

k=l 
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The number n of massive neutrinos is 3 for the Dirac mass term discussed in Section 1270 



and for the Majorana mass term discussed in Section |2.3| , in which cases there are only 
the three active flavour neutrinos. The number n of massive neutrinos is more than three 



in the case of a Dirac-Majorana mass term discussed in Section |2.2| with a mixing of 
both, active and sterile neutrinos. In general, the number of light massive neutrinos can 
be more than three. We enumerate the neutrino masses in such a way that 



mi < ni2 < rriz < . . . < m n 



(3.2) 



In this section we will consider in detail the phenomenon of neutrino oscillations in vacuum 
which is implied by the mixing relation Q3.1| ). 

If all neutrino mass differences are small, a state of a flavour neutrino v a produced 
in a weak process (as the tc + — > decay, nuclear beta-decays, etc.) with momentum 
p ^> ink is described by the coherent superposition of mass eigenstates (for a discussion 
of the quantum mechanical problems of neutrino oscillations see Refs. [|BT|, |62| , p~31| |1 46|| ) 



(3.3) 



fc=i 



Here is the state of a neutrino with negative helicity, mass and energy 



Ek 



p 2 + ml ~ p + 



2p 



(3.4) 



Let us assume that at the production point and at time t = the state of a neutrino 
is described by Eq. (|3.3|). According to the Schrodinger equation the mass eigenstates |z/ fc ) 
evolve in time with the phase factors exp(-iEkt) and at the time t at the detection point 
we have 



(3.5) 



k=l 



Neutrinos are detected by observing weak interaction processes. Expanding the state 
( ^.5[ ) in the basis of flavour neutrino states we obtain 



(t) W P ) 



where 



(3.6) 



(3.7) 



fc=i 



is the amplitude of v a — > v$ transitions at the time t at a distance L ~ t. Consequently, 
the probability of this transition is given by 



\A, 



e~ lEkt K k 



k=l 



(3.8) 
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This formula has a very simple interpretation. U* k is the amplitude to find the neutrino 
mass eigenstate \u k ) with energy E k in the state of the flavour neutrino \v a ), the factor 
exp(— iEkt) gives the time evolution of the mass eigenstate and, finally, the term Upk 
gives the amplitude to find the flavour neutrino state \vp) in the mass eigenstate \v k )- We 
want to remark that weak interaction processes responsible for neutrino production and 
detection involve active neutrinos. Therefore, strictly speaking, in the derivation of the 
Eqs. (|3.7| ) and fl3.8|) the indices a and (3 take only active flavour indices. However, transi- 
tions into sterile neutrinos can be revealed through neutral current neutrino experiments 
(disappearance of active neutrinos). In this sense, the formulas ( [3.71) and (|3.8| ) have a 
meaning also for transitions between active and sterile states. 

Notice that in order to have a non-negligible active-sterile transition probability the 
sterile fields must have a mixing with the light neutrino mass eigenfields the number 
of which must be more than three. Such a possibility is phenomenologically given by 
the Dirac-Majorana mass term ( |2.21| ), but it is not realized in the simple see-saw scheme 



discussed in Section |2.5| , where the scale of the right-handed Majorana mass term is large. 
However, the see-saw scenario can be modified by additional assumptions to include light 
sterile neutrinos ("singular see-saw" 1 147|| , "universal see-saw ||148| , |149|1 ). 



At this point a remark concerning the unitarity of the mixing matrix U is at order. 
If some of the mass eigenstates are so heavy that they are not produced in the standard 
weak processes then these mass eigenstates will not occur in the flavour state \v a ) Q3.3|) . 
Let us assume that the first n' mass eigenstates are light (n' < n). Consequently, only 
that part of U plays a role in neutrino oscillations where k < n' . In the following we will 
always assume that in the situation described here we can confine ourselves to an n' x n' 
submatrix of U which is unitary to a good approximation (see, e.g., Ref. [ 150(| ). This 



is realized in the see-saw mechanism with a sufficiently large right-handed scale. In the 
further discussion we will drop the distinction between n and n'. 

From the relation (|3.1|) it follows that the state describing a flavour antineutrino v a 
is given by 



\v a ) = y^U ak \v k ) . (3.9) 



k=l 



Thus the amplitude of v a — > vp transitions is given by 

n 

A^(t) = Y,U; k e-^U ak . (3.10) 



k=l 



Notice that the amplitude for antineutrino transitions differs from the corresponding 
amplitude ( |3.7| ) for neutrinos only by the exchange U — > U*. 
Using the unitarity relation 

n 

J2UpkU* ak = 5 a p, (3.11) 

k=l 
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Experiment 


L (m) 


E (MeV) 


Am 2 (eV 2 ) 


Reactor SBL 


10 2 


1 


io- 2 


Reactor LBL 


10 3 


1 


KT 3 


Accelerator SBL 


10 3 


10 3 


1 


Accelerator LBL 


10 6 


10 3 


10~ 3 


Atmospheric 


10 7 


10 3 


IO" 4 


Solar 


10 11 


1 


io- 11 



Tabic 3.1. Order of magnitude estimates of the values of Am 2 which can be probed in reactor short- 
baseline (SBL) and long-baseline (LBL), accelerator SBL and LBL, atmospheric and solar neutrino os- 
cillation experiments. Note, however, that the energies and distances of the various types of experiments 
can vary in a wide range and only some representative values are given in this table. 



the probability ( |3.8|) can be written as 



p 



+ u pk u* k 



k=2 



exp 



Am 2 ^ 
IE 



(3.12) 



where AmL = m^—m 2 and the ultrarelativistic approximation ( |3.4| ) has been used. Thus 
the probability of v a — > v@ transitions depends on the elements of the mixing matrix, on 
n — 1 independent mass-squared differences and on the parameter L/E, whose range is 
determined by the experimental setup. 

If there is no mixing (U = I) or/and /\m 2 kl L/ E <C 1 for all k — 2, . . . , n, there are no 
transitions (P Ua -, U/3 = 5 a p). Neutrino transitions can be observed only if neutrino mixing 
takes place and at least oneQ Am 2 satisfies the condition 



Am 2 > j . 



(3.13) 



In this inequality, Am 2 is the neutrino mass-squared in eV 2 , L is the distance between 
neutrino source and detector in m (km) and E is neutrino energy in MeV (GeV). Thus, 
the larger the value of the parameter L/E, the smaller are values of Am 2 which can 
be probed in an experiment. The inequality (|3.13j ) allows to estimate (for large mixing 
angles) the sensitivity to the parameter Am 2 of different types of neutrino oscillation 
experiments. These estimates are presented in Table |3.1|. Let us stress that this table 



gives only a very rough idea of the sensitivity to Am . For example, in the LSND 
short-baseline (SBL) experiment antineutrino energies are between 20 and 60 MeV, the 
distance is approximately 30 m and the minimal value of Am 2 probed in this experiment 
is 4 x IO" 2 eV 2 . 



3 We use the symbol Am 2 for a generic mass-squared difference. Hence Am 2 can be any of the Am^ 's. 
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The probability ( |3.8| ) and the corresponding one for antineutrinos are invariant under 
the phase transformation 



U, 



ak 



-up a 



U ak e 



(3.14) 



Therefore, it is clear that the probabilities of v a — > vp and v a — > vp transitions do 
not depend |28, |29|, pT06[1 on the Majorana CP- violating phases discussed at the end of 
Section and it is not possible to distinguish the Dirac and Majorana cases by the 
observation of neutrino oscillations. 

Comparing the expressions (3.7) and (|3.10| ) for the transition amplitudes of neutrinos 



and antineutrinos we see that A Voi ^ Vg {t) = Aj? g -+v a (t). Therefore, for the transition 
probabilities we have 

Pv a -^Vp — Pv^V a ■ (3.15) 

This relation is a consequence of CPT invariance inherent in any local field theory ||151| . 
From the equality (|3.15|) it follows that the neutrino and antineutrino survival probabil- 
ities are equal: 

P„ = P^ a • (3.16) 



On the other hand, the transition probabilities of neutrinos and antineutrinos in general 
are different for a ^ (3. They are equal only if there is CP invariance in the lepton sector. 
In fact, in the case of massive Dirac neutrinos the phases of the neutrino fields and of the 
charged lepton fields can be chosen in such a way that the mixing matrix U is real. In 
the case of Majorana neutrinos, CP invariance implies that (see Eq. (|2.50) )) 



TJ* 

U ak 



U a k Pk 



(3.17) 



with pk = —iVk F = where t]^ F is the CP parity of the Majorana neutrinos with mass 
mfc (see Eq.( p.53[ )). Thus, in the Dirac as well as in the Majorana case, CP invariance 
implies that 



P rj 



(3.18) 



3.2 Oscillations in the two- neutrino case 

The results of neutrino oscillation experiments are usually analysed under the simplest 
assumption of oscillations between two neutrino types. In this case, for the transition 
probability ( |3.12j ) we get 



Pu 



8aj3 + U '/32U* 2 



exp 



.Am 2 L 



(3.19) 



where Am 2 



ml 



ml and a, [3 are e, /i, or fi, r, etc. . . Thus in the simplest case of 



transitions between two neutrino types the probability is determined only by the elements 
of U which connect flavour neutrinos with u 2 (or Ui). It is obvious that phases drop out 
in the expression ( |3.19| ). This is an illustration of Eq.( |3.14j ) and of the fact that no 
information about CP violation can be obtained in the case of transitions between only 
two neutrino types. If we put U a 2 = sin^ then we have U$i — cost? and the transition 
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and survival probabilities are given by the following standard expressions valid for both, 
neutrinos and antineutrinos: 



R-) (-) 



sin 2-d ( 1 — cos 



Am 2 L 



2E 



H-) (-) 



P(-) <-) • 



The transition probability ( |3.20| ) can be written in the form 



(3.20) 
(3.21) 



(-) 



- sin 2 2$ I 1 - cos 2.53 



Am 2 L 
E 



(a ? (3) . 



(3.22) 



where L is the source - detector distance expressed in m (km), E is the neutrino energy 
in MeV (GeV) and Am 2 is the neutrino mass-squared difference in eV 2 . Thus, the 
transition probability is a periodic function of L/E. This phenomenon is called neutrino 
oscillations. The amplitude of oscillations is sin 2 2i9 and the oscillation length is given by 



Am 2 ' Am 2 (eV 2 ) 
The condition ( |3.13| ) can be rewritten in the form 

L osc < L . 



m . 



(3.23) 



(3.24) 



Therefore, neutrino oscillations can be observed if the oscillation length is not much larger 
than the source - detector distance L. 

The oscillatory behaviour of the transition probability ( E>.2U| ) with sin 2 2$ = 1 is shown 
in Fig. |3.1| , where we have plotted it as a function of Am 2 L/AirE = L/L osc . The grey line 
represents the transition probability ( |3.20| ), whereas, in order to demonstrate the effect of 
energy averaging, the black line represents the transition probability ( |3.20D averaged over 
a Gaussian energy distribution with mean value E and standard deviation a = E/10. The 
averaged probability is the measurable quantity in neutrino oscillation experiments. One 
can see that the averaging over the energy spectrum practically reduces the probability 
to the constant 1 - sin 2 2$/ 2 for L > L osc . 

The expressions ( |3.20|) and ( |3.21|) are usually employed in analyses of the data of 
neutrino oscillation experiments. In many SBL experiments with neutrinos from reactors 
and accelerators, no indication in favour of neutrino oscillations was found. The data 
of these experiments give an upper bound for the transition probability which implies 
an excluded region in the space of the parameters Am 2 and sin 2 2i9. A typical exclusion 



plot is presented in Fig. |3.2| [152]. This plot shows the exclusion curves in the 
channel obtained in the CDHS |53j, FNAL E531 pi 



CHARM II PSj , CCFR jTSEfl , 



CHORUS [157] and NOMAD [158] experiments. The excluded region lies on the right of 
the curves. 



The two most stringent exclusion curves in Fig. [3.2| have been obtained in the CHO- 
RUS [f57] and NOMAD ||158|] experiments, which are operating at CERN using the 
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Figure 3.1. Transition probability for sin 2 2-!? = 1 as a function of Am 2 L/4:irE = L/L osc , where L osc 
is the oscillation length. The grey line represents the transition probability ( |3.20D and the black line 
represents the same transition probability averaged over a Gaussian energy spectrum with mean value 
E and standard deviation a = E/10. 



neutrino beam from the SPS (with an average energy of about 30 GeV). 800 kg of emul- 
sions are used in the CHORUS experiment as target. The production and decay of r's 
in the emulsion is searched for. In the NOMAD experiment a magnetic detector is used 
and the production of r's is identified with kinematical criteria. 

Figure |0| shows the exclusion curves obtained in the CHOOZ [|159| , Gosgen | |160| |, 
Krasnoyarsk [ |161|| and Bugey ||162|| reactor v e — > v e experiments. The region allowed by 
the results of the Kamiokande atmospheric neutrino experiment (see Section |5.1| ) is 
also depicted in this figure. The Gosgen, Krasnoyarsk and Bugey experiments are SBL 
reactor experiments, whereas the recent CHOOZ experiment is the first LBL reactor 
experiment. In this experiment the detector (5 tons of liquid scintillator loaded with Gd) 
is at the distance of about 1 km from the CHOOZ power station, which has two water 
reactors with a total thermal power of 8.5 GW. The antineutrinos are detected through 
the observation of the reaction 

u e +p^e + + n (3.25) 

(the photons from annihilation of the positron and the delayed photons from the capture 
of neutron by Gd are detected). No indications in favour of neutrino oscillations were 
found in the CHOOZ experiment. The ratio R of the numbers of measured antineutrino 
events and of expected antineutrino events in the case of absence of neutrino oscillations 
in the CHOOZ detector is 



159 



R = 0.98 ±0.04 ±0.04. 



(3.26) 



Since the average value of L/E in the CHOOZ experiment is approximately 300 ((E) ~ 
3MeV, L ~ 1000 m), in this experiment it is possible to probe the value of the relevant 
neutrino mass-squared difference down to 10 _3 eV 2 . 
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|3.2[ Exclusion curves (90% CL) in the 
v T channel obtained in the CDHS [ 153 1 , 



FNAL E53TJ15J], CHARM II [155 1, CCFR |56| 
CHORUS |157| and NOMAD pf experiments. 



Figure |3.3| . Exclusion curves (90% CL) i n th e 
v e — ► v e cha nnel obtained in th e CHOOZ [1591, 
Gosgen [ 160 1 , Krasnoyarsk [161] and Bugey [162 
experiments. The shadowed region is allowed by 
the results of the Kamiokande atmospheric neu- 
trino experiment |49J. 



With the help of the expression ( |3.22j ), it is possible to understand qualitatively the 
general features of exclusion curves. In the region of large Am 2 such that the oscillation 
length is much smaller than the source - detector distance L, the cosine in the expression 
( p.22| ) oscillates very rapidly as a function of the neutrino energy E. Since in practice 



all neutrino beams have an energy spectrum and the neutrino sources and detectors are 
extended in space, only the average transition probability 

(Pi-) ( -,> = ^sin 2 2tf (3.27) 

can be determined in the region of large Am 2 . The average probability is independent 
from Am 2 and, therefore, from an experimental upper bound ( _))o on (i(_ 3 ( _ } ) one 

obtains the vertical-line part of the exclusion curve. 
At 

Am 2 ~ - — ^- , (3.28) 
2 1.27 (L) ' V ; 

where (E) is the average energy and (L) is the average distance, the parameter sin 2 2$ 
has the minimal value 

(sin 2 2tf) min = ( _)) (3.29) 
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on the boundary curve. Note that we are using here and in the rest of the section the 
same units as in Eq. (|3.22| ). 

Typically, the upper bound (if_) ( _))o is much less than one. Then, in the region where 

sin 2 2$ is large the expression ( |3.22|) for the transition probability can be approximated 
by 

Am 2 L x 2 



if_, (_, ~ sin 2 2$ 1.27 — — (3.30) 

Va^Vfj \ El 



and in this region the boundary curve in the exclusion plot is given by 



'{Pi-) C-)>0 

Am 2 ~ / . (3.31) 

1.27v/sin 2 2tf (L 2 )(E~ 2 ) 

Therefore, this part of the exclusion curve is a straight line in the log sin 2 2$ - log Am 2 
plot as can be seen from Fig. [3.2| . From Eq. (|3.31|) it follows that the minimal value of 
the parameter Am 2 that can be probed by an experiment is 



'(P(-) c-))o 

Am 2 „ v ""^ = . (3.32) 
1.27 ^(L 2 )(E~ 2 ) 



It corresponds to sin 2 2$ = 1. 



4 Neutrino oscillations and transitions in matter 
4.1 The effective Hamiltonian for neutrinos in matter 



It has been pointed out by Wolfenstein [|163j| and by Mikheyev and Smirnov [164] that 



the neutrino oscillation pattern in vacuum can get significantly modified by the passage 
of neutrinos through matter because of the effect of coherent forward scattering. This 
effect can be described by an effective Hamiltonian. Starting with neutrino oscillations 
in vacuum, one can easily check that the transition probability (|3.12| ) can be obtained by 



considering the evolution of the state vector ip of the neutrino types with the "Schrodinger 
equation" 

= H ™^w = 2^ u ^^ x ) > c 4 - 1 ) 

where if ™ c is the effective Hamiltonian for neutrino oscillations in vacuum with U being 
the mixing matrix, E the neutrino energy and rh the diagonal neutrino mass matrix. The 
presence of matter will give a correction to Eq. fl4.ip . Note that the variable in Eq.( [4.1[ ) 
is not time but space and the components of ip(x) are given by the amplitudes denoted 
by a a (x) for the neutrino types a = e,/i,r, s. The corresponding effective Hamiltonian 
for antineutrinos H™ c is obtained from if™ c by the exchange U — > U*. After the seminal 
work of Wolfenstein about neutrino oscillations in matter (see also Ref. [165]) and the 



discovery of the importance of if ™ at , the analogue of ii™ c in matter, for the solar neutrino 
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or 



problem [|164j| many rederivations of the effective matter Hamiltonian ||f 66| , |f 67|| have 
been presented using a coupled system of Dirac equations |168| , |169|| , field theory ||170 
studying coherent forward scattering in more detail [171|. 

We want to sketch a derivation using the Dirac equation and following Ref. 



169| . The 



starting point in most derivations is the expectation value of the currents for isotropic 
non-relativistic matter given by ||163| , |164| , |172j 



matter 



■NfS. 



(4.2) 



where Nf is the number density of the particles represented by the field /. Note that the 
75 term does not contribute for non-aligned spins ||173|| . Eq.( |4.2|) is a good approximation 
even for electrons in the core of the sun with a temperature T ~ 16 x 10 6 K [ JI74j , |175| and 
thus an average velocity of the electrons of around 10% of the velocity of light and correc- 
tions of order kT/m e ~ 2.6 x 1CT 3 . Starting from the weak interaction Lagrangians ( |2.1|) 
and ( |2.2j ) one gets for low-energy neutrino interactions of flavour £ with the background 
matter 

(rr 1 ^ s " ' e 2sin 2 ^Q/), (4.3) 



£mat J 



v\{\ 



7s)^^iV/(% + T 3fL 

where Gp is the Fermi coupling constant, 9w the weak mixing angle, T 3 f L the eigenvalue of 
the field /l of the third component of the weak isospin (T 3 f R = in the Standard Model) 
and Qf is the charge of /. In the matter Lagrangian (|4.3j ) the charged current interaction 
is represented by the Kronecker symbol Sgf saying that for neutrinos of flavour I the 
charged current only contributes when background matter containing charged leptons of 
the same flavour is present. Concentrating now on realistic matter with electrons, protons 
and neutrons which is electrically neutral, i.e., N t 
have T 3eL = -T< 



3pl 



3n L 



N p , and taking into account that we 
— 1/2 and Q e = —Q p = —1, Q n = for electrons, protons 



and neutrons, respect ively,[] we get an effective Hamiltonian 



H, 



D 



-ia 



■ V + (3(MP L + M^P R ) + v^GVdiag (iV e - ijV n , ~N n , ~N r )j P L 



(4.4) 



Pl,r = (1T 7s)/2, otj 



acting on the vector \l/(t, x) of the flavour neutrino wave functions, where we have defined 

7°7 : ' and (3 = 7 . M is the non-diagonal mass matrix. The 
is called the matter potential term. Several remarks concerning 
, which is the point of departure for deriving the effective matter 
T68| , [169[ 



last term in Eq. 
the Hamiltonian 
Hamiltonian in Refs. 



are at order: 



1. H„ has spinor and flavour indices, it leads thus to a system of Dirac equations coupled 
via the mass term. 

2. The neutral current contributions of electrons and protons cancel for realistic matter 
because of opposite T 3 f L quantum numbers and electric charges (see discussion above 
before Eq.(fO)). 



For protons and neutrons this is ensured by CVC (conserved vector current) for the T3 part and 
conservation of the electromagnetic current. 
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3. The Hamiltonian ( |4.4| ) is valid for Dirac neutrinos. For antineutrinos, M is replaced by 
M T and the matter potential term of has the opposite sign and the projector Pr 
instead of Pl- This is a consequence of proceeding along the same lines as before but 
using the charge-conjugate fields to describe antineutrinos {ptipPhVi = —v^pPRvf). 

4. In the Majorana case, in order to obtain the Hamiltonian H„ the antineutrino matter 
potential has to be added to and the neutrino field vector \1/ is subject to the 
Majorana condition \1/ = C^'fy* . This condition is compatible with the time evolution 
governed by H^ 1 [|169|| . 



5. Sterile neutrinos can easily be incorporated into HP, and H^ 1 by simply adding the 

entry to the matter potential term. Note that the probability of active - active 
transitions does not depend on the neutron density N n , which is, however, important 
for active - sterile transitions (see Table |4.1| ). 

The essence of deriving an effective Hamiltonian for neutrino oscillations in matter 
is to get rid of the spinor indices in and to obtain an equation involving only the 
indices a for the different neutrino types as in the vacuum case ( |4.1|) . To this end, let 
us now assume that the neutrino propagation proceeds along the x 3 = x axis, that the 
neutrino momentum p > corresponding to propagation in vacuum is much larger than 
the matter potentials and that the neutrinos are ultrarelativistic. Thus we consider a 
one-dimensional problem from now on. Defining a wave function via 

V(t,x) = e ipx <P(x,t) with t ^=( P a 3 + H^)(f ) , (4.5) 

this wave function changes little over distances of the order of the de Broglie wave length 
of the neutrino. With respect to 0:3 we can decompose the Hamiltonian ( f4.4|) into = 
H even + H odd which are the parts commuting or anticommuting with a 3 . H odd consists 
solely of the mass term and H even contains the rest of . With p being a large parameter 
we can perform a Foldy - Wouthuysen transformation where we truncate the series at 
1/p leading to the Hamiltonian [|169| , |176| 



Hvw,v = P a 3 + H cven + —a 3 H^ dd 

2p 



d 



j + V%Gp diag ^N e - ^N n , ~N n , ~Nnj P L 



+ ^-a 3 (M^MP L + MM^P B ), (4.6) 
Zp 

such that all terms in Hp W v commute with a 3 and, therefore, all matrices in Hp W v with 
spinor indices can be diagonalized at the same time in order to separate positive and 
negative energy states. Denoting this unitary diagonalization matrix by Uq we can achieve 
U a 3 Ul = diag(l, 1, -1, -1), Uo^pul = diag(0, 1, 1, 0) and U ^U j = diag(l, 0, 0, 1). 
Note that there are no transitions between left and right-handed states through H^ w u . 
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For neutrinos with magnetic moments (electric dipole moments) the method described 
above would also yield the appropriate left-right transitions ||169|| . 

Going back to the Foldy - Wouthuysen transform of the wave function \l/ instead of <p 
removes the p from the Hamiltonian (|4.6|) . The final step in deriving the effective matter 
Hamiltonian consists of considering stationary states and splitting off the plane wave part 
by 



(4.7) 



Then the wave function ip(x), taken for positive energies and left-handed neutrinos, fulfills 
the first order differential equation ( |4.1|) with H™ c replaced by [ |163| 



\ (M^M + 2V2 EG F dia^ 
IE \ 



N e - ^N n , --N n , --N n 



correct up to order 1/E. This effective Hamiltonian is valid for left-handed Dirac neutri- 
nos or Majorana neutrinos, whereas in the case of right-handed (Dirac) antineutrinos or 
right-handed Majorana neutrinos we have 



(4.9) 



#p mat = ^ (M T M* - 2V2 EG F diag (n c - ~iV n , -~N n , J , 

because for (Dirac) antineutrinos one has to replace M by M T in the Hamiltonian ( f4.6|) 
and for Majorana neutrinos M = M T holds (see also points || and Rafter Eq.( j4.4|) ). With 



M = U R mUl and U L = U 



(4.10) 



only the left-handed mixing matrix U enters and connection between Eqs. (|4.8| ) and ( f4. 1|) 
is made. Thus in the final results ( f4.8| ) and ( f4.9| ) there is no distinction between the 
Dirac and Majorana case for ultrarelativistic neutrinos. This is an illustration of the 
fact that with neutrino oscillation experiments the Dirac or Majorana nature cannot be 
distinguished |98| , |177fl . 

The Hamiltonians ( |4.8|) and ( [4.9D have been used to investigate neutrino oscillations 
in the sun, in the earth and in supernovae. In the following we will only be concerned 
with the first two subjects. Application limits of the neutrino evolution equations in 
matter have been discussed in Ref. ||178| |. Elastic and inelastic neutrino scattering in- 
troduces quantum damping into the evolution equations which is proportional to the 
neutrino interaction rate [|179|| . In the sun and the earth this effect is negligible, in par- 



ticular, for low energies, whereas in the early universe it is of crucial importance 1 179 



Density fluctuations have been found to influence considerably neutrino propagation in 



the sun [ 180 1 , however, more realistic considerations with helioseismic waves as density 
fluctuations show no observable effect on the solar neutrino problem discussed in terms 



of neutrino oscillations [181]. Solar neutrinos are also influenced by their passage through 



the earth [182 
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A 


N e N e -\N n -\N n 


2\f2EG F 



Table 4.1. The list of matter densities relevant for two-neutrino oscillations (for the definition of A see 
Eq.(ETf). 



4.2 The two-neutrino case and adiabatic transitions 



Let us now concentrate on left-handed neutrinos and specify the effective Hamiltonian 
( [4.8D to two neutrino types. Thus for two-neutrino oscillations in matter with the defini- 
tions 



N(v a ) = 5 ae N e - ^N n {a = e, fi, r) , N{v s ) = 



and 



A = 2V2G F E (JV( 



N{v p )) 



(4.11) 
(4.12) 



the differential equation 
. d 



dx 



dp 



1 

AE 
+ 



mat 



ap 



ml + m 2 2 + 2V2G F {N{v a ) + N{v p )) 



A — Am 2 cos 2i9 Am 2 sin 2^ 
Am 2 sin 2t9 —A + Am 2 cos 2$ 



1 
1 



a/3 



(4.13) 



has to be studied ||1 63| . |1 64| , |1 72j , |1 78| , |1 83| . |1 84|| , where a a , ap are the amplitudes for the 
neutrino types a, (3 (a, (3 = e,/i, r, s), respectively. The evolution equation for antineu- 
trinos is obtained from Eq. (|4.13|) by A — > —A. In the two-neutrino case there is one 



Am 2 



m 



2 — m\ and the mixing matrix 

U = 



cos $ sin $ 
sin $ cos $ 



(4.14) 



is real without loss of generality. The list of all possible matter densities which determine 
A and occur in the different oscillation channels is given in Table |4.1| . Evidently, z/ M <-> v T 
oscillations proceed as in vacuum because A = and the term proportional to the unit 
matrix in Eq.( (4.13j ) has no effect on transitions. As mentioned before, in the rest of this 



section we will have in mind that Eq.( 4.13 ) is applied to neutrino propagation in the sun 
and the earth. 

Let us first define the eigenf unctions of the effective Hamiltonian as 



i^mj = Ejip mj with Vmi 



cosi? m 
• sin d m 



m2 



sin dm 

cos -d m 



(4.15) 
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where the matter angle $ m and the energy eigenvalues Ej are functions of x. In the 
limit of vanishing matter density $ m is identical with the vacuum angle d ( 14.141) . The 
eigenvalues Ej and the matter angle are given by 

£1,2 = -^{m 2 1 + m 2 2 + 2V2G F (N(u a )+N(u (3 )) 

=Fa/ (A — Am 2 cos 2$) 2 + (Am 2 sin 2tf) 2 } (4.16) 



and 



tan 2$ 

tan2tf m = j , (4.17) 



Am 2 cos 2$ 

respectively. The general solution of the differential equation (|4 .13 ) can be represented 



as 

2 



ip(x) = aj{x)ip m j{x) exp ( — i / dx'Ej(x') j (4-18) 

j = l V ^XO / 

with the coefficients fulfilling 

± ( ai \ = ( exp (i J: q dx'AE(x')) \ / fll 

dx V «2 / I ^exp (-if* dx'AE{x')) I V a 2 



(4.19) 
and AE = E 2 - E x . 

The adiabatic solution is defined as an approximate solution where dj ~ (j = 1, 2).Q 
Postponing the discussion of the question under which condition adiabaticity is fulfilled, 
we consider temporarily the case of an arbitrary number n of neutrino flavours or types 
and define the mixing matrix in matter U m (x) via 

U m (x) ] H^U m {x) = diag(#i(a;), . . . , E n {x)) with U m {x) = (Vw, • • • , VwO , (4.20) 

where the eigenvectors ip m j generalize Eq.( |4.15 ). One readily obtains the generalization 
of the vacuum oscillation amplitude Eq. fl3.12| ) for neutrino oscillations in matter if the 
evolution of the neutrino wave function is adiabatic: 

= l>m(*i)»«P (~i \Sj + £ dx'Ejix')^ U* m (x ) aj . (4.21) 

In this formula neutrino production and detection happen at x and x±, respectively. The 
adiabatic phases 8j, which are defined by 

Sj = -i ! dx'^ m j{x')H m j{x') , (4.22) 



10 



We use the dot above a symbol as abbreviation for 4- 
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are necessary for the correct evolution of ip(x) in the adiabatic limit. Usually, they can be 
absorbed into the eigenvectors ipmj [|185|| , except for special matter densities, where these 
phases acquire a topological meaning, together with the presence of CP violation [186]. 
Note that for real vectors ip m j one has ipmjipmj = and therefore 5j = 0.0 This is so in 
the two-dimensional case Eq. (|4.15 ). 

Evaluating Eq. fl4.21|) for two neutrino types and assuming that an averaging over 
neutrino energies takes place such that 



exp —i 



x 1 



dx'AE(x r 







(4.23) 



it is easy to show that the averaged two-neutrino survival probability can be written in 
the adiabatic 



P Va -+u a = ^ + cos2$ m (x )cos2$ m (£i)) . 



(4.24) 



Let us now derive a condition for adiabaticity in the two-neutrino case. The evolution 
of the neutrino state in matter is adiabatic if the right-hand side of Eq. (|4.19 ) can be 
neglected. A formal solution of Eq.( 4.19 ) is given by 



ai(x) 
a 2 (x) 



J2 / dx'J(x') / dx"J(x") 
k=o 



■I'D 



da;W J(xW) ( ai \ X i ) , (4.25) 
V a 2 (x ) I 



where J denotes the matrix on the right-hand side of Eq. ( |4.19| ) . Defining the variable 



1 

a = - 
2 



dx'AE(x') 



and the adiabaticity parameter (as a function of x) 

AE 



rg^rgg] 



7W 



2\A 



(4.26) 



(4.27) 



one can determine an upper bound to the typical integral which occurs in all the terms 
of the sum ( [L25| ) on the right end. In an interval [yo, yi] where 7 is monotonous there is 
a value y c such that [|190|| 



± 



m 



dy $ m cos 2a 



!J'< 



da — cos 2a 



7(l/o) 



da cos 2a + 



7(2/i 



yi 



da cos 2a 
(4.28) 

according to the second mean value theorem of integral calculus (yo < y c < yi)- This 
consideration allows to put the bound 



dyd m e 



2ia 



.1/0 



2V2 

< with 7, 

7min 



, uin = min j(y) 
y£[yo,yi] 



(4.29) 



The scalar product in Eq.(4.22) is imaginary for arbitrary vectors ip m j{x) normalized to one Vx, 
therefore, for real vectors it must be zero. 
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Having found this bound, in the rest of the integrals in the terms in Eq.( f4.25|) one can 



simply take the absolute values of the integrands and one is left with integrations of the 
type dx|^ m | = |di? m |. With the refinement that the interval [xq,x] possibly has to be 
divided into several parts labelled by a such that in each part interval 7(2/) is monotonous 
we can define 



and obtain the exact bound 



\a 



[x)-aj{x )\ < 



2V2 ( coshAB sinhA© \ / |ai(x )| 



7 \ sinhA© cosh A© J \ 1 02(^0 



(4.31) 



where the symbol \j denotes the j-th element of the vector on the right-hand side of this 
inequality and A© has a contribution |A$ m | for every interval where d m has a definite 
sign. ( Ad m is the difference between the angle d m in the initial and the final point of such 



an interval.) Thus with the exact bound (|4.31|) we have found the following sufficient 
condition: 

7 ^> 1 =>- adiabatic evolution . (4.32) 

It is fulfilled if on the scale of the oscillation length in matter given by 1/ AE the matter 
angle tf m changes very little, i.e., |i? m | <C AE. Eq.( |4.31| ) allows to get an upper bound 
on the crossing or jumping probability from ip m i to ip m 2- Assuming the initial conditions 
oi(^o) — 1 an d 02(^0) = we get 

2 a/2 

I oa (a;) I < sinhA© (4.33) 

7 

which illustrates once more the importance of 7 for adiabaticity. 



4.3 The resonance 

In the context of the solar neutrino problem the discovery that a resonance in the pas- 
sage of v e through the sun is possible [|164|| gave a major boost to the investigation of the 



propagation of neutrinos in matter. The possibility of a resonance is most easily under- 
stood in the adiabatic approximation by looking at Eqs.( 4.17|) and ( [4.24 ). For different 



neutrino masses we can always label them in such a way that Am 2 > 0. If on the way 
from the creation point xq in matter of v e to a point x\ in vacuum the neutrino passes 
through a point x res where the resonance condition 

A{x rcs ) = Am 2 cos 2d (4.34) 

is fulfilled then $ m (xi) = d has to be between 0° and 45° (for v e the quantity A is 
positive!) whereas according to Eq. ( |4.17|) the matter angle ^ m (x ) is found between 45° 
and 90°. Consequently, the product of cosines in the v e survival probability Q4.24j ) is 

12 Without loss of generality we consider 0° < •& < 90°. 
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negative and the probability to find a v e after the passage through the sun is less than 
1/2. The interesting phenomenon is that this can happen even for relatively small mixing 
angles, provided the matter potential at the production point is large enough such that 
there is a point a; res where the condition Q4.34Q is fulfilled. Indeed, for d close to 0°, 
one can even have $ m (xi) close to 90° and thus P Ue ^ Ue — ( |4.24j ). It has been shown 
that the resonance is also effective in a certain area in the Am 2 -sin 2 2$ plane where the 
evolution of the neutrino state is non-adiabatic. For antineutrinos a resonance is possible 
if 45° < i? < 90°. 

If there is a resonance, then for reasonable matter densities the adiabaticity parameter 
( [4.271 ) is smallest at the resonance and thus adiabaticity is most likely violated there. 

Therefore, considering 7 at the resonance, from Eq.( |4.17|) one easily computes $ m = 

res 

v4 res /(2 sin 2$ Am 2 ) . A suitable measure of adiabaticity in the case of a resonance is thus 
given by P7j , |l88fl 



AE 



7r. 



2\A 



Am 2 sin 2 2$ 
2E cos 2&(\A\/A) t 



(4.35) 



res U--I/ - -/res 

In the solar interior the electron density is maximal in the center with N ( 

6.022 x 10 23 cm" 3 and N n \ ~ ± N e \ 



100 x iV avo where 7V avo 
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max 

This leads to 



2\f2G F E N e 



1.5 x 10" 



1 MeV/ 



(4.36) 



Considering the resonance condition Q4.34 ) one can read off from this equation for which 
Am 2 a resonance is possible. Whereas the density variation of N e and N n in the sun 
is smooth, in the earth due to its layered structure (lithosphere - mantle - core) the 
densities can be approximated as step functions |191| , [L92[| . Furthermore, approximately 
N n ~ A^ e is valid in all layers. In the lithosphere, which is relevant for LBL neutrino 
oscillation experiments, the average mass density p is around 3 g/cm 3 . It increases to 
around 13 g/cm 3 in the core. Thus we can conveniently write 

3^) (tL) ■ (437) 



2V2G F EN e ~ 2.3 x 10~ 4 eV 2 



4.4 Non-adiabatic neutrino oscillations in matter and crossing 
probabilities 

The analogue of the adiabatic formula ( 4.21Q in terms of probabilities is given by 



with 



+U/3 



^2 U m (xi)f3jB(xi, xo) jk U* m {xo) ak 

3,k 



(4.38) 



B{ Xl ,x ) = Ul(x 1 )Pesxp^-iJ" 1 dxH^ t (x)ju m (x ) 

= Pexp|-z^ dx (E -iUl(x)U m (x)^ , (4.39) 
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where E is the diagonal matrix of energy eigenvalues present in Eq.( }4.20D (see Ref. [|183|| ). 
B(x\, xq) is a unitary nxn matrix which is diagonal in the adiabatic limit corresponding 
to -i W m (x)U m (x) ~ diag (5 h . . . , 5 n ) (see Eq.(ggg)). 

Confining ourselves now to two neutrino types and having in mind neutrino production 
in matter and detection in vacuum we use the notation $ m (xo) = $m and $ m (xi) = 

With the crucial assumption that averaging over neutrino energies and the neutrino 
production region all terms other than probabilities can be dropped we get 



(B jk B*, k ,) £ 



djj/dkk' 



\ B jk\ 2 



and \B 



12 



\B- 



21 



Pr, \B 



11 



\B- 



22 



l-Pn 



(4.40) 

where P c is the crossing probability from ip m i to ipm2 and in the second part of Eq. (|4.40|) 
unitarity has been used. Such averaging procedures have been shown to be effective in 
the context of solar neutrinos, e.g., in Refs. ||193| , |194|j . Inserting the relations fl4.40| ) into 



Eq. ( |4.38|) specialized to the survival probability of v e one obtains after some algebra with 
trigonometric functions the Parke formula 



P,, 



1 

2 + 



--PA cos 2tf cos 



(4.41) 



which is the generalization of Eq. (|4.24[) [183, 189fl . This formula was first derived in the 
context of an exact solution of Eq. ([4.13 ) for matter densities linear in x 
also Refs. p7jp6| ). 



188] , [1951 (see 



The probability P c can be estimated with the Landau - Zener method ||196| , |197 
Following the derivation of Landau, the idea is to make an analytic continuation of the 
matter densities and therefore of the effective matter Hamiltonian (|4.8[ ) into the complex 
plane by x — > z where 2 is a complex variable. Then also the neutrino state %p ( [4.18|) is 
analytic. Considering the energy eigenvalues E\^{x) QL16Q as functions of z we find two 
branching points zq and Zq {lmz > 0) defined by the equation 



E^zq) = E 2 (z ) A(z ) = Am'e 



(4.42) 



We consider Ei(z) and i) m \{z) for definiteness and follow its evolution along a curve C(z) 
which starts at a value z in the vicinity of the real axis and goes along a negatively 
oriented loop around zq such that we end at the same point z but now located on the 
second sheet.[^] Along this path E\(z) changes into E 2 (z) ||197|| and ip m i(z) into 'f/W^-H 



The crucial observation is that, though Ej(z) and ijjmj(z) (j = 1,2) are two- valued 
analytic objects, the analytic continuation ip(z) of the solution (|4.18| ) of the differential 



equation ( |4.13| ) is single- valued J198 
that M 



From this fact it can easily be shown with Eq. (|4.18[) 



aAz) 



cxp 



Ei(z)dz I a 2 (z 



C(z) 



(4.43) 



ioo in the complex plane. Crossing 



13 It is convenient to imagine cuts from zq to zq + ioo and Zq to Zq 
these cuts, one moves from one sheet to the other. 

14 In general this change is accompanied by a factor of geometrical origin. It has no effect in our case 
for the jumping probability [196]. 
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where a%{z) is the analytic continuation of ai(x) in the vicinity of the real axis and 
then along C{z). The Landau - Zener crossing probability is derived from the following 
consideration. Starting on the real axis at <C x rcs , where the evolution of ip is adiabatic 
and where the boundary conditions ai(x-) = 1, a 2 (x_) = hold, and going along a path 
in the complex plane which passes above zq and returning to the real axis at a point 
x + ^> a; res on the second sheet this path can be chosen such that the evolution in the 
complex variable z is adiabatic because we have passed the resonance in safe distance 
and we have thus |ai(a;+)| ~ 1. On the other hand, we can go from x + on the first sheet 
to the point x + on the second sheet via C(x+) and with Eq. (|4.43|) we thus obtain 



P c = \a 2 (x + )\ 2 



cxp 



E\(z)dz 



C(x+) 



(4.44) 



neglecting deviations of |ai(rr+)| from 1. From this formula we can read off that integra- 
tions along the real axis do not contribute to P c . Because of analyticity, the path C(x + ) 
can be deformed such that it goes along the real axis from x + to x res from where it leads 
to zq, circles around the branching point with infinitely small radius and then goes back 
to x res and x + on the second sheet. In this way, we obtain from Eq.( [4.44D the final result 
for the Landau - Zener crossing probability 



InP, 



E 



-0 



Im / dz [(A — Am cos 2$y + {Am 2 sin 2$) 



211/2 



(4.45) 



To get this formula we have used the explicit form of the energies of the adiabatic states 

With the variable transformation A = A(z) and dz = dAj^ it is easy to evaluate 
Eq. ( f4.45| ) for a linear density [ |172| , |187| - |189[ | and is also possible for an exponential density 
|199|| . These calculations can be summarized by ||183| , |189|| 



P c = exp ( --7 res -F 



7T 



(4.46) 



with F = 1 for the linear and F = 1— tan 2 d for the exponential case. For an exponentially 
varying matter density the dependence of P c on $ is particularly simple because 7 res F = 
45 sin 2 $ with 5 = 7 res cos2'i9/ sin 2 2t? which is independent of the vacuum mixing angle. 
For a more general discussion of the crossing probability ( |4.45| ) see Ref. |200| 



Exact solutions of the differential equation (|4 . X 3|) for neutrino oscillations in matter 
exist not only for the linear case [ 187 . |188| , 195 , 196| in terms of Weber functions but also 



for the exponentially varying matter density in terms of Whittaker functions [193, 201 



This case is of particular importance for solar neutrinos because it approximates the 
real density variation in the sun. Further exact solutions are known for A varying with 
tanhx ||202|| and 1/x ||189|| . In Refs. [|183j , |189|| a list of the factors F for all these cases is 
given, calculated with the Landau - Zener formula ( ^4.45| ). 

The survival probability (|4.41| ) with P c in the Landau - Zener approximation does 
not reproduce well the exact survival probability in the extremely non-adiabatic region 
193| , |195|| (see also Refs. ||172| , |183j , |189|| ). In these references the example of an extremely 
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dense medium (A — > oo and thus = ft/ 2) with a sharp boundary to the vacuum is 
given. In this case the neutrino does not oscillate in the medium and therefore P Ve ^ Ve = 
1 — 5 sin 2 2$ stems purely from vacuum oscillations. This has to be compared with P c — 1, 
because of the jump in density one gets 7 res = 0, inserted into Eq. (|4.41| ). Obviously, the 
resulting expression P Ve ^ Ve = cos 2 $ does not agree with the previous one. 

A remedy of this deficiency was found in the framework of the exact solution for an 
exponentially varying matter density leading to the following modification ||193| of the 
Landau - Zener crossing probability P c (|4.45|) : 

exp ( -^7 res F ) - exp ( -^7 res - 



„ /res--- I l „ /res . 2 a 

\ 2 I V 2 sin v, , ,_ 

P <=— 7 u F ^ 

1 " exp {-^-^ 

Numerical calculations for solar neutrinos with an exponential electron density using this 
formula agree very well with numerical integrations of the differential equation (|4.13 ), 



typically the agreement is in the percent range for relevant mixing parameters ||203|| . It 
has been conjectured in Ref. [ 189 ] that the form ( 4.47Q of P c holds for a wide class of 



matter density profiles. In any case, in the limit 7 res — > one gets P c — * cos 2, $ which, 
when inserted into formula ( f4.41|) , correctly describes the survival probability in the above 
example of an extremely non-adiabatic evolution. 

Concluding this section, we consider again an arbitrary number of neutrino flavours 
or types and envisage the interesting case where one of the eigenf unctions ( |4.2U| ) of H™ at 
labelled by the index £q has an adiabatic evolution whereas the other part of ip has an 
arbitrary evolution. We assume exact adiabaticity for ?/w which amounts to Bj£ = 
P>£ k = for j, k £ . This allows to write 

Pu a ^ vp = (1 - ItU^Wna - \U m (x )^\ 2 )Pgi Vp + \U m (xi)pe \ 2 \U m (x U \ 2 , (4.48) 
where we have defined 



Pit,, = [(i-|c/^iKi 2 )(i-i^(xoUI 2 )] 



2m-1 



2 

(xi)/3jB(xi, Xo)jkU^ n (xo)ak 

< 1. (4.49) 

That P^l U/3 is smaller than 1 follows from the Cauchy - Schwarz inequality and from 
the above assumption for B because after dropping the row and the column labelled by 
C.q in B the remaining matrix is again unitary. Eq.( |4.49| ) and similar cases are useful if 
neutrino masses differing by orders of magnitude occur. 

5 Indications of neutrino oscillations 

5.1 Atmospheric neutrino experiments 
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5.1.1 The atmospheric neutrino flux 



In 1912 it was discovered by V.F. Hess [|204|| in a manned balloon flight that the intensity 
of the ionizing radiation in the atmosphere as a function of the altitude did not conform 
with the idea that this ionization was caused by radioactive elements in the surface of the 
earth but rather pointed to an extraterrestrial origin. In the following decades, before 
the advent of accelerator physics, this radiation, which was called first "ultraradiation" 
and later baptized "cosmic rays" by R.A. Millikan, proved to be one of the most fruitful 
means for doing particle physics experiments. At the end of the first half of the 20th 
century such experiments had lead to the discovery of the positron, the pion and the 
muon and also the first particles with strangeness were found with cosmic rays |205 



Eventually, in the beginning of the fifties, proton beams from accelerators replaced the 
cosmic proton flux as an experimental tool. However, after many years where accelerator 
physics was dominating in particle physics, at the end of the 20th century cosmic rays 
play again a major role through atmospheric neutrinos which allow to use the whole globe 
as a neutrino physics laboratory and to probe neutrino mass-squared differences down to 
a few 10~ 4 eV 2 . In this way convincing evidence for the existence of neutrino oscillations 
and thus for non-zero neutrino masses has been obtained ||. 

In a simplified picture, the production of atmospheric neutrinos ||206H21S|| proceeds 
in three steps ||219|| . In the first step the primary cosmic rays |[219| , [220| 1 hit the nuclei 
in the atmosphere, thereby producing charged pions and kaons, either directly or via 
intermediate particles. In the second step, the decay of these particles gives rise to part 
of the atmospheric z/ M and neutrino fluxes: 



7T 



V + V. 



n ■ 
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+ Un and K + — > ji + + z/„ , K — > // + v. 



(5.1; 



In the third step, the v e and v e fluxes and further i/„ and z/„ fluxes are produced by 



e + + v e + Un and \i —> e + u e + z/„ . 



There is also a contribution to the neutrino fluxes from the decays 
K L -> 7r + + t + v t and K + 



7T° + i + + V t 



(5.2) 



(5.3) 



with £ = e,/i and the charge conjugate processes ||215| , [219| , |221|| which does not conform 
with the simple 3-step picture. At low energies the most important process in Eq.( |5.1| ) 
is the pion decay. The contribution of decays Q5.3Q to the neutrino fluxes is small 



222| . |223j| . Let us now describe some useful details of the production of the atmospheric 



neutrino flux and its properties. 



Cosmic rays: Galactic cosmic rays enter the solar system as an isotropic flux of par- 
ticles. There is convincing evidence that the bulk of the radiation with energies less 
than 10 6 GeV comes from our galaxy. Though there is no conclusive proof of the origin 
of cosmic rays yet, plausible mechanisms range from material ejected by supernovae to 
interstellar medium accelerated in supernova shock waves ||224|| . Solar cosmic rays are 
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emitted irregularly by major solar flares on the sun. Apart from electrons, primary cos- 
mic rays consist of protons and bare nuclei (mostly He). Roughly speaking, the chemical 
composition of the galactic cosmic ray particles is given approximately by 90 % H, 9 % He 
and less than 1 % heavier nuclei. However, the chemical composition varies with energy. 
At energies of around 100 MeV per nucleon the particle number ratio H/He is less than 
5, it increases to 10 at 1 GeV and is around 30 at 100 GeV ||225| , [226|1 . The absolute flux 
of cosmic protons is not very large: it is of the order of 1000/m 2 x sec x sr for energies 
of a few GeV and above the atmosphere. As a function of the energy the proton flux 
above a few GeV is well described by the power law E~ 2J until 10 6 GeV [ 225|| . Above an 



energy of 100 GeV per nucleon the cosmic ray fluxes are less precisely known [225, 226 



Solar modulation: The solar cosmic rays (solar wind) are important in so far as they 
weaken and modulate the flux of galactic cosmic ray particles with the solar activity. The 
stronger the solar wind the more difficult it is for the low energy galactic cosmic rays to 
enter the solar sphere of influence. This modulation is noticeable for kinetic energies of 
around 10 GeV per nucleon or less |215| . [224|| . This variation of the primary cosmic ray flux 
with the approximate 11-year cycle of solar activity can be parameterized as a function 
of the neutron monitor at Mt. Washington | |215| |. It induces a corresponding modulation 
of the low energy atmospheric neutrino flux but its influence is small. The effect of the 
solar wind on cosmic ray particles depends on the rigidity of the nuclei which is defined 
by momentum/ charge. It is strongest for cosmic ray particles with small rigidity. For 
neutrino energies E v ~ 2 GeV it is practically negligible, for energies around 1 GeV only 
appreciable at high geomagnetic latitudes because the geomagnetic cut-off admits cosmic 
ray particles with smaller rigidity there. 



Geomagnetic cut-off: 



The geomagnetic field prevents primary cosmic particles with 

In the atmospheric neutrino 



low rigidity from entering the atmosphere [215, 218, 227,228 



flux calculations it is assumed that the momenta of the neutrinos have the same direction 
as the primary cosmic ray particles responsible for their production ||215| , [217|1 . This is a 
good approximation for E u > 200 MeV as shown in Ref. |[229| | . In a first approximation, 
for a given direction, there is a cut-off in rigidity for the primary cosmic ray particles 
below which such a particle cannot reach the top of the atmosphere in a vertical altitude 
of around 20 km. The value of the geomagnetic cut-off can be obtained by a computer 
simulation with the help of the "backtracking technique" . In this technique, to establish 
if a particle with a given charge and momentum coming from a certain direction can 
reach a point on top of the atmosphere, one integrates the equations of motion for a 
particle with opposite charge and reflected momentum starting from its final position. 
If the backtracked particle reaches infinity it is assumed that the trajectory is allowed. 
As infinity one can take a distance of around 10 times the radius of the earth where 
the geomagnetic field has decreased to the level of the interstellar magnetic field with a 
strength ~ 3 x 10~ 8 Tesla ||215 |. More refined calculations replace the rigidity cut-off by a 
probability distribution |218| . [227| . [228| 1. In the rigidity cut-off approximation one can give 



for every experimental location a contour map depicting curves with constant cut-off in 
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a plot with the coordinate axes representing azimuth and zenith angles of the primary 
cosmic ray particle direction. For Kamioka such a plot is found in Ref. |[215|1 . At the 
geomagnetic latitude of 90° the geomagnetic cut-off disappears, therefore experiments 
situated at high geomagnetic latitudes have a higher flux of low energy atmospheric 
neutrinos. The magnetic field of the earth is not exactly of dipole form but also has a 
considerable multipole contribution which affects the local rigidity cut-off maps and thus 



the atmospheric neutrino flux in particular for down-going neutrinos [|215 



Hadronic interactions: After entering the atmosphere the cosmic ray particles collide 
with the nuclei in the air. Though the fractions of other particles than protons in the 
cosmic rays are small, their importance for the atmospheric neutrino flux gets enhanced 
because the predecessors of the neutrinos, pions and kaons, are created through the 
hadronic interactions of the cosmic rays with the air nuclei. These processes depend rather 
on the number of nucleons than on the number of nuclei. The ratio of the charged kaon 
to the charged pion flux depends sensitively on the description of hadronic interactions. 
In this context the quantities 



dx^^±, (5 . 4) 
ax 



where M. = n, K and 7 ~ 1.7 comes from the power law of the primary cosmic pro- 
tons, x is the fraction of the proton momentum carried by the meson and dN/dx is the 
distribution of the charged mesons produced by collisions of protons with nuclei in the at- 
mosphere, are important indicators for the fraction of the muon neutrino flux originating 
from the M. meson flux. The ratio Z p k±/Z p ^± ranges from 0.10 to 0.15 in different calcu- 



lations pi7|| . As a consequence, at neutrino energies below 100 GeV the pions dominate 



as neutrino sources [|216| , |22l| . 



Characteristics of the atmospheric neutrino fluxes: Some of the characteristic 
properties of the atmospheric neutrino fluxes are simple consequences of the produc- 
tion mechanisms (|5.1J ) and (|5.2|) , whereas others follow from a Monte Carlo or analytic 



calculation ||221|| of the neutrino fluxes. As mentioned above, the flux calculations are 



one-dimensional and energy losses in the air, in particular, for the muons, have to be 
taken into account |219|. 



a. Approximately, the fluxes have a power law energy spectrum E v 3 for 1 < E v < 10 3 
GeV whereas the v e fluxes decrease like -E~ 3 ' 5 [ 215fl . 



b. Fixing a neutrino energy E u , the largest contribution to the flux of this energy origi- 
nates from cosmic ray particles with energy E cr ~ 10 x E v [[215 , 216 ]. 



c. Denoting neutrino fluxes by (u) then it follows immediately from Eqs. (|5.1| ) and ( |5.2| 
that 



r p ^2. (5.5) 
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However, it can easily be estimated that, at neutrino energies larger than 1 GeV, 
muons from the reactions ( |5.1|) start to reach the surface of the earth before they 

decay. Therefore the neutrino fluxes from muon decay decrease, in particular, the { v e 



c. 



flux, and the ratio (|5.5|) begins to rise [215, 219|| . 



d. For the same reason muons cease to be the dominant source of [y e 
E u ~ 100 GeV and the reactions (O) take over [pT9|,p2T . 



+ (v e ) above 



(-) 



Pions cease to be the dominant source of the flux for neutrino energies above 100 
GeV p21|| . The reason is given by the different masses and lifetimes of the parent 
mesons M. = it, K in the decays M. — ► \w. For relativistic parent mesons and for a 
given energy E u , the lower bound on the energy of the parent meson in such a decay 
is given by E M > E„/(l — m 2 ^/m 2 M ). Thus, neutrinos with energy E u from the decay 
of the heavier parent meson need in average a lower parent meson energy where the 
corresponding meson flux is higher. This together with the shorter kaon lifetime leads 
to the prevailing of kaons as the source of muon neutrinos at sufficiently high neutrino 
energies gT| , gl| , g21| , g30 |. 



f. The ratio (p e ) / (u e ) is smaller than 1 because among the primary cosmic ray nucleons 
the protons by far dominate which leads to an enhancement of positively charged 
pions over negative ones in the hadronic interactions and, consequently,^ (7r + ) / (7r~) ~ 

> i mis. 



g. On the other hand, for E u <1 GeV, where the muons have time to decay in the atmo- 
sphere, one has {p^/iy^) — 1 because the muons supply the antineutrinos (neutrinos) 
to the neutrinos (antineutrinos) from pions and kaons. However, for neutrino energies 
above a few GeV, when part of the muons reaches the surface of the earth before they 
decay, this ratio becomes smaller than one because the ir + are more numerous than 
the 7r~ (see above). 

h. For E v > 5 GeV the geomagnetic effects become negligible. 

The largest error in the atmospheric neutrino flux calculations arises from an overall un- 
certainty of primary cosmic ray measurements of the order of ±15% ||215| , [217|| . Therefore, 
it is mandatory to use quantities which are ratios such that this uncertainty cancels or, 
in the case of a fit to the data, to include the overall normalization of the primary cosmic 
ray flux in the set of parameters to be determined by the fit. Among the flux calcula- 
tions of the different groups the major source of differences comes from the treatment of 
pion production in collisions of protons with light nuclei ||217|. [231|| . Thus, the calculated 



absolute fluxes have uncertainties of the order of ±20% and even larger uncertainties at 
neutrino energies below 1 GeV, whereas flux ratios differ in general only by a few percent 
(see the comparison between different flux calculations in Refs. [[217| , |23lfl ). 



15 



In general, we denote particle fluxes by putting the particle symbol within parentheses. 
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The angular dependence of the atmospheric neutrino flux: The direction of the 
atmospheric neutrino flux at the location of an experiment is described by the zenith 
angle 9 and the azimuth angle 0. Neutrinos going vertically downward have 9 = 
whereas those coming vertically upward through the earth have 9 = it. The azimuth 
angle indicates the direction of the flux in the horizontal plane. There are two natural 



causes for an angular dependence [218] of the atmospheric neutrino flux: 



a. The development of cosmic ray showers in the atmosphere depends on the density of 
the air. At zenith angles 9 ~ tt/2, where along the line of flight the increase in density 
is less steep, pion and kaon decay is enhanced compared to vertical directions. It is 
easily checked that the situation is symmetric with respect to 9 — > tt — 9. Therefore, 
the generation of pions and kaons in the atmosphere caused by primary cosmic rays 
induces a dependence of the atmospheric neutrino flux on | cos#|. In other words, this 
mechanism does not introduce an up-down asymmetry in the flux. 

b. The geomagnetic field, which acts on the primary cosmic ray flux, is the cause for a 
4> and 9 dependence of the atmospheric neutrino flux. Because of the positive charge 
of the primary cosmic rays the neutrino flux is highest (lowest) for directions coming 
from the west (east). The nature of the geomagnetic effects also generates an up-down 



asymmetry for low energy neutrinos [|215| , |218 



A third cause for a dependence on the zenith angle is possibly given by neutrino 
oscillations which arises because varying 9 from to tt the neutrino path length varies 
from around 10 km [|232| to around 13000 km. Clearly, this dependence is not up-down 



symmetric. If it is disentangled from the up-down asymmetry caused by the geomagnetic 
effects, e.g., by using the high energy component of the neutrino flux, it provides us with 
valuable information on neutrino masses and mixing. 

5.1.2 Experiments with atmospheric neutrinos 



Early efforts for detecting atmospheric neutrinos (see Ref. [ 2"3"3 ] for a summary) concen- 



trated on neutrino-induced upward muons, i.e., muons with zenith angles 90° < 9 < 180°, 
or horizontal muons, i.e., muons with a zenith angle around 90°, using the process 

{ v\ +N^^ + X, (5.6) 



where N is a nucleon in the rock surrounding the detector located underground [ (234|235|| . 



This method was proposed to distinguish muons generated by atmospheric neutrinos 
from muons originating in the meson decays (|5.1|) in the atmosphere. There were two 
experiments, the Kolar Gold Field experiment in India [|236| and an experiment in South 
Africa ||237|| , reporting the first evidence for atmospheric muon neutrinos. Since these 
experiments could not distinguish up and down directions they used horizontal muons 
crossing the detector which lay, due to the great depth of the experimental sites, in the 
large zenith angle intervals 60° < 9 < 120° and 45° < 9 < 135° for the Indian and 
the South African experiment, respectively. These early experiments were accompanied 
already by atmospheric neutrino flux calculations |230| , [238| - |240 . 
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Recent detectors ||241|| are divided into two classes: water Cherenkov detectors where 
the neutrino target is a large volume of water surveyed by a huge array of photomulti- 
plier tubes sitting on the surface of the volume (the Kamiokande [l"9, 242. p" 



Super- 



Kamiokande [[| and 1MB | f)0| , |244| . |245|| collaborations) and iron plate calorimeters where 
neutrino-induced charged particles ionize the gas between the plates and the particle 
paths are reconstructed electronically (the Frejus ||246|| , NUSEX ||247|| and Soudan-2 flSlfl 
collaborations). In contrast to the early detectors the recent detectors are sensitive to 
the direction of tracks and can thus distinguish between up and down through-going 
tracks. However, they cannot measure the charge of the leptons £ and thus cannot dis- 
tinguish between vg and vg. Therefore, at low energies they are approximately sensitive 
the flux combination (yi) + \(vi) (£ = e,//) because at low energies quasi-elastic scatter- 
ing is predominant and the ratio of quasi-elastic antineutrino to neutrino cross sections is 
approximately 1/3 ||217 , 233| , 248]. In atmospheric neutrino physics it is important to dis- 
tinguish e-like and /z-like events which is accomplished by distinguishing between showers 
or diffuse Cherenkov rings (e + ,e~,7) and tracks or sharp Cherenkov rings (fi + ,fi~ and 
charged pions, kaons, protons etc.). The separation between e-like and /i-like events is 
very good, e.g., for Super-Kamiokande its efficiency is estimated to be 98% or better j2[. 
For the deep underground detectors two different event classes are defined 
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Events in which the neutrino interacts with the material inside the detector and where 
all particles from the neutrino interaction deposit their energies inside the detector are 
called contained events. The second class refers to events where muon neutrinos interact 
with the material surrounding the detector via charged current interactions such that the 



high energy muons enter the detector 



muons and stopping muons. Recent measurements of the upward muon flux were per- 
formed by the Baksan [ |250| , Kamiokande ||251| ,|252 

collaborations. 



In this way one distinguishes through-going 
vard muon flux were ] 
Frejus PEfl , MACRO 



1MB 



and Super-Kamiokande |L| 

Super-Kamiokande - and before also Kamiokande - has an inner detector surrounded 
by an outer detector. This allows to further subdivided the contained events into fully 
contained events (FC) with all energy of an event deposited in the inner detector and 
partially contained events (PC) which have exiting tracks detected also in the outer 
detector. In Kamiokande and Super-Kamiokande, FC events are separated into those 
having a visible energy E vis < 1.33 GeV, the sub-GeV events, and those with E v - m > 1.33 
GeV, the multi-GeV events. Among the contained events in Cherenkov detectors the 
single-ring events are well understood, they are predominantly produced by quasi-elastic 
scattering of electron and muon neutrinos. In Kamiokande and Super-Kamiokande for 
the analyses of FC events only single-ring events are used with the additional criteria 
p e > 100 MeV and p^ > 200 MeV for electron and muon momenta, respectively, in 
the case of sub-GeV events. To quote the numbers of Super-Kamiokande @, a Monte 
Carlo simulation has shown that 88% (96%) of the sub-GeV e-like (//-like) events are 
charged current interactions whereas for multi-GeV events the number is 84% (99%). 
The remainder is given by neutral current events. The PC events were estimated to be 

98% ^-induced events for single and multi-ring configurations and, therefore, all PC 
events are used for the analyses. 
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For the analyses of atmospheric neutrino data it is very important to have a good 
understanding of the neutrino interactions in the detector. The deep-inelastic scattering 
(DIS) formulas are only valid for sufficiently high momentum transfer Q 2 from the leptons 
to the hadrons. For neutrinos with a few GeV interacting with the detector material, 
the lowest multiplicity exclusive channels represent an important fraction of the cross 



section. Therefore the following decomposition of the neutrino cross section ||253|| has 
been proposed: 

a u C = a QEL + + gdis , (5.7) 

where QEL indicates the quasi-elastic f£ + n — > t~ +p and i>i + p — » £ + + n cross sections 
and Itt single pion production. In the case of the latter it is assumed that A(1232) 
dominates. To avoid double counting in a\ n + cxdis, the maximal mass of the pion - 
nucleon system in the final state is taken to be W c = 1.4 GeV for single pion production 
whereas the kinematical region for DIS is bounded by W > W c , where W is the mass of 
the hadronic final state. 



5.1.3 The atmospheric neutrino anomaly 



The first observable to be measured in recent atmospheric neutrino experiments was the 
ratio of /i-like to e-like events denoted by (/i/e)data- As discussed earlier, in flux ratios the 
large uncertainty in the overall normalization of the primary cosmic ray flux cancels and 
there is also some cancellation of errors in the theoretical calculation. However, the above 
ratio is only a limited reflection of the corresponding ratio of atmospheric neutrino fluxes 
because of detector efficiencies and event selection criteria. Thus for the expected ratio 
one has to fold the theoretical flux calculations with the cross sections of the neutrino 
interactions in the detector and the detection efficiencies and apply the event selection 
criteria. Quoting the result of the Monte Carlo calculation of Super-Kamiokande as an 
example, this collaboration obtains (/i/e)Mc — 1-50 and 2.83 for the sub-GeV ratio of 
FC events and the ratio considering multi-GeV FC and PC events, respectively [0], using 
the neutrino flux calculations of Ref . ||215|| . Note that these numbers significantly deviate 
from the naive expectation 2. Therefore, the actual physically relevant quantity is given 
by the double ratio 

R = (/Vejdata (5g) 

WejMC 

The first indication that this ratio is smaller than 1 was reported more than ten years 
254J. In the meantime the most impressive measurements of R are represented by 



ago 



the Kamiokande and Super-Kamiokande results:^ 



R 



0.60 



0.57 



+0.07 
-0.06 

+0.08 
-0.07 



± 0.05 sub-GeV [243] 



± 0.07 multi-GeV 



0.63 + 0.03 + 0.05 
0.65 + 0.05 + 0.08 



49] 



sub-GeV 
multi-GeV 



Kamiokande 



Super-Kamiokande 



(5.9) 



3 For recent results on atmospheric neutrinos see Note added. 
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The PC events have been added to the multi-GeV data. The results R = 0.54 ± 0.05 ± 
0.11 of the 1MB Collaboration pj,p55j and R = 0.61 ± 0.15 ± 0.05 of Soudan-2 l256j 
are in agreement with those of the experiments in Kamioka. However, the early iron 
calorimeter experiments found values of R compatible with 1, namely R = 1.00 ± 0.15 ± 

0.08 (Frejus Coll. p46| ) and R = 0.96 +o"il ( NUSEX ColL S)- A P art from the latter 
two experiments, all others hint at a reduction of R compared to the expectation R = 1 . 
Kamiokande and, in particular, Super-Kamiokande have enough statistics to study 

the zenith angle dependence of the measured v e and fluxes. To this end, the cos 9 
interval [—1,1] is divided into five bins of length 0.4. Kamiokande [49|] has observed a 
zenith angle variation of R for the FC multi-GeV + PC events with indications that 
the zenith angle variation rather comes from /i-like events than from the e-like events. 
All this is amply confirmed by Super-Kamiokande with much more statistics and with a 
significant zenith angle variation in the /i-like events for both, sub-GeV and multi-GeV. 

It is important to study the zenith angle variation independent of the double ratio R 
in order to disentangle the zenith angle dependencies of the electron and muon neutrino 
fluxes. To this end, in addition to the oscillation parameters, also the normalization of 
primary cosmic ray flux has to be fitted. Super-Kamiokande has performed a statistical 

analysis of the data under the assumption of v^— >i/ T oscillations. With e-like and /i- 
like events in five cos# bins and seven momentum bins there are altogether 70 data 
points and three quantities to be determined: the mixing angle the neutrino mass- 
squared difference Am 2 and the overall neutrino flux normalization. The best fit gives 
2-d = 1 and Am 2 = 2.2 x 10~ 3 eV 2 with xLn = 65 - 2 for 67 DOF. The regions in the 



sin" 



(-) 



(-) 



sin 2 2t?-Am 2 plane allowed at 68%, 90% and 99% confidence level in the case of z> 
oscillations are shown in Fig. |5.1| 0. At 90% CL the mass-squared difference lies in the 
interval 5 x 10~ 4 eV 2 < Am 2 < 6 x 10~ 3 eV 2 . For the simulation of neutrino oscillations 
the profiles for the neutrino production heights of Ref. [232] were used. An analogous 

oscillations, taking into account matter effects 
87.8 for 67 DOF. 



(-) 



procedure with the hypothesis of z/ M - 
in the earth, gives a poor fit with Xm 



(-) 



It is interesting to note that, as shown in Fig. |5.1|, a fit for Kamiokande with 



(-) (-) 

oscillations gives 5 x 10" 3 eV 2 < Am 2 < 3 x 10" 2 eV r with 90% CL || for the multi-GeV 
data whereas the sub-GeV data show no indication for a zenith angle dependence of the 
number of events. However, Kamiokande has a lower statistics than Super-Kamiokande. 
In Refs. ||257H26lH analyses of all available experiments were performed using the three 

possible oscillation hypotheses v^*^ ^ e ,T,s with similar results as described above and 
showing in addition that also the hypothesis of oscillations into sterile neutrinos gives a 
satisfactory fit. For analyses with three neutrinos and including all neutrino oscillation 
experiments see ||261| - p66|] . In Ref. ||245|| (1MB) no zenith angle variation of R was seen 
though with rather small statistics. 

In the Super-Kamiokande experiment a significant up-down asymmetry for the /i-like 
events was found. The measured value of the asymmetry (U — D)/(U + D) as a function 
of momentum for e-like and /A-like events is shown in Fig. [5.2| [0. Here U is the number 
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Figure 5.1. The 68%, 90% and 99% confidence 



regions in the sin 2'd-Am plane for — > v T os- 
cillations obtained from 33.0 kton year of Super- 
Kamiokande data Q . The 90% confidence region 
obtained by the Kamiokande experiment is also 
shown. 
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Figure 5.2. Up-down asymmetry measured in the 
Super-Kamiokande atmospheric neutrino experi- 
ment as a function of the momentum of e-like 
and /x-like events [||. The hatched region shows 
the theoretical expectation without neutrino os- 
cillations, with statistical and systematic errors 
added in quadrature. The dashed line for /i- 
like events is the theoretical expectation in the 



case of two-generation 



oscillations with 



Am 2 = 2.2 x 10" 3 eV 2 and sin 2 2d = 1.0. 



of upward-going events with zenith angles in the range — 1 < cos 6* < —0.2 and D the 
number of downward-going events with 0.2 < cos 9 < 1. The value of the asymmetry for 
FC and PC multi-GeV /z-like events is [0 

= ( jjq^ j = -0.296 ±0.048 ±0.01. (5.10) 

The geomagnetic effect which is one of the causes of an up-down asymmetry contributes 
less than ±0.01 for multi-GeV events ||. For the PC events with a mean neutrino energy 
of 15 GeV 0] this effect is even less important. Note that for multi-GeV events the 
average angle between the incoming neutrino direction and the charged lepton seen in 
the single-ring events is of the order of 20° or less. Since cos# = ±0.2 corresponds to 
90° ± 11.5°, the asymmetry A of the //-like events is in effect an up-down asymmetry 



C-) 



of the v y. flux. The result ( |5.10|) constitutes the best indication in favour of neutrino 
oscillations found so far with A^ being 6 standard deviations away from 0. This is 
corroborated by the plots of A £}ll (p) as a function of the lepton momentum p in the case 



of the single-ring events shown in Fig. |5.2| 0. For e-like events, A e (p) is a flat function 



consistent with no Ve^v^ oscillations. The asymmetry for e-like events analogous to 



Eq.( |5.10D is compatible with zero: A e = —0.036 ± 0.067 ± 0.02. For the /i-like events, 



A^(p) starts with zero and decreases, thus indicating that at small energies both up and 
down-going neutrinos oscillate with averaged u^— >z/^ probabilities 1/2 whereas at muon 
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energies in the GeV range the probability of v 



(-) (-) 



survival for the down-going neutrinos 



approaches 1. This behaviour is in agreement with the L/ ' E v dependence of the oscillation 
probabilities. It is interesting to note that the experimental values of the asymmetries 
A £ifi can be used to discriminate between different oscillation scenarios ||267| - |270| 1. 

The importance of the results of the atmospheric neutrino oscillation experiments 
requires further scrutiny to test the interpretation in terms of neutrino oscillations. It 
has been proposed for Super-Kamiokande to use ratios of charged current events (CC) 
to neutral current (NC) events [271 , 272] in the spirit of the SNO experiment ||273|| in the 
context of solar neutrinos. The basic idea is that in Super-Kamiokande NC events could 
be seen through 



[ v\ +N -^v\ +vr° + N with N = p,n 



e,/i,r. 



whereas CC reactions are likely to produce a single charged pion via 



{ v\ +N 



£ T + n ± + N with 



(5.11) 



(5.12) 



To produce a test with the reactions ( |5.11| ) and ( |5.12| ) it is necessary to discuss their 
experimental signatures and the contaminations of these signature with other processes. 
Neutral pions are selected by taking events with two diffuse rings from the two decay 
photons with an invariant mass between 90 and 180 MeV and momenta smaller 400 MeV 
in order to separate the two rings. Such events were already observed in Kamiokande 



274 1 . Because of the selection criteria for it s, single pion production is important only 



for neutrino energies of order 1 GeV, therefore it involves the sub-GeV events and no 
tau lepton will produced in the CC reaction. Defining N n o via this procedure then the 
two-ring ratios |[272|| 

A - ' " - A - (5.13) 



and TZ. 



act as measures of the CC to NC event ratios because in the ideal situation one can make 
the identification N n o = N%°, N DS = N™ ± and N ss = N°g ± since electrons produce 



diffuse and muons and charged pions sharp rings. This picture is blurred [ 272 , 275] 
by detector efficiencies (e.g., it is 0.77 for the identification of neutral pions via two 
diffuse rings), misidentifications (like misidentifying a 7i~° whose two photons cannot be 
resolved with an electron) and contaminations (most notably, NC events with 7t + 7t~ in 
the final state contribute to Nss, electronic CC events with a tt° whose photons cannot 
be resolved add to N n o and muonic CC events with the same tt° configuration add to 
Nds)- Taking all this into account ||272j , |275|| , one has nevertheless reason the expect 
that with increasing statistics in Super-Kamiokande the two-ring events can be used to 
obtain information on distinguishing between oscillation of the muon neutrino into tau 
or sterile neutrinos and to discriminate between different regions of the parameter space 
of neutrino mixing ||272| , |275| , |276| |. Furthermore, it was suggested to use the asymmetry 
= (U w o — D w o)/ (U n o + D w o) for up and down-going 7r°'s coming from the NC reaction 
( |5.11| ) as an observable to distinguish muon neutrino oscillations into active neutrinos 



from those into sterile neutrinos ||277|| . One has to assume in this case that up and down 
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7r s originate from up and down neutrinos, respectively, with high probability. This is 
not so obvious, however, since for the identification of the 7r° the events should rather 
have low energy as explained above. 

The further tests discussed here concern the stopping and through-going muon events 
where special efforts have been made to calculate the fluxes ||249|,[253|,|278||. Whereas the 



FC events have neutrino energies of around 1 GeV the stopping muon events correspond 
to a mean neutrino energy of 10 GeV and the through-going muons to 100 GeV ] [23l| , |249| . 
Thus we are now discussing a different range of energy compared to the discussion above 
(with the exception of the PC sample). 

Obviously, also the zenith angle distribution of upward stopping or through-going 
muons can be used to test the neutrino oscillation hypothesis [279]. Among other exper- 



iments (see above) upward through-going muons have been studied by Kamiokande ||252|| 
and upward muons by MACRO ]52|. Kamiokande has 372 such events above an energy 

threshold of 1.6 GeV. Fitting the data to the ( z/^-Vz/j. oscillation hypothesis yields a best 
fit with Am 2 = 3.2 x 10~ 3 eV 2 agreeing rather well with the Super-Kamiokande result. 
The analysis of the MACRO Coll. based on 479 events gives a similar result for Am 2 , 
however, the zenith angle distribution does not fit very well with the oscillation hypothe- 
sis into tau neutrinos. An attempt has been made to explain the zenith angle distribution 

of the MACRO experiment with ► v s oscillations where matter effects in the earth 
play a crucial role ||280| . 

In the earth, the density profile can approximately be represented by constant densities 
in the mantle and the core, respectively (see end of Section 4.3). Such a profile can lead 
to an enhancement of neutrino transitions due to the matter effect^ in the earth if 
atmospheric neutrinos cross the coreQ such that the phase picked up by a neutrino wave 
function traversing the mantle for the first time and the phase acquired by traversing the 
core are each approximately equal to 7r. Such an effect for atmospheric neutrinos was 
recently considered in detail in Refs. |280| , [282| - 1285| . It has also been proposed to exploit 



(-) 



this effect to discriminate between v 



(-) 



fno matter effects" 



and Vp- 



transitions of 



atmospheric neutrinos ||280|| and to explain an excess of e-like events [ [286|| possibly seen 
in the Super-Kamiokande experiment [H. 



It was suggested in Ref. [244 to use the observable 



N. 



stop 



thru 



(5.14) 



where iV s top and iVthm are the numbers of stopping and though-going muons, respec- 
tively, as an indicator for neutrino oscillations because this ratio is reduced for neutrino 



oscillations with respect to the no-oscillation hypothesis [253, 287]. Also the ratios of 
"horizontal" to "vertical" muon events for stopping and through-going muons defined by 



V, 



hor 



stop, thru 



vert 



(5.15) 



stop, thru 



17 Note that for the case of a periodical matter density the effect of enhancement was considered in 



Refs. [ g81 |. 

18 Neutrinos with cos ( 



< -0.837 
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with 

f° dN f cosdc dN 

Nhor = / dcosO- and N vcrt = / dcos9- (5.16) 

J cos8 dcosd 7-i d COS u 



are useful observables as shown in Ref. ||287|| . In this work a study was made for Super- 



Kamiokande, taking cos6 c = —0.5 as the boundary between "horizontal" and "vertical", 
showing that the prospects are good for confirming atmospheric neutrino oscillations 
found with the FC single-ring and PC events by using stopping and through-going muons 
and the observables r (|5.14j) , S stop and Sthru fl5-15| ). However, for a precise determination 



of Am 2 these variables are not suitable. 

5.1.4 Long-baseline experiments and tests of the atmospheric neutrino os- 
cillation parameters 

The explanation of the atmospheric neutrino anomaly in terms of neutrino oscillations can 
be checked with long-baseline neutrino oscillation experiments. The first long-baseline 
reactor experiment CHOOZ ||159|| (see Section 3.2 and Fig. 3.3) has already excluded 
atmospheric v e oscillations with a large mixing angle for Am 2 tm > 10 _3 eV 2 . Two 

other long-baseline reactor experiments are under construction: Palo Verde [|288| , |289| 



and Kam-Land |[290| , |291|| . The Palo Verde experiment has a setup and a sensitivity 



similar to the CHOOZ experiment, whereas the Kam-Land experiment, which is the 
result of the conversion of the old Kamiokande detector to a liquid scintillator detector, 
will detect z/ e 's produced by Japanese reactors 150-200 km away and will be sensitive to 
Am 2 > 10~ 5 eV 2 and a large mixing angle. The Borexino experiment (see the end of the 
Section 5.2 and |[291| , |292fl ) will allow to perform a similar measurement. 

Accelerator long-baseline experiments will study the oscillation channels — + v e ^ T . 



The K2K [ 271 1 experiment, with a baseline of about 235 km from KEK to Super- 
Kamiokande and a neutrino beam with 1.4 GeV average energy, will be sensitive to 
disappearance and v ^ — > u e transitions with Am 2 >2x 10~ 3 eV 2 . A near 1 kton water- 
Cherenkov detector will be placed at a distance of about 1 km from the beam dump and 
will allow to measure the initial flux and energy spectrum of z/^'s. This experiment is 
under construction and is planned to begin taking data in the year 1999. 

Also the MINOS | |293| | experiment is under construction. This experiment will have a 
near detector at Fermilab and a baseline of about 730 km from Fermilab to the Soudan 
mine, where the far detector will be placed. The neutrino beam will be produced by 
protons from the new Main Injector at Fermilab and will have an average energy of 
about 10 GeV. The far detector is an 8 kton sampling calorimeter made of magnetized 
iron and scintillators. This experiment will be sensitive to disappearance and z/ M — > z/ e , 
Vp, v ti v \i v s transitions, with the possibility to distinguish the different channels, for 
Am 2 > 10~ 3 eV 2 (the possibility to extend the sensitivity to Am 2 > 5 x 10~ 5 eV 2 lowering 
the neutrino energy is under study). In particular — > v s transitions can be revealed 
through the measurement of a deficit in the NC/CC ratio. The MINOS experiment is 
scheduled to start data-taking around the year 2003. 

The ICARUS experiment [294| in Gran Sasso, constituted of a 0.6 kton liquid argon 



detector is scheduled to start in the year 2000. In the future three new modules with a 



49 



total mass of 2.4 kton will be installed. This detector will be sensitive to atmospheric and 
solar neutrinos and will allow to reveal long-baseline — > v T oscillations using a neutrino 
beam produced at CERN about 730 km away. Since the average energy of the neutrino 
beam is rather high, about 25 GeV, in order to allow the detection of v T through the 
CC production of a r, this experiment will be sensitive to Am 2 > 10 _3 eV 2 . Four other 
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detectors for future LBL CERN-Gran Sasso experiments, OPERA R295|| , NOE 
AQUA-RICH [[29711 and NICE R298H , have been proposed and are under consideration 



(see Ref. [ 299 ]), together with the possibility of a new atmospheric neutrino detector 
consisting of a large-mass and high-density tracking calorimeter [[300] . 

In Ref. |P01|1 a comparison is made between the possibilities using atmospheric neu- 
trinos and LBL neutrino experiments for the determination of the oscillation parameters. 
The feasibility to distinguish between atmospheric — > v T and z/ M — > v s oscillations 
using a combination of the results of future atmospheric, LBL and SBL experiments is 
discussed in Ref. | |302| . 



5.2 Solar neutrino experiments 

The earliest indication in favour of neutrino oscillations was obtained about 30 years 
ago in the radiochemical solar neutrino experiment by R. Davies et al. ||36|| . The flux 
of solar electron neutrinos measured in this experiment was significantly less than the 
predicted one. This phenomenon was called solar neutrino problem. The existence of 
this problem was confirmed in all five solar neutrino experiments (Homestake p6|-|38|, 
Kamiokande f5|-gT|, GALLEX |E],|4]|, SAGE g§g§ and Super-Kamiokande g§-g8] 



GALLEX |E],g3|, SAGE g§ 
which measure a flux of electron neutrinos significantly smaller than the one predicted by 
the Standard Solar Model (SSM) pj |ITl|iT5]|3^3H3T0 



The solar neutrino problem (see, 
for example, [|60| , |311| - pl6|| ) arose in the Homestake experiment by the low counting rate 
showing that the flux of the high-energy 8 B neutrinos {Es-q < 15 MeV) and of the medium- 
energy 7 Be neutrinos (Etq c = 0.862 MeV) is suppressed by a factor of about 1/3 with 
respect to the SSM prediction. In 1988 the solar neutrino problem was confirmed by the 
results of the real-time water-Cherenkov Kamiokande experiment [39| which measured a 
flux of 8 B neutrinos of about half of the SSM flux. The measurements of the Kamiokande 
experiment proved that the observed neutrinos arrive at the detector from the direction of 



the sun. In 1992 the radiochemical GALLEX [42] and SAGE experiments [45] succeeded 
in measuring the neutrino flux with a low energy threshold E t h = 233 keV, which allowed 
to detect low-energy pp neutrinos produced by the fundamental reaction of the pp cycle. 
Also these experiments measured a neutrino flux of about half of the one predicted by 
the SSM. Finally, the Super-Kamiokande experiment has recently confirmed ]^6| -|48|] with 
high statistics the suppression of the 8 B neutrino flux with respect to the SSM one by a 
factor of about 1/2. 

The energy of the sun is produced in the reactions of the thermonuclear pp and CNO 
cycles shown in Figs. |573| and |5]4] (see, e.g., Ref. |K]]). The overall result of both cycles is 
the transition 

4p + 2e~ ^ 4 He + 2z/ e + g, (5.17) 
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Figure 5.3. The pp cycle. 
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Figure 5.4. The CNO cycle. 
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Table 5.1. Sources of solar neutrinos [317-31E] 



where Q = 4m p + 2m e — m4 Hc = 26.73 MeV is the energy release.^ Hence, the production 
of energy in the sun is accompanied by the emission of electron neutrinos. The main part 
of the solar energy is radiated through photons and a small part (about 2%) is emitted 
through neutrinos. The sources of solar neutrinos are listed in Table |5.1| . The pp, pep, 
7 Be, 8 B and hep reactions belong to the pp cycle (see Fig. [5.3|) , whereas the 13 N, 15 



and 17 F reactions belong to the CNO cycle (see Fig. |5.4| ), which produces only about 
2% of the solar energy. The average and maximum neutrino energies listed in Table |5~I 
are taken from Refs. | |317||318[] . The neutrino fluxes and predictions for the neutrino 
capture rates in the chlorine Homestake experiment and in the gallium GALLEX and 
SAGE experiments given by the Bahcall-Pinsonneault 1998 (BP98) ||304|| SSM are listed 



in Table 5.2 



As it is seen from the Tables 154] and the major part of solar neutrinos are low 
energy neutrinos coming from the pp reaction. Monoenergetic neutrinos with intermediate 
energy are produced in the capture of electrons by 7 Be and in the pep reaction. High 
energy neutrinos are produced in the decay of 8 B (the flux of hep neutrinos is so small that 
its contribution to the event rates of solar neutrino experiments is negligible). The flux 
of 8 B neutrinos is much smaller than the fluxes of pp, 7 Be and pep neutrinos. However, 
as we will see later, these neutrinos give the major contribution to the event rates of 
experiments with a high energy detection threshold. The CNO 13 N, 15 0, 17 F reactions 
are sources of intermediate energy neutrinos with a spectrum that extends up to about 
1.7 MeV. Their contribution to the event rates of solar neutrino experiment is small but 
not negligible. 



19 Here m p = 938.272 MeV is the proton mass, m e — 0.511 MeV is the electron mass and m4 Hc = 
2m p + 2m n — E^L is the mass of the 4 He-nucleus, where m n — 939.566 MeV is the neutron mass and 
Ei Ke = 28.296 MeV is the binding energy of the 4 He-nucleus. 
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Table 5.2 . Standard Solar Model [304] neutrino flux es, average neutrino cross sections |pC, 317 , 318 and 
SSM predictions for the neutrino capture rates 1 304 ] in the chlorine (CI) Homestake experiment and in 
the gallium (Ga) GALLEX and SAGE experiments. 



The neutrino flux coming from each source as a function of the neutrino energy E can 
be written as 

<j) r (E) = $ r X r (E) (r = pp,pep, 7 Be, 8 B,/iep, 13 N, 15 0, 17 F) , (5.18) 

where $ r is the total flux and X r (E) is the energy spectrum (J dE X r (E) = 1). The 
energy spectrum X r (E) for each source r is known with negligible uncertainties [|317| , |320|| 
because it is determined by the weak interactions and it is practically independent from 
solar physics. On the other hand, the total flux $ r of each source r must be calculated 
with a solar model and the resulting uncertainties represent one of the main problem 
for the interpretation of the experimental results. However, there are some relations 
that allow to extract model-independent information on the neutrino fluxes from the 
experimental data. The main one is the luminosity constraint, which is based on the 
assumption that the sun is in a stable state (the energy is produced in the central region 
of the sun and for its electromagnetic part it takes more than 10 4 years to reach the 
surface, whereas neutrinos escape the sun in about two seconds). Let us consider a 
solar neutrino with energy E. The luminous energy released together with this neutrino 
is Q/2 — E. Multiplying this quantity with the total flux of neutrinos ^2 r 4> r (E) and 
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integrating over the neutrino energy one obtains the luminosity constraint^ 



(E)r $ 



Co 







AitR 2 



(5.19) 



where (E) r = J E X r (E) dE is the average energy of the neutrinos from the source r (see 
Table jO|), £ = 2.40 x 10 39 MeV s" 1 [|5j is the luminosity of the sun, R = 1.496 x 10 13 cm 



[ 35|] is the sun 
form 



earth distance. The luminosity constraint can be rewritten in the compact 

Q r $ r = K e , (5.20) 



where Q r = Q/2 - (E) r and K Q = C e /AirR 2 = 8.54 x 10 11 MeV cm" 2 s" 1 is the solar 
constant. Let us emphasize that the luminosity relation is valid under the assumption 
that solar z/ e 's on their way to the earth do not transform into other states. Neglecting 
the small fraction of energy carried away by neutrinos, the luminosity constraint gives 
the approximate value of the total solar neutrino flux $ = ^ r $ r : 



2K C . 







Q 



6.4 x lO^cnr^V 1 



(5.21) 



Since the 3 He nuclei necessary for the formation of 7 Be and 8 B are created by the 
pp or pep reactions, there is another model-independent constraint for the solar neutrino 
fluxes of the pp cycle (see ||321|| ): 



'Be 



+ 0«B < 



+ 



'pep 



(5.22) 



Let us now consider the experimental data. The results of five solar neutrino experi- 



ments are available at present and are listed in Table |5~3 
Homestake 



36 3 



GALLEX ||,(4| and SAGE |0[|45| are radiochemical exper- 
iments. In the pioneering chlorine Homestake experiment of R. Davis et. ai, which 
started in 1967, the detector is a tank with a volume of 6 x 10 5 liters filled with C2CI4. 
Radioactive atoms of 37 Ar are produced by solar electron neutrinos through the reac- 
tion 



324,325 



+ 37 C1 



+ 37 Ar, 



(5.23) 

which has an energy threshold E t ^ = 0.81 MeV. The radioactive 37 Ar atoms that are 
created during the time of exposition of each run (about two months) are extracted from 
the detector by purging it with 4 He and counted in small proportional counters which 
detect the Auger electron produced in the electron-capture of the 37 Ar nuclei. About 0.5 
atoms of 37 Ar are produced every day by solar neutrinos and about 16 atoms are extracted 
in each run (this number is smaller than 30 because of the 37 Ar lifetime of about 35 days 
and because of the extraction efficiency of about 90%). Since the energy threshold is 



20 The corrections to this relation due to the generation of gravitational energy [174] and due to the 
fact that the abundance of 3 He nuclei is out of equilibrium in the outer regions of the solar core [321 



are estimated to be less than 1% and thus negligible at the present level of accuracy. 
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Experiment 



Result 



Theory 



Result 
Theory 



Homestake 38| 




2.56 ±0.16 ±0.16 
(2.56 ±0.23) 


7 7+I.2 
' • ' -1.0 


n oq +0.06 
uo ° -0.05 


GALLEX p22] 




77.5 ± 6.2 ^ 

(78 ±8) 


129 +| 


0.60 ±0.07 


SAGE [323] 




66 6 + 6 - 8 +3 - 8 

UU.O _ 71 _ 4 o 

(67 ± 8) 


129 ±| 


0.52 ±0.07 


Kamiokande 41 


1 


2.80 ±0.19 ±0.33 
(2.80 ±0.38) 


5 15+ 1 - 


0.54 ±0.07 


Super-Kamiokande [48 


2.44 ± 0.05 
(2-44 ^ 9 )' 


5-15 Iq!? 


47 +0 ' 07 

u "^' -0.09 



Table 5.3. The results of solar neutrino experiments confronted with the corresponding theoretical 
predictions [304]. The results of the Homestake, GALLEX and SAGE experiments are expressed in terms 
of event rates in SNU units (1SNU = 10 _36 eventsatom" 1 s _1 ), whereas the results of the Kamiokande and 
Super-Kamiokande experiments are expressed in terms of the 8 B neutrino flux in units of 10 cm s . 
The first experimental error is statistical and the second is systematic. The experimental values in 
parenthesis have the statistical and systematic errors added in quadrature. 



above the end-point of the pp neutrino spectrum and the cross section of the detecting 
process ( |5.23|) grows with the neutrino energy, the main contribution to the counting rate 
in the Homestake experiment comes from 8 B and 7 Be neutrinos. According to the BP98 
SSM ||304|| the event rate in the Homestake experiment should be0 7.7 ij'o SNU, of which 
5.9 SNU come from 8 B neutrinos, 1.15 SNU are produced by 7 Be neutrinos and 0.7 SNU 
are due to pep and CNO neutrinos (see Table [5.2|) . The Homestake event rate presented 
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in Table |T3] is the event rate averaged over 108 runs 

In the radiochemical gallium experiments GALLEX and SAGE electron neutrinos 
from the sun are detected through the observation of radioactive 71 Ge that is produced 
in the process 

(5.24) 



± 71 Ga 



± 71 Ge. 



The GALLEX detector is a tank containing 30.3 tons of 71 Ga (100 tons of a water solution 
of gallium chloride), whereas the SAGE experiment uses about 57 tons of Ga in metallic 
form. The event rate measured in the GALLEX experiment is 0.699 ± 0.069 events per 
day. (Since the 71 Ge lifetime is about 16.5 days, only around 7 atoms of 71 Ge are extracted 
in each 3- weeks run.) 



Since the threshold of the process ( [5.24|) is E. 



th 



= 0.233 MeV, neutrinos from all 
sources are detected in gallium experiments. According to the SSM, the contributions 
to the total predicted event rate from pp, 7 Be and 8 B neutrinos are ||304|| 54%, 27% and 
10%, respectively. 

As can be seen from Table |5.3j, both gallium detectors measure an event rate that is 



One Solar Neutrino Unit (SNU) is defined as 10 ' ib events atom 
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about one half of the SSM prediction. The weighted average of the GALLEX and SAGE 
results yields the event rate 

= 72.5 ± 6 SNU . (5.25) 

Taking into account the experimental and theoretical errors, this rate differs from the 
SSM prediction by about seven standard deviations! 

Both the GALLEX and SAGE detectors have been calibrated using an intense 51 Cr 
neutrino source. The ratio of observed and expected events is 0.93 ± 0.08 for GALLEX 



326|| and 0.95 ±0.12 for SAGE [fl5]| . These results demonstrate the absence of unexpected 
systematic errors at the 10% level in both experiments. In addition, the GALLEX Col- 
laboration calibrated the detector by introducing a known number of radioactive 71 As 
atoms in the target solution [ |326] , |327|| . The atoms of 71 Ge resulting from 71 As decay have 
been extracted in the usual way, and the As tests prove, at the 1% level, the reliability 
of the technique [ |328| | (the number of 71 Ge atoms produced in these tests is of the order 
of 10 5 , whereas the number produced in the 51 Cr tests is of the order of ten per day). As 



emphasized by the GALLEX Collaboration [ 326 , 327f , these results rule out the presence 



of unexpected radiochemical effects that could explain the deficit of solar z/ e 's measured 
by the gallium experiments. 

It is necessary to emphasize that the gallium experiments are not only very important 
for the assessment of the solar neutrino problem but also for the theory of thermonuclear 
energy production in the sun: They have provided the first observation of low-energy 
solar neutrinos produced in the pp reaction that is the basic reaction of the pp cycle and, 
therefore, the first direct experimental confirmation of the theory of thermonuclear origin 
of solar energy production. 

In the Kamiokande [|3Ti|-|4T|j and Super-Kamiokande |46|-f48"| experiments water- 



Cherenkov detectors are used for the detection of solar neutrinos through the observation 
of the Cherenkov light emitted by the recoil electrons in the elastic-scattering process 

v + e~ -» v + e~ . (5.26) 

Since the direction of the recoil electron is peaked in the direction of the incoming neu- 
trino, water-Cherenkov detectors measure the direction of the neutrino flux and the results 
of Kamiokande and Super-Kamiokande have proven that there is a flux of high-energy 
neutrinos coming from the sun with an intensity about half of that predicted by the SSM. 

The energy threshold of water-Cherenkov detectors is given by the threshold for the 
detection of the recoil electron in the reaction Q5.26|) . It is higher than in other solar 



neutrino experiments because of the large background at low energies: E^ m ~ 7 MeV in 
the Kamiokande experiment and ~ 6.5 MeV in the Super-Kamiokande experiment. 
Therefore, only 8 B neutrinos can be detected in water-Cherenkov experiments. If nothing 
happens to electron neutrinos during their trip from the core of the sun to the earth (i.e., 
no neutrino oscillations or other transitions), the results of the Kamiokande and Super- 
Kamiokande experiments provide a measurement of the total flux $sb of 8 B neutrinos. 
Therefore, the results of the Kamiokande and Super-Kamiokande experiments are usually 
presented in terms of the measured flux of 8 B neutrinos. In the Super-Kamiokande 
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experiment it was found 



= (2.44lg;$) x !0 6 cm 



,-2 s -l 



(5.27) 



which is about one half of the flux predicted by the SSM (see Table [5~3|).P| 

From the experimental results and the most updated theoretical predictions ||304|| 
listed in Table |5.3| one can see that the observed event rates in all solar neutrino ex- 
periments are significantly smaller than the predicted rates (see also ||329|| where the 
comparison with other theoretical predictions is discussed). This discrepancy constitutes 
the solar neutrino problem. The SSM is robust and it has been recently tested in a con- 
vincing way by comparing its predicted value for the sound speed in the interior of the 
sun with precise helioseismological measurements (see ||304| , |315| , |330| - |333| ) . However, it 
is clear that a model-independent proof of the existence of the solar neutrino problem 
would be more convincing. Such a model-independent approach has been discussed in 
several papers ||308| , |312j . |321| . |329| . P3 -4H34311 using the model- independent luminosity con- 
straint ( |5.20| ) and the fact that the energy spectra X r (E) Eq. ( |5.18| ) of the various neutrino 
sources are practically independent from solar physics ||317| , [320|| . 

From the luminosity constraint it is possible to obtain a model-independent lower 
bound on the gallium event rate. Indeed, since pp neutrinos have the smallest average 
energy and the u e - 71 Ga. cross section o-Q a (E) increases with the neutrino energy E, we 
have 

(E) <f>r( E ) dE = D(^Ga>r $r > (<T Ga ) pp $ , (5.28) 
r r 

where (o"G a )r = / o'Ga(-E') X r (E) dE is the average cross section of neutrinos from the 
source r (see Table |5.2|) . Since the luminosity constraint ( f5.20| ) implies that 



Q PP $>K e 



j 



(5.29) 



we obtain 



S Ga > 



\0'Gs) P pK (l 



Qpp 

= 11.7 x 10" 46 



76 ± 2 SNU , 



(5.30) 



cm 2 given in Table |5.2| and Q 



'pp 



where we used the value of (cxGa)pp 
13.1 MeV. 

The lower bound ( 5.30|) is just compatible with the combined result ( [5.25 ) of the 
gallium experiments. This means that the results of the GALLEX and SAGE experiments 
can be explained if practically only pp neutrinos are emitted by the sun. This possibility 
is incompatible with any solar model constrained by the helioseismological data [[3Qj| . |315| • 

More stringent model-independent conclusions on the existence of the solar neutrino 
problem can be obtained by comparing the results of different solar neutrino experiments. 
If the survival probability of solar z/ e 's is equal to one, the result of the Super-Kamiokande 
experiment gives the value (|5.27|) for the flux of 8 B neutrinos, whose contributions to the 



2 For recent Super-Kamiokande results on solar neutrinos see Note added. 
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event rates of the chlorine and gallium experiments are 



5 c f ' SK = (a a >8 B = 2.78 ± 0.27 SNU , (5.31) 

4a SK = (^Ga)BB *Sg = 5.9 ^ SNU . (5.32) 

Subtracting the contribution ( |5.31| ) of 8 B i/ e 's from the event rate measured in the 
Homestake experiment (see Table |5.3|) , one obtains 

Sc? ~ S S' SK = -°- 22 ± °- 35 • ( 5 - 33 ) 

Since the quantity Sci — SqP represents the contribution to the chlorine event rate due 
to pep, 7 Be, hep and CNO neutrinos, it must be positive (or zero) but the result in 
Eq.( |5.33D is only marginally compatible with a positive value (the probability is less than 



26%). Furthermore, the result Eq. ( |5.33|) shows that the fluxes of intermediate energy 
neutrinos (pep, 7 Be and CNO neutrinos) are strongly suppressed with respect to the 
SSM prediction (see Table |5.2j). This is incompatible with any solar model constrained 



by the helioseismological data ||304| , pl5|1 and with the fact that, since both the reactions 
that produce 7 Be and 8 B neutrinos originate from 7 Be nuclei (see Fig. |5.3| ), it is very 
difficult to have a suppression of the 7 Be neutrino flux (with respect to the SSM value) 
that is stronger than the suppression of the 8 B neutrino flux.Q 

Let us consider now the gallium experiments. Subtracting the contribution of 8 B 
neutrinos from the luminosity constraint (]5.20p, we have 



K Q - Qs B $8 B = & $r < Q PP ® r ■ ( 5 - 34 ) 

r^ 8 B r^ 8 B 



Hence, following the same reasoning as in Eq. (|5.28 ) we obtain 



= S£' SK + (-Ga) PP K ®-^^ = 82 ±1 SNU , (5.35) 

and the contribution of pep, 7 Be, hep and CNO neutrinos to the gallium event rate cannot 
be bigger than 

^Ga P -^ n = -9.5±6. (5.36) 

This result is compatible with a positive value with a probability smaller than 6%. Again, 
the result ( |5.36| ) implies that the fluxes of intermediate energy neutrinos (pep, 7 Be and 
CNO neutrinos) are strongly suppressed with respect to the SSM prediction, in contra- 
diction with any solar model constrained by the helioseismological data and with the 
moderate suppression of the 8 B neutrino flux. 

Further model-independent methods for proving the existence of the solar neutrino 
problem on the basis of the data of solar neutrino experiments are nicely discussed in 

23 Non-standard solar models with low central temperature yield a suppression of the 8 B neutrino flux 



that is stronger than that of the 7 Be neutrino flux (see [312] and references therein). 
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Figure 5.6. Vacuum v e — > or v e — ► v T oscilla- 
tions: result of the fit of the experimental event 



rates listed in Table 5.3 . Th e shadowed regions 
are allowed at 99% CL |32£|. The dot indicates 
the best- fit point. 



Ref. [[3i2j| . In the following we will assume that there is a solar neutrino problem that is 
caused by neutrino oscillations. 

The solar neutrino data have been analysed in many papers under the assumption of 
two- neutrino mixing (see ||329| . |342| . |344| 1 and references therein) and in a few papers under 



the assumption of three- neutrino mixing (see ||345| - p48|] and references therein). Here we 



will briefly review the results of the most updated two-generation analysis | p29j ] of all 
the solar neutrino data, including those obtained in the first 504 days of operation of the 
Super-Kamiokande experiment (see Table |5.3| ) . 

The deficit of solar z/ e 's can be explained in terms of two- generation neutrino mixing 
either through vacuum oscillations or through MSW resonant transitions in matter ||163 



164 . Furthermore, the two mixed neutrinos can be the electron neutrino and another 
active iv^ or u T ) neutrino or the electron neutrino and a sterile neutrino. There are two 
differences between transitions of solar z/ e 's into active and sterile neutrinos: 

1. The probability of MSW transitions is different in the two cases, because sterile neu- 
trinos do not have the weak-interaction potential due to elastic forward scattering that 
is present for active neutrinos. 

2. The active neutrinos (z^'s or ^r's) that are produced by the oscillation mechanism 
(either in vacuum or in matter) contribute to the event rates measured in the Kamio- 
kande and Super-Kamiokande water- Cherenkov experiments (notice, however, that the 
cross section of v^-e~ and v T -e~ scattering is about six times smaller than the v e -e~ 
cross section). 

The formulas for the survival probability of solar z/ e 's in the case of vacuum oscillations 
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and MSW transitions are given in Eq.( |3.21| ) (with a = e) and in Eq.( |4.41| ) together with 
Eq.( f4.47|) | 193|| ) respectively. These formulas depend on the two mixing parameters Am 2 
and sin 2 2$. The allowed regions in the sin 2 2$-Am 2 plane, obtained in Ref. [329] from 



the fit of the measured event rates listed in Table |5.3| by using the BP98 SSM |[304j| and 



the analytic formulas in Ref. [ |193|| for the MSW case, are shown in Figs. |5.5h fT6 

Figure |5.5| ||329|| shows the three allowed regions in the case of MSW v e 
v e — > v T transitions. They are the small mixing angle (SMA-active) region at 



Am 2 ~5x l(T 6 eV' 



sin 2 2d ~ 6 x 10" 



and the large mixing angle (LMA) region at 

Am 2 ~ 2 x 10~ 5 eV 2 , sin 2 2^9 ~ 0.76 



(SMA-active) , 



(LMA) 



or 



(5.37) 



(5.38) 



The best fit, with a confidence level of 19%, is obtained in the SMA region, whereas the 
LMA region has a confidence level of 4%. In Ref. ||329|1 there is also the so-called low 
mass (LOW) region at Am 2 ~ 8 x 10 _8 eV 2 and sin 2 2$ ~ 0.96, however, the LOW region 
is only marginally acceptable, with a confidence level of 0.7%. 

The allowed regions in the sin 2 2i9— Am 2 plane in the case of v e — * or u e — > u T 



vacuum oscillations are shown in Fig. |5.5| ||329|| . These regions extend over large ranges 
of sin 2 2$ and Am 2 around the best fit values 



Am 2 ~ 8 x 10 _11 eV 



sin 2 2$ ~ 0.75 



(Vac. Osc.) . 



(5.39) 



The confidence level of the fit is 3.8%. There is no allowed region in the case of v e — > v s 
vacuum oscillations (the fit has a confidence level of 0.05%) ||329| , |349|| . 

As shown in Fig. |5.7| ||329|| , only the small mixing angle (SMA-sterile) region at 



Am 2 ~ 4 x 10" 6 eV 2 



sin 2 2t?~7x 10" 



(SMA-sterile) 



(5.40) 



is allowed in the case of MSW v e — > v s transitions, with a confidence level of 19% (the 
large mixing angle and low mass solutions have a confidence level of 0.001% and 0.003%, 



respectively) [[329| . |350 



The authors of Ref. ||329|| performed also global fits of the solar neutrino data including 
the energy spectrum and zenith-angle distribution of the recoil electrons measured in 



the Super-Kamiokande experiment |f48| . As noted in Ref. [[329|1 , since these data are still 
preliminary, the results of this analysis are less robust than those obtained fitting only the 
global rates. However, it is interesting to note that in the case of MSW v e — > or v e — > v T 



transitions only the SMA-active region remains allowed, as shown in Fig. |5.8j [ |329|| , with 
a confidence level of 7%. 

During the night, solar neutrinos pass through the earth and the matter effect can 



cause a regeneration of z/ e 's (see Ref. ||316|| and references therein) and, therefore, a zenith 
angle dependence of the solar neutrino flux. The size of this effect depends on the 
values of Am 2 and sin 2 2$ and is sizable only for large values of the mixing angle. The 
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preliminary value of the day-night asymmetry of solar neutrino events measured in the 
Super-Kamiokande experiment is [[48| 



D-N 
D + N 



-0.023 ±0.020 ±0.014 



(5.41) 



compatible with zero. Hence this experimental result tends to exclude solutions of the 
solar neutrino problem with a large mixing angle. 



As was shown in Refs. ||283| , |351]| the step-like profile of the matter density of the earth 
could lead to an enhancement of neutrino transitions not only for atmospheric neutrinos 
(see Section 5.1.3) but also for solar neutrinos. 

In conclusion of this section, we would like to emphasize that the results of solar 
neutrino experiments provide a rather strong indication on favour of neutrino mixing and 
several experiments are under construction |[273j , |322| , |352|| or in project (see ||353| , |354|| and 
references therein). In particular, the measurement of the electron neutrino spectrum 
at SNO may allow to obtain model-independent information on the neutrino oscillation 
probability and on the flux $s B of 8 B neutrinos [p55| - |359| , the measurement of the flux 
of 7 Be neutrinos on the earth in the Borexino experiment and the results of Super- 
Kamiokande, SNO, Borexino and GNO will allow to distinguish the different possible 



solutions of the solar neutrino problem (see Refs. [ |350| , |360| - |363| | and references therein). 

The Sudbury Neutrino Observatory (SNO) ||273| , |364|| is located 6800 feet below ground 
in the Creighton mine, near Sudbury in Ontario (Canada). In the SNO experiment a 
Cherenkov detector with 1 kton of heavy water (D2O) contained in an acrylic vessel of 
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12 m diameter will be used. The Cherenkov light is detected with a geodesic array of 
10 4 photomultiplier tubes surrounding the heavy water vessel. The heavy water detector 
is immersed in normal water in order to reduce the background. Solar neutrinos will be 
observed in real-time through the charged-current (CC), neutral-current (NC) and elastic 
scattering (ES) reactions 

u e + d^e~+p + p (CC), (5.42) 

u t + d->Vi+p + n , £ = e,fi,T, (NC) , (5.43) 

u t + e~ -»• Pi + e~ , £ = e,fi,r, (ES) . (5.44) 

Since the energy threshold for the observation of the recoil electrons in the CC and ES 
processes is about 5 MeV and the neutrino energy threshold for the NC reaction is 2.2 
MeV, only 8 B neutrinos can be observed. The event rates predicted by the SSM are 
around 23 per day for the CC reaction, 7 per day for the NC reaction and 3 per day for 
the ES reaction. From the hit pattern of the photomultipliers, the direction and energy 
of the neutrino in the CC reaction can be reconstructed, allowing a direct measurement 
of the spectrum of the solar electron neutrino flux on the earth. The observation of a 
distortion of this spectrum with respect to the one calculated without neutrino oscillations 
will represent a model-independent proof of the occurrence of neutrino oscillations. The 
NC reaction will allow to detect all active neutrinos with the same cross section (because 
of the e-\L—r universality of weak interactions), whereas the cross sections of and v T 
scattering in the ES reaction is about six times smaller than that of v e . The measurement 
of the neutral current reaction will allow to determine the total flux of active neutrinos 



from the sun ||355| , |359|| , which coincides with the produced flux if there are only active- 
active neutrino oscillations, which can be revealed trough a measurement of a deficit in 
the NC/CC event rate. On the other hand, if the NC/CC event rate will not show any 
deficit and a distortion of the electron neutrino spectrum will be observed, it will mean 
that there are transitions of solar u^s into sterile states ||356|| (of course, complicated 



scenarios with transitions solar z/ e 's into active and sterile states are not excluded and 
must be taken into account). The SNO experiment is expected to release the first results 
in the year 1999. 

The Borexino experiment | |352| , |365| is designed to detect low-energy solar neutrinos in 
real-time through the observation of the elastic scattering process ( |5.44j ) with an energy 



threshold for the recoil electron of 250 keV. It will use 300 tons of liquid scintillator 
in an unsegmented detector with 2000 photomultiplier tubes. The event rate predicted 
by the SSM for a fiducial volume of about 100 tons is about 50 events per day, mostly 
generated by the 0.86 MeV monoenergetic line of Be solar neutrinos. Since this line 
gives a characteristic spectral signature in the ES process, the flux on the earth of 7 Be 
solar neutrinos will be determined and it will be possible to check if it is suppressed 
with respect to the one predicted by the SSM as suggested by the results of current 
experiments. The neutrino oscillation solutions predict a measurable flux: ~ 10 events 
per day for the SMA MSW solution, ~ 30 events per day for the LMA MSW solution 
and ~ 35 events per day for the vacuum oscillation solution (see |360|). The Borexino 



experiment is under construction in the Laboratori Nazionali del Gran Sasso in Italy and 
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is scheduled to start data taking around the year 2001. 



5.3 The LSND experiment 



The Liquid Scintillator Neutrino Detector (LSND) experiment |5^ , |5^ , |366l - p71|1 at the Los 
Alamos Meson Physics Facility (LAMPF) is a SBL neutrino oscillation experiment with 
a baseline of about 30 meters. It utilizes the LAMPF proton beam with a kinetic energy 
of 800 MeV which produces pions by hitting a 30 cm long water target located about 1 m 
upstream from a copper beam stop. Most of the 7r + 's are stopped in the target and decay 
into muons which come to rest and decay in the target as well. The decay at rest of the 
positively charged muons allows to investigate — > v e oscillations p3| , |366| , |367| , |370| , |371| . 
A small fraction of the positively charged pions (~ 3.4%) decays in flight in the 1 m 
long space between the target and the beam stop and is used for the study of — ► v e 
oscillations f54] , [369H"371|| (only ~ 0.001% of the /i + 's decay in flight and produce a small 
contamination of z/ e 's in the beam). The detector is an approximately cylindrical 
tank 8.3 m long by 5.7 m in diameter, filled with a scintillator medium consisting mainly 
of mineral oil (CH 2 ). In the following we review the results obtained by the LSND 
Collaboration from the analyses of the data from fi + decay at rest and 7r + decay in flight. 



fi + decay at rest |5l||366| , |3"67l , [370| , |371|| . In this case one has the decay chain 



7T 



e + + v e + v u . 



(5.45) 



Note that in these decays no v e appears, which makes the study of — > v e transitions 
possible. In /i + decay at rest the upper bound on the energy is m M /2 = 52.8 MeV 
(m M is the muon mass) and the energy spectrum of the flux is very well known, being 
determined by the kinematics of muon decay. The electron antineutrinos are detected by 



V 



n . 



(5.46) 



The v e signature is a coincidence between e + and a delayed 7 from the capture n + p —>■ 
d + 7 (2.2 MeV). 

There is a background of v e from the decay chain symmetrical to Eq.( |5.45[ ) starting 
with 7r~. This v e flux is suppressed relative to the flux from the tt + chain by a factor 
of 7.8 x 10~ 4 for three reasons: ir production in the beam stop is a factor of 8 smaller 
than 7r + production, only 5 % of the tt~ decay in flight (all stopped tt~ are absorbed) and 
of the resulting /i~ only 12 % decay before they are captured by nuclei. In the latter case 
they give rise to but not to v e . The LSND apparatus does not distinguish between e + 
and e~. Therefore, the reaction v e + 12 C — > + n + 11 N gives a signal like a u e but can 
only occur for E e < 20 MeV. 

LSND has seen an excess of e + events with energies between 20 and 60 MeV of 
51.0 1^95 ±8.0 corresponding to a — > v e transition probability of (3.1 ±1.2 ±0.5) x 10~ 3 
In Refs. ||370| , |371]| preliminary results can be found for the data taken from 1993 to 
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Figure |5.10| . The region between the two solid 
lines is allowed at 95% CL by the ir + de cay in 
flight analysis of the LSND data ||,|69|. The 
regions enclosed by dashed lines are the 99% CL 
LSND favoured regions obtained with the 1996 



muon decay at rest analysis |53, 367 1 



1997. The value for the transition probability agrees with the previous one, except for 
a reduced statistical error of 0.9 x 1CT 3 instead of 1.2 x 1CT 3 . The plot of the favoured 



region in the sin 2t?-Am plane obtained by the LSND collaboration is shown in Fig. |T9 
In the positron energy range between 36 and 60 MeV the background is reduced because 
in this energy range there are no electrons from v e + l2 C — > e~ + 12 N induced by electron 



neutrinos from /z + decay at rest [53]. LSND has observed 22 events with such positron 



energies compared to a background of only 4.6 ± 0.6 events [53|. 



7r + decay in flight [|54| , [369H371I . The n + — > ji + + decay in flight generates a 
beam with which — > v e oscillations are studied. The electron neutrinos are detected 
by electrons from the inclusive reaction 

v e + 12 C -> e~ + X , (5.47) 

where the electrons have an energy E e between 60 and 200 MeV. In this case the signal 
is a single electron. The electron neutrinos from fi + — > e + + v e + decay in flight 
are suppressed by the longer muon lifetime and the kinematics of the three-body decay, 
whereas those from 7r + — > e + + v e by the small branching ratio of 1.23 x 10 -4 . The lower 
bound of 60 MeV on E e lies well above = 52.8 MeV, the endpoint of the electron 
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spectrum from muon decay at rest. Such muons can be induced by beam-related 



(-) 

events. Also u e from fi + decay at rest (|5.45|) generate electrons below this energy. LSND 
has seen 40 electron events compared to a background of 21.9 ± 2.1. This corresponds 
to a — > v e transition probability of (2.6 ± 1.0 ± 0.5) x 10~ 3 [54|, consistent with the 



one determined by muon decay at rest. The favoured region in the sin 2??-Am plane 



obtained by the LSND collaboration is shown in Fig. |5.10| as the area between the two 
solid lines ||54|| . The regions in Fig. 5.10 enclosed by dashed lines are the LSND favoured 
regions obtained with the 1996 muon decay at rest analysis ]|53"| , |367|| . One can see that 
the two results are perfectly compatible. 

Recently, the LSND collaboration has reported ||376|| new evidence in favour of — > v e 
oscillations using the beam produced by tt + decay in flight and detecting the z/ e 's 
through the observation of the reaction 



Ve + 12 C 



12 



Ng.8. + e 



(5.48) 



where 12 N g . s . is the ground state of 12 N, which decays with a lifetime of 15 ms into 
12 C + e + + v e . A temporal (< 15 ms) and spatial (< 30 cm) correlation between the 
observations of the e~ and e + allows an almost background-free measurement. The LSND 
collaboration reported [376] the observation of 5 events, with an expected background of 
0.5 events. These numbers are in agreement with the results of the other two analyses of 
the LSND data summarized above. 

The region of neutrino mixing parameters indicated by the results of the LSND ex- 
periments is presently explored also by the KARMEN experiment ||375| , |377j - |382|| . The 
KARMEN experiment is located at the ISIS spallation neutron facility of the Rutherford 
Laboratories (UK). At ISIS a pulsed proton beam is used: two 100 ns wide proton pulses 
separated by 330 ns are produced every 20 ms. The proton pulses are directed on a target 
where they produce positive pions (the negative pions are absorbed in the source before 
they decay) which decay at rest according to the chain ( |5.45| ) producing an equal number 



of v e and v^. The KARMEN Collaboration searches for v e at a distance of about 18 
m. The time structure of the neutrino beam is important for the identification of the 
neutrino induced reactions in the KARMEN detector and for the effective suppression 
of the cosmic ray background. The KARMEN experiment started in 1990 and ran until 
1995 as KARMEN 1 
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In 1996 the experiment was upgraded to KARMEN 2, elim- 
inating the main cosmic ray induced background component in the search for — > u, 
oscillations ||379|| . So far, the KARMEN 2 experiment measured no events [|375] . |380| - |382 



also no background events, of which 2.88 ± 0.13 events were expected. The resulting 
exclusion curve obtained with a Bayesian analysis ||375|| is shown in Fig. |5]9|. One can 
see that the null result of the KARMEN 2 experiment excludes part of LSND-allowed 
region in favour of — > v e oscillations, but at present the sensitivity of the KARMEN 
experiment is not enough to check all the LSND-allowed region.P| 



24 The KARMEN 2 result has also been ana lysed using the Unified Approach [|383[ , which is a frequentist 
method with proper coverag e (se e Re f. [384]). The resulting exclusion curve seems to forbid almost all 



of the LSND-allowed region |380, 382], but this result depends on the non-observation of the expected 
background events. The analysis of the KARMEN 2 null result with an alternative frequentist method 
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The KARMEN experiment is expected to reach a sensitivity that will allow to cover 

If the KARMEN experiment will find 



the LSND-allowed region in one or two years [|37S 



an oscillation signal compatible with the one observed in the LSND experiment, the LSND 

evidence in favour of oscillations will be confirmed. On the other hand, if the 

KARMEN experiment will not find any neutrino oscillation signal, then the KARMEN 
exclusion curve will lie on the left of the LSND allowed region in the sin 2 2$-Am 2 plane, 
but not far away from it [37G]. Hence, it is important to test the LSND result with a 
new experiment with much higher sensitivity. Four such experiments have been proposed 
and are under study: BooNE |387| at Fermilab, 1-216 |388| at CERN, ORLaND ||9| at 
Oak Ridge and NESS at the European Spallation Source ||39Q|| . Among these proposals 
BooNE is approved and will start in 2001. 

We would like to finish this section by emphasizing the importance of the results 
obtained in the LSND experiment for neutrino physics. As we will discuss in the next 
section, the LSND evidence in favour of SBL ^v^—^-Vg, oscillations taken together with the 
indications in favour of solar and atmospheric neutrino oscillations require the existence 
of at least one sterile neutrino, which represents a manifestation of new physics beyond 
the Standard Model. Hence it is important that the LSND indication is checked as soon 
as possible by another experiment with high sensitivity. 



6 Analysis of neutrino oscillation data 

As we have discussed in the previous section, at present there are three indications in 
favour of non-zero neutrino masses and mixing. The analysis of the experimental data 
(based on the simplest assumption of two neutrino mixing) indicate that there must be 
at least three different Am 2 , each one relevant for the oscillations of neutrinos in one 
type of experiments: Am 2 ~ 10~ 5 eV 2 (MSW) or Am 2 ~ 10 _10 eV 2 (Vac. Osc.) is 
relevant for oscillations of solar neutrinos, Am 2 ~ 10~ 3 eV 2 is relevant for oscillations 
of atmospheric neutrinos, and Am 2 ~ 1 eV 2 is relevant for SBL neutrino oscillations 
according to the LSND experiment. Three different Am 2 require the mixing of at least 
four massive neutrinos. Here we consider two possible scenarios: 

1. Only the indications in favour of neutrino mixing obtained in atmospheric and solar 
neutrino experiments are confirmed by future experiments. In this case the mixing of 
three massive neutrinos can accommodate all data. 

2. All the existing indications in favour of neutrino mixing, including that obtained in 
the LSND experiment, are confirmed by future experiments. In this case, the mini- 
mal mixing scheme that allows to accommodate all data is one with mixing of four 
neutrinos.^]. 



shown in Fig. |J fg| , g8e | 



with proper coverage |385| leads to an exclusion curve compatible with the Bayesian exclusion curve 



25 From time to time special schemes with three neutrinos have been proposed in order to accommodate 



all neutrino oscillation data [ p91 - 397 1 However, none of these schemes can fit all data in a satisfactory way 
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6.1 Mixing of three massive neutrinos 

Let us consider the mixing of three massive neutrinos, i.e., n — 3 in Eqs. (|3.3| ) and (|3.9| ). 



'a) = ^U* k \v k ) , \v a ) = ^U ak \v k ) (a = e,fi,r) 



k=i 



k=l 



and assume the neutrino mass hierarchy^ 

mi « m 2 « m 3 



(6.1) 



(6.2) 



that is indicated by the see-saw mechanism (see Section pT5|) and by the models beyond 
the Standard Model (see Section |2.6|) . In this case, the only available possibility for the 
explanation of the solar and atmospheric neutrino anomalies through neutrino oscillations 



is that Afrijj = m 2 — m? is relevant for the transitions of solar neutrinos and Am 2 



.31 



nr. 



is relevant for the transitions of atmospheric neutrinos: 



Am 2 ! ~ 



10~ b eV 2 (MSW) or Am 2 21 ~ 10~ lu eV 2 (Vac. Osc 
10" 3 eV 2 . 



(6.3) 
(6.4) 



Let us consider the oscillations of atmospheric and LBL neutrinos and assume that 



E 



< 1. 



(6.5) 



The distance L travelled by atmospheric neutrinos ranges from ~ 10 — 20 km for 
downward-going neutrinos to ~ 1.3 x 10 4 km for upward-going neutrinos and the energy E 
relevant for the oscillations of atmospheric neutrino experiments ranges from ~ 100 MeV 
to ~ 100 GeV. Therefore, if the solar neutrino problem is due to MSW transitions, the 
condition ( |6.5|) is not satisfied by atmospheric neutrinos with energies below ~ 1 GeV. 
In this case, the following discussion applies only to high-energy atmospheric neutrinos. 
On the other hand, the inequality ( |6.5| ) is valid in any case for reactor LBL experiments 
with a source - detector distance L < 1 km and for accelerator LBL experiments with a 
source - detector distance L < 10 3 km and neutrino energy E > 1 GeV. 

If the inequality fl6.5p is satisfied, from Eq.( |3.12| ), we have for the probability of 
v a v fi transitions 



p(atm,LBL) 



5 a p + U*z U, 



a3 U P3 



exp 



Am?iL 

>, 

2E 



1 



(6.6) 



|p98 , 399]. Schemes with more than four neutrinos have been considered in the papers in Refs. |130| , 
(see the review Ref. [401]). 

26 Actually, the phenomenology of neutrino oscillations discussed in this section depends only on the 
hierarchy Am^ <C ATO32 of mass-squared differences. Hence, the discussion of neutrino oscillations 
presented in this section applies, for example, also to the case of three almost degenerate neutrinos with 
masses in the eV range that give a contribution to the dark matter in the universe. Furthermore, with a 
cyclic permutation of the indices 1, 2, 3, the discussion of neutrino oscillations presented in this section 
can be applied to all cases in which Am| 2 <C Am^, as, for example, the case of an inverted mass 
hierarchy mi <C m,2 — m^. 
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Therefore, a mass hierarchy of three neutrinos implies that the transition probabilities 
in atmospheric and LBL experiments depend only on the largest mass squared difference 
AmJi and on the elements U a3 that connect flavour neutrinos with the heaviest neutrino 



v 3 . From Eq.flOD we obtain p2j|402| -M 



Pt^ BL) = \ A* (l - cos ^§^) (a * 0) , (6.7) 
p ( f ,,lbl) = 1 _l ^ ^ _ cog Am|L^ (g g) 

with the oscillation amplitudes 

A a;/3 = i\U a3 \ 2 \Up 3 \ 2 , (6.9) 
B a - a = A\U a3 \ 2 (l~\U a3 \ 2 ) . (6.10) 

The expressions (|6.7| ) and (|6.8|) have the same dependence on the quantity Am 3 l 1 L/2E as 
the standard expressions (|3.20|) and ( fj.21|) that describe neutrino oscillations in the case 



of two flavours, with Am 2 replaced by Am^. This is due to the fact that only one mass- 
squared difference is relevant for the oscillations of atmospheric and LBL neutrinos in the 
case of the mass hierarchy ( |6.2j ). The expressions ( |6.9p and ( ft. ICQ describe, however, all 
possible oscillations between three flavour neutrinos: v^ t> v e , v^^ v r and v e ^ v T . 

The unitarity of the mixing matrix implies that Y2 a I 1 2 — 1 an d LBL and at- 
mospheric neutrino oscillations in the scheme under consideration are described by only 
three parameters: Am| 1; \U e3 \ 2 and l^^l 2 . 

Let us also stress that the transition probabilities (|6.7| ) are independent on the phases 
of the mixing matrix U. Thus, in the case of a hierarchy of neutrino masses the relation 
( P-18Q is automatically satisfied and effects of CP violation in the lepton sector cannot 



be revealed in neutrino oscillation experiments. 

As it is seen from Eq.(ft.8[), the v e survival probability in reactor LBL experiments 



depends on the parameters Am^ and |t/ e 3| 2 . This fact allows to obtain information on 
the value of \U e3 \ 2 in LBL reactor experiments and, in particular, from the results of the 



first of them, CHOOZ ||159| , in which no indications in favour of neutrino oscillations were 
found. The exclusion curve in the plane of the two-neutrino mixing parameters sin 2 2$ 
and Am 2 that was obtained in the CHOOZ experiment is shown in Fig. |3.3| . 

Comparing Eqs. (|3.21|) and (|6.8| ), one can see that in the three-neutrino scheme under 
consideration in this section 

sin 2 2d = B e . e = 4 \U e3 \ 2 (1 - \U e3 \ 2 ) . (6.11) 

Therefore, we can obtain information on the parameter \U e3 \ 2 from the exclusion plot 
obtained in the CHOOZ experiment. From Eq.( 6.11| ) we have 



|[/ e3 | 2 =^(l±Vl-£e;e) , (6.12) 

and the exclusion curves shown in Fig. |3l^ imply the constraint 

Be;e<Bi e , (6.13) 
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with a value of B®. e that depends on Am^. From Eqs.( |6.12| ) and ( |6.13| ) it follows that 
|[/ e3 | 2 must satisfy one of the two inequalities 



|2 / ^0 



Uj z <a}i or \UA'>l-d 



e3| 



(6.14) 



with 



a° e = ~(l-Jl-Bi, 



(6.15) 



From the CHOOZ exclusion curve in Fig. |3]| (90% CL) we find B e . e < 0.18 for Am^ > 
2 x 10~ 3 eV 2 and therefore 



[7 e3 | 2 < 5 x 10~ 2 or |£/, 3 | 2 > 0.95 



c3 | 



(6.16) 



6.1.1 Solar neutrinos 



Large values of |£7 e3 | 2 as those implied by the second inequality in Eq.( |6.16|) are excluded 
by the data of solar neutrino experiments. In order to prove this statement, let us derive 
the expression for the averaged survival probability of solar electron neutrinos in the case 
of vacuum oscillations. From Eq.( |3.12D with n = 3, taking into account that for solar 
neutrinos Am^L/E 3> 1, the survival probability averaged over the unobservable fast 
oscillations due to Am^ is given by 



p(sun) 



\U el \ 2 + \U e2 \ 2 exp 



Am^L 
2E 



\U, 



(3 1 



This expression can be rewritten in the form 

Pt^=(^\Ues\ 2 ) 2 P^l + \U, 

(1 2) 

where Po e U Ve is the two-generation survival probability 



4 

e3 1 i 



_p(1.2) 



1 - - sin 2 2$ s 



cos 



IE 



(6.17) 



(6.18) 



(6.19) 



(see Eq.flimp) with A 



cos ^ s 



Am 2 x and the mixing angle ^ sun determined by 



sin 2 $ sun 





u 


el 


2 


1 - 




u e3 


2 





u 


c2 


2 


1 - 




u e3 


2 



(6.20) 



The relation ( 6.18 ) is also valid if matter effects are taken into account (see Eq.( 4.48 ) and 
Ref. |406|1 ). In this case, pj^ Ve is the two-neutrino survival probability in matter given 
in Eq.( [4.41| ), with Am 2 replaced by Am^ and by $ sun . 

Considering now the constraints (|6.16|) on |£/ e3 | 2 , if \U e3 \ 2 is large (|£7 e 3| 2 > 0.95), from 
Eq.( |6.18|) one can see that Pul™l e > 0.90. Such a large value of P Ue ~,u e is incompatible 
with the deficit of solar z/ e 's observed in all solar neutrino experiments. Therefore, from 
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the results of the CHOOZ experiment we reach the conclusion that the element \U e s\ 2 is 
small: 

\U e3 \ 2 < 5 x KT 2 for Am 2 31 > 2 x l(r 3 eV 2 . (6.21) 
This limit implies that 

Pl?X - P£'\ and sin 2 2tf sun ~ 4 \U e2 \ 2 (l - \U e2 \ 2 ) . (6.22) 

Hence, the oscillations of solar neutrinos are practically independent of the value of 
\U e s\ 2 and are decoupled from the oscillations of atmospheric and LBL neutrinos, which 
depend only on the elements U a3 of the mixing matrix (see Eqs. (|6.9|) and ( |6.10|) ). This 
means that the two-generation analyses of the solar neutrino data are appropriate also 
in the three-neutrino scheme and they give information on the values of Am^ = Am 2 
|Z7 e i| — cos^ sun and \U e2 \ ~ sin$ sun ||407|| . 



sun' 



According to a recent analysis of the solar neutrino data [|344|| , which include prelimi- 
nary data from Super-Kamiokande |4"6] , |408| | , the ranges of the mixing parameters allowed 
at 90% CL for the small and large mixing angle MSW solutions and for the vacuum 
oscillation solution are, respectively, 

4 x 10" 6 eV 2 < Am s 2 un < 1.2 x 10~ 5 eV 2 , 3 x 10~ 3 < sin 2 2tf sun < 1.1 x 10" 2 ,(6.23) 
8 x 10- 6 eV 2 < Am s 2 un < 3.0 x 10" 5 eV 2 , 0.42 < sin 2 2tf sun < 0.74 , (6.24) 
6 x 10~ n eV 2 < Am 2 un < 1.1 x 10~ 10 eV 2 , 0.70 < sin 2 2tf sun < 1 . (6.25) 

Therefore, taking into account that \U e3 \ 2 is small and assuming that \U e2 \ < \U e i\ (this 
choice is necessary only for the MSW solutions in order to fulfill the resonance condition 
Q4.34j )), we have 



\U e i\ ~ 1 , \U e2 \~ 0.03 - 0.05 (small mixing MSW) , (6.26) 

\U el \ ~ 0.87 - 0.94, \U e2 \ ~ 0.35 - 0.49 (large mixing MSW) , (6.27) 

\U e i\ ~ 0.71 - 0.88 , \U e2 \ ~ 0.48 - 0.71 (vacuum oscillations) . (6.28) 

Notice that these ranges are statistically rather stable. For example, the range of 

sin 2 2-# sun allowed at 99% CL in the case of the large mixing angle MSW solution is 

0.36 < sin 2 2$ sun < 0.85 which implies \U el \ ~ 0.83 - 0.95, \U e2 \ ~ 0.32 - 0.55 
(compare with Eq.( |6.27D ). 



6.1.2 Atmospheric neutrinos 

The atmospheric neutrino data of the Super-Kamiokande experiment give an additional 
indication that \U e3 \ 2 is small. Indeed, the up-down asymmetry of electron events A e that 
could be generated by — > v e or v e — > z/g (with /3^e) transitions is compatible with zero 
(see Section 5.1.3) §: 

A e = -0.036 ±0.067 ±0.02. (6.29) 
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In the three-neutrino scheme under consideration, the Super-Kamiokande atmospheric 
neutrino data can be explained by — > v T oscillations p|, p61| , |409| . The analysis of the 
data under the assumption of two-generation mixing yields the following constraints: 

= sin 2 2tf atm > 0.82 , 5 x 10" 4 eV 2 < Am 2 < 6 x 10" 3 eV 2 , (6.30) 

from which it follows with Eqs.( |6.7| ) and ( p.9|) that \U^ 3 \ 2 and \U r3 \ 2 must be close to 0.5 
and therefore \U e 3\ 2 small. 

Coming back to Eq.( |6.21| ) which implies that \U e s\ < 0.22 we see that the CHOOZ 
upper bound for \U^\ is not very small. It is suggestive, however, to assume that 

|f/ e3 |<l, (6.31) 
i.e. that the element \U e ^\ is so small that it is unimportant for neutrino oscillations 



407| . [410|| . Such an assumption has been recently adopted by several authors ||411 



Let us infer some consequences of the assumption ( |6.31|) . First of all, it is easy to 
see that if the inequality ( |6.31| ) is satisfied, matter effects are not essential for — > 
v T transitions of atmospheric and LBL neutrinos. Indeed, the evolution equation for 
neutrinos in matter in the flavour representation is (see Section 4) 

i-^-i/j = — (uAM 2 U ] + iW (6.32) 
ax 2E \ J 

with 

^ ee ^ a M j , AM 2 ee diag(0, Am 2 1; Am 2 31 ) , A = dmg(A cc , 0, 0) , (6.33) 

and 

A cc = 2EV CC , V cc = V2G F N e . (6.34) 
Remember that for anti-neutrinos the matter potential term Acc must be replaced by 
Acc — —Acc- 

The inequality ( |6.5| ) implies that the phase generated by Am^ can be neglected for 
atmospheric neutrinos and AM 2 can be approximated by 

AM 2 ~ diag(0, 0, Am 2 x ) . (6.35) 

In this case (taking into account that the phases of the matrix elements U a % can be 
included in the charged lepton fields) we have 

(UAM 2 rf) ai3 ~ Ar/4 \U a3 \ \Ups\ . (6.36) 

If the condition ( |6.31| ) is satisfied, the elements of the matrix U AM 2 U^ with indices 
a = e or/and (3 = e can be neglected. This means that the equation for the wave 
function a e is decoupled from those of a M , a T : 

• dcte Acc (a qv\ 

*te = -2E a " (6 - 37) 



. da a Amjj 
? 1L7 ' IE 



Uc&\ ^2\Up 3 \a p (a = /i,r). (6.38) 
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These equations imply that Pj^,'^ BL ^ = 1 and Pj^?^ BL ) is given by the two-generation 
expression for neutrino oscillations in vacuum, 

P^i BL) = \ sin 2 2tf atm (l - cos ^f^) (6-39) 

(see Eq.( |3.20| )), with Am atm = Am^ and the mixing angle $ a t m determined by 

sin 2 tf atm = \U^\ 2 , cos 2 tf atm = \U r3 \ 2 . (6.40) 



Therefore, if the condition ( |6.31| ) is satisfied, there are no ^5 v e oscillations in atmo- 



spheric and LBL neutrino experiments and matter effects are not important for z/ M — > v T 
transitions. This means that the two-generation analyses of the atmospheric neutrino 
data in terms of — > v T are appropriate in the three-neutrino scheme under con- 
sideration and yield information on the values of the parameters Am^ = Am 2 tm and 
|t^3| = sin$ atm . 



According to a recent analysis [IDS] of the atmospheric neutrino data, the ranges 
of Am atm and sin 2 2-# atm for — > v T oscillations allowed at 90% CL by the Super- 
Kamiokande multi-GeV data and by the data of atmospheric neutrino oscillation exper- 
iments are, respectively, 

4 x 10" 4 eV 2 < Am 2 tm < 8 x 10~ 3 eV 2 , 0.72 < sin 2 2tf atm < 1 , (6.41) 

4 x 10" 4 eV 2 < Ami < 6 x 10~ 3 eV 2 , 0.76 < sin 2 2tf atm < 1 . (6.42) 



Taking into account the comments after Eq.( |6.5|) , the ranges ( |6.42 ) of the two-generation 
mixing parameters can be used in order to constrain the values of |t^3| and \U T -&\ only 
in the case of the vacuum oscillation solution of the solar neutrino problem. If the solar 
neutrino problem is due to MSW transitions, we must use the less restrictive ranges 
( 5.41] ). Thus, assuming that \U^\ < \U t3 \, we have 



11^1^0.49-0.71, \U t3 \ ~ 0.71 - 0.87 (solar MSW) , (6.43) 
|C/ M3 | ~ 0.51 - 0.71 , \U t3 \ ~ 0.71 - 0.86 (solar vacuum oscillations) . (6.44) 



As in the case of the ranges ( |6.26|) -( |6.28|) , also the ranges (|6T43D C|6T44 D are statistically 
rather stable. For example, the range of sin 2 2$ sun allowed at 99% CL by all atmospheric 
neutrino data is 0.66 < sin 2 2$ sun < 1 ||409|| , which implies |{7 m3 | ~ 0.46 — 0.71, \U r3 \ ~ 
0.71 - 0.89 (compare with Eq.(ggg)). 



6.1.3 The mixing matrix in the case \U e3 \ 1 

We have seen that, if we assume that \U e3 \ is very small, the results of the two-generation 
analyses of the data of solar and atmospheric neutrino experiments allow to determine 
the moduli of the elements of the first line and of the third column of the mixing matrix 
U . The moduli of the other elements can be determined from the unitarity of U . The 
simplest way to do it is to start from the parameterization ( ^.2(J| ) of the 3x3 mixing 
matrix U. 
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A very small \U e s\ implies that | sx3 1 <C 1- Since in the parameterization ( |2.20| ) the 
CP-violating phase <5 13 is associated with s 13 , it follows that CP violation is negligible in 
the lepton sector (see also footnote 5) and we have 



(6.45) 



Using the information on | S12 1 — \U e 2\ and |s23| — | | given by the two-generation 
analyses of the results of solar and atmospheric neutrino experiments, we obtain for the 
moduli of the elements of the mixing matrix |[407|[ : 





( Ci 2 
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< 1 


U~ | 


~ s 12 023 
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Large mixing MSW: 
Vacuum oscillations: 
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0.94 
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-0.49 
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-0.71 
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-0.76 


0.51-0.71 


0.24- 


0.50 
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-0.62 


0.71-0.86 



(6.46) 
(6.47) 
(6.48) 



In the case of the small mixing angle MSW solution of the solar neutrino problem \U e z\ <C 
1 could be of the same order of magnitude as |t/ e2 |. 

The best fit of Super-Kamiokande data corresponds to maximal mixing, $23 = 7r/4 
and the vacuum oscillation solution of the solar neutrino problem includes the maximal 
mixing case $12 = vr/4. The combination of these two possibilities gives the bi-maximal 
mixing matrix 



/ 4s ± 0\ 



U 



V 



1 

5 
2 
1 
2 



1 

2 
1 
■ 2 



1 
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(6.49) 



Notice that bi-maximal mixing 



that has been assumed recently by several authors ||411 
is not compatible with the large mixing angle MSW solution of the solar neutrino problem 

EH- 

It is interesting to note that, because of the large mixing of and v T with z/2, the 
transitions of solar u e '8 in u^s and zz/s are of comparable magnitude. However, this 
phenomenon and the values of the entries in the (1/^, i/ T )-(i/i, v 2 ) sector of the mixing 
matrix cannot be checked with solar neutrino experiments because the low-energy z/^'s 
and v T 's coming from the sun can be detected only with neutral-current interactions, 
which are flavour-blind. Moreover, it will be very difficult to check the values of {U^l, 
\Ufi2\j \U T i\ and \U T 2\ in laboratory experiments because of the smallness of Am; 



2 

21- 



In the derivation of Eqs. fl6.46"l) -( |6.48| ) we have assumed that \U^\ < \U e i\ and \U t 



fi3 



< 



\U T s\. The other possibilities, \U e 2\ > |t/ e i| an d \Uu3\ > \U T 3\, are equivalent, respectively, 
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to an exchange of the first and second columns and to an exchange of the second and 
third rows in the matrices Q6.46 )— (|6.48 ). Unfortunately, these alternatives are hard to 



distinguish experimentally because of the above mentioned difficulty to measure directly 
the values of {U^l, {U^l, \U T \ \ and \U T 2\- Only the choice \U e 2\ < \U e \\, which is necessary 
for the MSW solutions of the solar neutrino problem, could be confirmed by the results 
of the new generation of solar neutrino experiments (Super-Kamiokande, SNO, ICARUS, 
Borexino, GNO and others ||353|| ) if they will allow to exclude the vacuum oscillation 
solution. 



6.1.4 Accelerator long-baseline experiments 



Future results from reactor long-baseline neutrino oscillation experiments (CHOOZ [|413 



Palo Verde |[288|| , Kam-Land ||290|| ) could allow to improve the upper bound ( |6.21| ) on 



\U e $\ . In this section we discuss how an improvement of this upper bound could be 
obtained by future accelerator long-baseline neutrino oscillation experiments that are 
sensitive to ^ -> v e transitions (K2K CT , MINOS p93| , ICARUS and others 



295| - ^98|1 ; for a review see Ref. 



If matter effects are not important, in the scheme under consideration the parameter 
Me measured in 



sin 2 2^^ measured in — > v e long-baseline experiments is given by (see Eq. (|6.9|) ) 



sin 2 2^ e = 4|f/ e3 | 2 |^3r. (6.50) 

If accelerator long-baseline neutrino oscillation experiments will not observe z/ M — > v e 
transitions and will place an upper bound sin 2 2-d^ e < sin 2 2-#j^ Lax ' ) , it will be possible to 
obtain the limit 

sin 2 2^r x) 



W < f , (6-51) 

4 l t7 M3|(min) 

where |t^i3| 2 min ) is the minimum value of \U^\ 2 allowed by the solution of the atmospheric 
neutrino anomaly and by the observation of —>■ v T long-baseline transitions. For 
example, taking |^3| 2 min ) = 0.25 (see Eq. ( |6.44j )) we have \U e ^\ 2 < sin 2 2i9^ ax ' ) . If a value 

of sin 2 2-$[ ( ™ ax ^ ~ 10~ 3 , corresponding to the sensitivity of the ICARUS experiment for 
one year of running ||294|| , will be reached, it will be possible to put the upper bound 
|f/ e3 | < 3 x 10- 2 . 

The observation of — > v T transitions in long-baseline experiments will allow to 
establish a lower bound for (L^l 2 because the parameter sin 2 2$^ T is given in the scheme 
under consideration by (see Eq.(|6.9j)) 

sin 2 2?V = 4|^ 3 | 2 |^ 3 | 2 . (6.52) 

From the unitarity relation \U e ^\ 2 + |f/ M 3| 2 + \U T s | 2 = 1 it follows that an experimental 
lower bound sin 2 2$ MT > sin 2 2-&^ m ' allows to constrain \U^\ 2 in the range 



2 



' ( ' / l-sin 2 2^r ) ) < |[/ M3 | 2 < \ ( 1 + y/l - sin 2 2^r } ) . (6.53) 
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It sm 2v, 



'fiT 



is found to be close to one, as suggested by the solution of the atmo- 



spheric neutrino problem (see Eqs. fl6.41p and fl6.42|) ), the lower bound |t^i3^ min ) 



sin 2 24r in) ^ is close to 1/2. 



If matter effects are important, the extraction of an upper bound for |?7 e3 | 2 from the 
accelerator long-baseline experiments is more complicated. In this case 



data of — 

the probability of z/ M — > v e oscillations is given by ||192 | 



4|[/ e3 l 2 lt/, 



/ / 3 1 



sm 



2 Acc 
Am 31 



Am^L 
AE 



1 - 



Acc 

Am ll 



+ 4|C/ e3 | 2 ^ 

I ej l Ami, 



(6.54) 

where E is the neutrino energy and L is the distance of propagation. This probability 
depends on the neutrino energy not only through the explicit E in the denominator 
of the phase, but also through the energy dependence of A cc = 2EVcc- For long- 
baseline neutrino beams propagating in the lithosphere of the earth the charged-current 
effective potential Vcc — \f2GpN e is practically constant: N e ~ 1.5iV avo cm -3 and 
V cc ^ 1-1 x 10" 13 eV (see Eq.(F37p). 



If long-baseline experiments will not observe 



transitions (or will find that 



they have an extremely small probability) for neutrino energies such that Acc A m !i> 
it will mean that \U e s\ 2 is small and a fit of the experimental data with the formula 
( |6.54| ) will yield a stringent upper limit for \U e 3\ 2 (taking into account the lower limit 
\U^3\ 2 > |£^i3| 2 min ) obtained from the solution of the atmospheric neutrino anomaly and 
from the observation of — » v T long-baseline transitions). On the other hand, the non- 
observation of — > v e transitions for neutrino energies such that Acc ^ ^ m \i does not 



provide any information on | C/ e 3 1 because in this case the transition probability ( |6.54j ) is 



suppressed by the matter effect. Hence, in order to check the hypothesis \U e %\ <C 1, as 
well as to have some possibility to observe — > u e transitions if this hypothesis is wrong, 
it is necessary that a substantial part of the energy spectrum of the neutrino beam lies 
below 



Am 2 x 



2V< 



30 GeV 



cc 



A ^31 

10- 2 eV 2 



(6.55) 



This requirement will be satisfied in the accelerator long-baseline experiments under 
preparation (K2K J27IJ, MINOS ICARUS and others P^hM ; see |TTJ) if 



Am^ is not much smaller than 10 2 eV 2 . 



6.2 Mixing of four massive neutrinos 



Let us consider now the case of mixing of four neutrinos [ 415 -417], i.e., n = 4 in Eqs.(p^3] 
and Q, 



| "a) = y^^qfc Wk) , Wa) = ^2U ak \v k ) (a = 6,fl,T,s) 



k=l 



k=l 



(6.56) 
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Figure 6.1. The six types of neutrino mass spectra that can accommodate the solar, atmospheric and 



LSND scales of Am 2 . 



which allows to accommodate all the three existing indications in favour of neutrino 
mixing. Here \v s ) is a sterile neutrino state (see Section 2.2). 

We do not make any assumptions about the mass spectrum of four neutrinos and will 
consider all possibilities. The six possible four-neutrino mass spectra with three different 
scales of neutrino mass-squared differences (atmospheric, solar and LSND) are shown in 
Fig. |5TT| . In all these mass spectra there are two groups of close masses separated by 
the "LSND gap" of the order of 1 eV. The small mass-squared differences correspond 
to Am;? un (the smallest one, Am 21 in schemes I and B, Am§ 2 in schemes II and IV, 
Am.4 3 in schemes III and A) and Am^ tm (Am| x in schemes I and II, Am^ in schemes III 
and IV, Am 21 in scheme A, Am| 3 in scheme B) and the largest mass squared difference 
Amjj = Am^sND i s relevant for the oscillations observed in the LSND experiment. 

Let us consider SBL oscillations in the case of neutrino mass spectra presented in 



Fig. p.l| . Taking into account that 



/Ami L 

^— < 1 

E 



and 



Am^ tm L 
E 



< 1 



from the general formula ( 3.12j ) we obtain that the probabilities of v a 
and of v a survival in SBL experiments are given by 



p(SBL) 



p(SBL) 



2^ 



cos • 



Am^ L 



IE 



1 — - B a . a ( 1 



cos • 



Am^ L 
2E 



(6.57) 
up transitions 

(6.58) 
(6.59) 



The oscillation amplitudes A a .p and B a . a depend on the elements of the mixing matrix 
and on the type of neutrino mass spectrum: 

2 



a;/3 



Sk 




(6.60) 



(6.61) 
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where the index k runs over the indices of the first or (because of the unitarity of the 
mixing matrix) of the second group of neutrino masses (see Fig. |6.1| ). The expressions 
( |6.5S ) and (|6.59|) have the same form as the corresponding two-neutrino formulas ( |3.20| ) 
and (|3.21|) and are a direct generalization of the expressions (|6.7|) and ( |6.8|) presented in 
Section |0 for the case of three massive neutrinos and a neutrino mass hierarchy. From 



Eq. ( |6.60| ) it is obvious that the amplitude of v a — ► vp oscillations is equal to A a .p and 



,(SBL) 



p(SBL) 



(6.62) 



Thus, CP violation in the lepton sector cannot be revealed in SBL neutrino oscillation 
experiments if only one mass-squared difference is relevant for SBL oscillations. 

The expression (|6.61|) allows to obtain information on the quantities Ylk \U e k\ 2 and 
J2k \ U[ik\ 2 from the results of SBL v e and disappearance experiments, whose exclusion 
curves imply the constraint 



B a - a < B 



[a 



e,/i) . 



(6.63) 



The upper bound B®. a depends on the value of Am^ 



that each Y^ fc |£/o,fc| 2 for a 



E 



Equations (|6.59|) and (|6.63| ) imply 
e, \i must satisfy one of the two inequalities 



U ak \ 2 <a° n 



or 



E 



u. 



ak 



> 1 



(6.64) 



with 



5" 



1 - B° 



(6.65) 



The most stringent limits on B e]e and B^ are given by the exclusion curves of the 
Bugey [|162|1 reactor v e — > v e experiment and of the CDHS |1 53| and CCFR R4181 accelera- 
tor v„ — > v a experiments. The values of a° e and a° corresponding to the 90% CL exclusion 

]. One can see that 



curves of these experiments are shown in Fig. 



a° e < 4 x 10 
a° < 2 x 10 



2 
-1 



for 
for 



^O.leV 2 , 
Am 2 ! > 0.3 eV 2 . 



(6.66) 
(6.67) 



Hence, from Eq.( |6.64| ) it follows that the quantities J2k \ Uek\ 2 and {U^l 2 are either 
small or large (close to one). We will show now that from the four possibilities for the 
two quantities \ U e k\ 2 and |t/ M fc| 2 (small-small, small-large, large-small, large-large) 



for each neutrino mass spectrum in Fig. |6.1| only one possibility is compatible with the 
results of solar and atmospheric neutrino experiments [[416| . |119| . |420| . 
Let us consider, for example, the hierarchical mass spectrum 



mi < m 2 « m 3 m 4 , (6.68) 
which corresponds to type I in Fig. O. In this case, the inequalities (|6.64j) become 



\U a 4 2 <a° a or \U ai \ 2 >l-a° a (a = e,/i). (6.69) 
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Figure |6.' 

as functions of Amfi 



Values of a? e and a° (see Eq.( 3.65 )) 
obtained from the exclusion 
and CCFR 



curves of the Bugey |6g], CDHS |1S 
[ 418 1 SBL disappearance experiments 



> 
u 



g 




Figure 6.3. The region on the right of the solid 
line passing through the black circles is excluded 
at 90% CL by the results of SBL disappearance 
experiments for small \U e 4\ 2 and {U^l 2 in the 
model with mixing of four neutrinos and a mass 
hierarchy. The shadowed regions are allowed at 
90% CL by the results of the LSND experiment. 
Also shown are the exclusion curves of the Bugey 
1161, BNL E734 §7|, BNL E776 p7f , CCFR 



[ 374( | and KARMEN (Bayesian analysis) |7§ ex- 
periments. 



The survival probability of solar z/ e 's averaged over the fast, unobservable oscillations 
due to Am^ and Am^ is given by [[421] 



pt»l = ( i - £ i^i 2 ) p £l + E i^i 4 > ( 6 - 7 °) 

V fc=3,4 / fc=3,4 

(1 2) 

where Pv e ^-Ju e is the v e survival probability due to the mixing of v e with v\ and z/ 2 . 
This probability has the two-generation form (|3.21|) and depends on the neutrino energy, 



Am 2 = Ami„ = Am^ and 



4 




2 


u e2 


2 


{Wei 


2 + 


u e2 


2\2 



sin 2 2tf = sin 2 2tf sun = 1 f e2 . (6.71) 



Though Eqs.( |b.7C| ) and ( |6.71| ) were derived here for vacuum oscillations, the considera- 
tions in Section 4.4 show that they are also valid including matter effects. If |?7 e 4| 2 > 1— aP e , 
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from Eqs. fl6.66"D and ( |6.70| ) one can see that Pj^l e practically does not depend on the 
neutrino energy and Pjf™2 e <L 0.92 for all solar neutrino energies. This is not compatible 
with the solar neutrino data [ 329 , 349| , |422|| . Therefore, we reach the conclusion that from 
the two possibilities ( |6.69|) for \U e ^\ 2 only 



\U e4 \ 2 < a° e 



(6.72) 



is compatible with the results of solar neutrino experiments. 

Let us now consider the survival probability of atmospheric i^'s. From the general 
formula ( |3.12j ) we obtain 



p>(atm) 



£im 2 +i^iW-;^) 

fe=l,2 ^ ' 



+ \U, 12 



ha\ exp -i 



2E 



(6.73) 



where we have taken into account that for atmospheric neutrinos^] Am^iL/E <C 1. The 
probability averaged over the fast unobservable oscillations due to Am^ is given by 



= ( E i^i 2 ) +i^i 4 +2 f e i^i 2 ) 

\fc=l,2 / \jfc=l,2 / 



u, 



fi 3 I 



COS 



Am^I 
2£ 



+117, 



(6.74) 



Using the unitarity relation 



Ei^i 



1 - \u, 



/i4| 5 



(6.75) 



fc=i 



the survival probability (|6.74|) can be written as ||416| , [420| (see also Section 4.4) 

P$£l = (1 " MfP^l + IM* , (6.76) 

where Pi^-i Vjl is the survival probability of atmospheric u^ s due to the mixing of with 
U3 and V2, vi- This probability has the two- generation form 

1 



p(l,2;3) 



~ sin 2 



cos • 



Am^L 
2E 



(6.77) 



(see Eq.(|3~2il)), with 



sin 2 2$ 



atm 



4 


^3 


2 (1^1 


2 + 




2 ) 


(l^i 


2 


+ 




2 + 




2^ 2 



(6.78) 



The zenith angle dependence of /i-like events in atmospheric neutrino experiments is due 



to 

Pvp^vm which is multiplied by (1 — |t/^4| 2 ) in the expression 

am for p^i. if 



2\ 2 



)(atm) 



27 As explained after Eq.( |6.5| ), in the case of the MSW solutions of the solar neutrino problem this 
inequality and the following equations are valid only for multi-GeV neutrinos. 
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I £7^4 1 2 is large, as given by the second inequality in Eq. ( |6.69| ), |t/ M 4| 2 > 1 — a°, we have 
(1 - \U^ A \ 2 f < 4 x 10~ 2 for Am^ > 0.3 eV 2 (see Eq. (|B37D ). A zenith-angle variation 
of //-like events smaller than 4 x 10~ 2 is incompatible with the up-down asymmetry of 
multi-GeV /i-like events observed by the Super-Kamiokande experiment (see Section [5.1| ). 
Therefore, from the two inequalities ( |6.69| ) for |£7 M 4| 2 only the first one is compatible with 



the results of the atmospheric neutrino experiments: 

|^ 4 | 2 <a° (6.79) 



The same conclusion can be reached for all the neutrino mass spectra in Fig. |6.1| in 
which one mass is separated by the other three by the LSND gap, i.e., for the spectra 
I, II, III, IV. Hence, in the framework of these spectra the results of SBL disappearance 
experiments and of solar and atmospheric neutrino experiments imply that 

\U ej \ 2 <a° e , |f/ Mi | 2 <aJ, (6.80) 

with j = 4 in the schemes I and II and j — 1 in the schemes III and IV. These inequalities 
imply that i^— >f e transitions in SBL experiments are strongly suppressed. Indeed, from 
Eqs.( |6.60| ) and (|6.80|) we have 

^ ;e = 4|f/ e /|f^| 2 <4a°aJ. (6.81) 

Thus, the amplitude of SBL v^— >v e oscillations is quadratic in the small quantities \U e j\ 2 
and \U^j\ 2 . The upper bound ( |6.81| ) obtained from the values of a° e and a° presented 
in Fig. |5.2| is shown in Fig. |6.3| by the solid line passing through the black circles. The 



solid, dotted, dashed, dash-dotted and dash-dot-dotted curves in Fig. IO are the 90% 



CL exclusion curves obtained in the Bugey M , BNL E734 p72| , BNL E776 |373|, 



CCFR [ 374j| and KARMEN (Bayesian analysis) ||375|| experiments. The shadowed regions 



are allowed by the results of the LSND |370| , |371|1 experiments at 90% CL. From this figure 



one can see that the result of the LSND experiment, taken together with the results of SBL 
disappearance experiments and solar and atmospheric neutrino experiments, disfavours 
the schemes I, II, III and IV of Fig. |Q|-H 

Let us consider now the spectra A and B in Fig. |6.1| . In this case the inequalities 
( |6.64| ) apply with the index k running over 1,2 or 3,4 (the two choices are equivalent 
because of the unitarity of the mixing matrix). We consider explicitly only the spectrum 
A, but all the following results are valid also for scheme B and can be obtained with the 
index exchange 1,2 ^ 3, 4. 

In the case of spectrum A, the survival probability of solar z/ e 's is given by ||421 



PtX = E l^i 4 + 1 " E \ U *\* PllX ■ (6-82) 



fc=l,2 \ fc=l,2 




The presence of a small allowed region in Fig. 6.3 for 0.2 eV 2 < Am^ < 0.3 eV 2 is due to the fact that 
for Am 2 ! < 0.3 eV 2 there are no constraints on B®.^ from the results of disappearance experiments. 
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where Pj£$u e is the survival probability of solar u e 's due to the mixing of v e with u 3 and 



z/ 4 . This probability has the two-generation form ( |3.21|) and depends on the neutrino 
energy, Am 2 = Am 2 u 



Am.43 anc I 



sin 2 2d = sin 2 2tf, 



4 


^e3 


to 


£/"e4 


2 


(1^3 


2 + 


t/e4 


2^) 2 



f/ e fc| 2 > 1 — a°, i.e., the second inequality ( |6.64| ) for a 



If £,=1,2 

is satisfied, from Eq. (|6.66 ) we have 



(6.83) 



e and k — 1,2 



\U, 



< 1.6 x 10 3 , which implies 



£fc=l,2 \ u ek\ 

that Pvl^le — X^fc=i 2 l^efcl 4 is constant and Pvl^l e > 0.5. Since these constraints are 
incompatible with the solar neutrino data, we reach the conclusion that only the first 
inequality (|6.64j) for a = e and k = 1, 2 is compatible with the results of solar neutrino 
experiments, i.e., £ fc=1 2 \ U e k\ 2 < oP e . 

Now we consider the survival probability of atmospheric z/ M 's. In scheme A it is given 
by |42l 



p(atm) 



E 1^1 

fc=3,4 



\fc=3,4 / 



(6.84) 



(1 2") 

where the probability Pv^'-Jv^ has the two-generation form ( |3.21| ) with A 



m 



Am ltm 



Am^ and 



sin 2 2d = sin 2 2d 



at m 



4 




to 




2 




2 + 




2^2 



(6.85) 



\U, 



tj,k\ > 1 _ which corresponds to the second of the two inequalities ( |6.64|) 



3,4, the upper bound (|6.67|) for a implies that Puf™l^ is practi 



^ £,=3,4 

for a — n and k 

cally constant and incompatible with the zenith-angle dependence of multi-GeV /i-like 
events observed in the Super-Kamiokande experiment ||. Hence, only the first inequality 
( p. 64 ) for a — \i and k = 3,4 is compatible with the results of the Super-Kamiokande 



U, 



atmospheric neutrino experiments and we have £ fc=3 4 

Summarizing, in the framework of scheme A the results of SBL disappearance exper- 
iments and of solar and atmospheric neutrino experiments imply that 



E i^i 2 ^ a ° 



and El^l 2 ^! ( A )' 



(6.86) 



fc=l,2 fc=3,4 

The corresponding bounds in scheme B are obtained with the exchange 1, 2 

E l^fel 2 < a " aIld E M < a l (B) • 



3,4: 



(6.87) 



fc=3,4 



fc=l,2 



The constraints (|6.86| ) and ( |6.87| ) imply the following upper bound for the amplitude 



of — > v e transitions in SBL experiments: 



E u ^ u u 



< 



4 (?H(? 



U, k \ 2 <4Min(a>^ 



(6.88) 
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(here the index k runs over 1, 2 or over 3, 4). This bound is linear in the small quantities 
<2g and a° and turns out to be compatible with the results of the LSND experiment. 
Actually, this bound is worse than the unitarity bound 



A Ke < Mm(B° e;e , flj. ) 



(6.89) 



which causes the exclusion of the large sin 2 2t? part of the LSND allowed region (see 
Fig. |BT3| , where B®. e is represented by the Bugey curve). 



We want to note that the bound to exist shortly below Am^ = 0.3 eV 2 

(see Fig. |6.2|) . On the other hand, it follows from the LSND and Bugey experiments that 
Am^ > 0.2 eV 2 (see Fig. Thus, for the small interval 0.2 eV 2 < Am 2 x < 0.3 eV 2 

our arguments leading to schemes A and B have to be considered with caution and it 
would be desirable to have experimental information on B KIX in this interval. 

Let us now discuss some consequences of the schemes A and B, taking into account 
the constraints (|6.86|) and ( |6.87| ). First of all, we will show that these bounds lead to 
constraints on the neutrino mass measured in tritium /3-decay experiments and on the 
effective Majorana mass that determines the probability of neutrinoless double-/? decay. 

The electron spectrum in the decay 3 H — > 3 He + e~ + v e is given by 



^ = Cp e E e (Q-T)F(E e 



J2\U ei \y(Q-T) 



(6.90) 



where Q is the energy release, p e and E e are the electron momentum and energy, T = 
E e — m e , F(E e ) is the Fermi function and C is a constant. Taking into account that 
mi — m2 "C — ttl^, the first inequality in Eq. ( |6.86D implies that in scheme A the 
effective neutrino mass determined from the measurement of the high-energy part of the 
/3-spectrum of 3 H is practically equal to the "LSND mass" m 4 : 



m( 



TL) ~ 1^ 



k \ 2 m A 



k=3A 



I - > \U e k\" I "' I - m i ■ 
fc=l,2 



(6.91) 



If massive neutrinos are Majorana particles, the matrix element of neutrinoless double- 
P decay is proportional to the effective Majorana mass [|95l , |106|| (see Appendix A) 



ek m k 



k=l 



In the case of scheme A, the first inequality in Eq.( |6.86D implies that 



(6.92) 



I Ml 



k=3A 



m 4 ~ m 4 \ 1 — 4 |[/ e4 | 2 (1 — \U eA \ 2 ) sin z 



(6.93) 



where is the difference of the phases of U e s and U e ±. Since is unknown, from Eq. ( |6.93| ) 
we obtain the constraints 



|l-2|f/ e4 | 2 |m 4 < |(m) I <m 4 



(6.94) 
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where the upper and lower bounds correspond to = or n and 
Taking into account Eqs. fl6.83p and fl6.86|) , we have 



\U, 



e3 | 



COSl9 s 



and 



\U, 



el 



and Eq.( |6.94[ ) can be written as 



m 4 Vl - sin 2 27? aun < \ (m) \ < 



sin??.. 



m 4 . 



±7r/2, respectively.^ 



(6.95) 



(6.96) 



' and of the conservation of CP. In the 
sun , which correspond to the MSW solution 



If sin 2 2$ sun <C 1, which corresponds to the MSW solution with a small mixing angle, we 
have \(m)\ ~ m 4 , independent of the value of 
case of a large value of the parameter sin 2 2$ 
with a large mixing angle or to the vacuum oscillation solution, in the future it will be 
possible to obtain information about the violation of CP in the lepton sector if both 
the 3 H /3-decay and neutrinoless double-/? decay experiments will obtain positive results. 
Indeed, from the measurements of m( 4 H) ~ m 4 , \(m)\ and sin 2 2$ sun with the help of 
Eq. Q6.96p it will be possible to determine the value of sin 2 (if sin 2 2$ sun is large). Let us 
emphasize that in the case of CP conservation the relative CP parities of and z/ 4 can 



be determined: the CP parities are equal if sin z 



and opposite if sin 



Considering now scheme B, the first inequality in Eq.( |6.87| ) implies that both m( H) 
and | (m) | are strongly suppressed: 



m 



w. 



ek 



uii < a e m 4 -C m 4 . 



\k=3A 



IWI 



fc=3,4 



ek 



m 4 < 




m 4 < a m 4 <C m 4 



(6.97) 



(6.98) 



Therefore, only if scheme A is realized in nature, there is some possibility to measure the 
effect of the LSND mass in future tritium decay experiments and neutrinoless double-/^ 
decay experiments. 

We discuss now the implications of the schemes A and B for LBL neutrino oscillation 



The smallness of J2 k=12 \ U e k\ 2 in scheme A and of J2k=3 4 Wek 



m 



experiments 

scheme B implies that the electron neutrino has a small mixing with the neutrinos whose 
mass-squared difference is responsible for the oscillations of atmospheric neutrinos (i.e., 
Vi, v<i in scheme A and Us, z/ 4 in scheme B). Hence, as shown in the following part of 
this section, the transition probability of electron neutrinos and antineutrinos into other 
states in atmospheric and LBL experiments is suppressed. 

Taking into account that in scheme A we have Am^L/i? <C 1 and Am^L/E 3> 1 in 
LBL experiments, we obtain from the general formula (|3.12|) the following expression for 



29 If CP is conserved, we have (f> = or it if the CP p aritie s of V3 and 1^4 are equal and <fi 
they are opposite. Indeed, if CP is cons erved from Eq.( 2.5C ) we get U* k = U e kPk where p k 



the CP parity of v k is r\ k — ipk (see Eq.( 2.53 )) 
if Pk = +1 and 
if P3 = -pi- 



= ±7r/2 if pk = — 1. Hence, we have 



Writing U e k = \U e k\e l<pk , we obtain 
t>3 — 4>A = or 7r if pz = Pa and <fo — 



= ±tt/2 if 
= ±1 and 
k = or 7r 
64 = ±7r/2 
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the measurable probability of u a — > up transitions averaged over the fast oscillations due 



to Am^: 



p(LBL;A) 



Up 1 U* al + Up 2 U* a2 exp - 



IS.m\ x L 
2E 



£ 

fc=3,4 



(6.99) 



3, 4 and the 



The corresponding probability in scheme B is given by the exchange 1, 2 
probability of u a — > up transitions is obtained changing Upk U* k — > U% k U a k- 

Let us consider reactor LBL experiments which measure the survival probability 
Pp^y . We will show that the negative results of SBL reactor experiments and the results 
of solar neutrino experiments imply strong constraints on the transition probability of 
LBL reactor z/ e 's into other states. Indeed, from Eq. ( |6.99|) and the corresponding one in 
scheme B, for the survival probability of LBL reactor z/ e 's we have the lower bounds 



(6.100) 
(6.101) 




in schemes A and B, respectively. From the inequalities ( |6.86 ) and ( 6.87 ) it follows that 
in both schemes 

P£S>(l-a°) 2 ■ (6.102) 

Since a° e is small (see Fig. |0| and Eq. ( |6.66| )), the lower bound ( J6.102| ) implies that the 
survival probability Pp^& e is close to one. For the transition probability of u e 's into any 
other state in LBL experiments, 



Pi 



(LBL) 



P^e 



(LBL) 

Up > 



we have the upper bound 



1-P 



(LBL) 



< a° e (2 - al) 



(6.103) 



(6.104) 



The value of this bound, which depends on Am| x through the dependence of a° e on Am^ 
(see Fig. |5T^ ), is shown in Fig |T3| (the solid line), where it is compared with the upper 
bound for 1 — P^^j, obtained in the CHOOZ experiment ||159|| (dash-dotted line) and 
with the final sensitivity of the CHOOZ experiment (dash-dot-dotted line). The shadowed 
region in Fig ICTl corresponds to the range of Am^ allowed at 90% CL by the results 



of the LSND and other experiments. One can see from Fig 5.4 that if Am 4 2 , lies in the 



LSND range, the upper bound for 1 — Pj£^p e is much less than the sensitivity of the 
CHOOZ experiment and it will be very difficult to measure the disappearance of reactor 
u e in future LBL experiments. 

u e transitions, which will be investigated in the near future 
ICARUS g9§] and/or other p5H 298| experiments 



Let us consider now u^ — 
MINOS 



by the K2K [271 



293] 
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Figure ^J. The value of the bound ( 13.104 ) on 
the transition probability of D e into all other neu- 
trino types in LBL experiments as a function of 
Am| x (solid line). The dash-dotted and dash- 
dot-dotted lines represent the upper bound ob- 
tained in the CHOOZ experiment [159 and the 
final sensitivity planned in this experiment, re- 
spectively. The shadowed region corresponds to 
the range of Am^ allowed at 90% CL by the re- 
sults of LSND and other experiments. 
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Figure 3.5 
ity of 

as functions of Amfi 



Upper bounds on the probabil- 
Vf, transitions in LBL experiments 
The dashed line repre- 



sents the bound (6.109), whereas the dash-dotted 
and dash-dot-dotted curves include matter cor- 
rections to this bound for the K2K [L = 250 km) 
and MINOS and ICARUS (L = 730 km) exper- 
iments, respectively. The vertical lines show the 
planned sensitivities of these experiments. The 
dotted curve depicts the bound (6.111). The 
range of Am^ allowed by LSND and other ex- 
periments is also indicated. 



sec 



414|| ). With the help of the Cauchy— Schwarz inequality, from Eq. (|6.99|) and Eq.( |6.60|) 



we obtain 



< (e (e i^i 



_ A 

4 A l*;e > 



(6.105) 



where the index k runs over 1, 2 in the scheme A and over 3, 4 in the scheme B. From the 
constraints (|6.86|) and (|6.87|) , in both schemes A and B we have 



p(LBL) < o , 1 A o 



(6.106) 



where v4°. e is the upper bound for the oscillation amplitude A^. e measured in SBL exper- 
iments. 

Since the CPT relation ( |3.15| ) and the conservation of probability, relation ( |6.103| ), 
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imply that 



,(LBL) 

!/ e — *V<L 1 



(6.107) 



Ep(LBL) p(LBL) _ -I _ p_( 

another upper limit for Pj^k) can be obtained from the inequality (|6.104|) : 

PS^X <al(2- al) . (6.108) 
Combining the upper bounds (|6.106|) and ( |6.108 ) and taking into account that they are 



valid also for antineutrinos 423 , we have 



^(lbl) < Min ^ fl 



1 



(-) (-) 



+ 4 K\e > 



A (2 - al) 



(6.109) 



The curve corresponding to this limit obtained from the 90% CL exclusion plots of the 
Bugey pi experiment for a° and of the BNL E734 g7|, BNL E776 |373|, CCFR |37l 
and KARMEN (Bayesian analysis) ||375|| experiments for A°. e is represented by the dashed 
line in Fig. |6.5| . For comparison, the expected sensitivities of the LBL accelerator neutrino 
experiments K2K ^71)] , MINOS |2U3J and ICARUS fH§ are also indicated (the three 
straight vertical lines). The dash-dotted and dash-dot-dotted lines in Fig. |6T5| represent, 



respectively, the upper bound for P^^X in the K2K and MINOS, ICARUS experiments 
corrected for the matter effects [423 due to the propagation of the neutrino beam in 
the crust of the earth. These effects are different for the K2K and MINOS, ICARUS 
experiments because of the different value of L in this experiments: L ~ 250 km in K2K 
and L ~ 730 km in MINOS and ICARUS. (It is obvious that the bound for other LBL 
experiments with a neutrino beam from CERN to Gran Sasso is the same as that for the 
ICARUS experiment.) 



The region between the two horizontal solid lines in Fig. |6.5| corresponds to the range 
of Am^j allowed at 90% CL by the results of the LSND and other experiments. One 
can see that if Am^ is in the LSND region, the sensitivity of the K2K experiment is not 
sufficient for the observation of u e transitions, whereas experiments with higher 

sensitivities like MINOS and ICARUS could reveal these transitions. 

Finally, another upper bound on B^ B ^\ can be obtained from Eq.( |6.106l ) noting that 



U e k\ 2 



(6.110) 



with k — 1, 2 in scheme A and k = 3, 4 in scheme B. Taking into account the inequalities 
( |6.86l ) an d Q6.871) , we obtain 



(6.111) 



This bound is stable against matter effects [423] and its value is represented by the dotted 
curve in Fig. |B~5] . For a <C aP e 1 this bound is about half of that given by Eq. ( 16.108]) . 
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Figure 3.6. The two types of neutrino mass spectra that can accommodate the solar, atmospheric and 



LSND scales of Am and the mixing schemes that emerge if the Big-Bang Nucleosynthesis constraint on 
the number of light neutrinos is less than 4. 



However, since a° is only small in the same range of Am 2 where A^. e is small, numerically 
the bound ( |6.111| ) turns out to be worse than the bound ( |6.109D (dashed line in Fig. 



and the corresponding matter-corrected bound (dash-dotted and dash-dot-dotted lines in 
Fig. £3). 

Summarizing the results presented in this section, we have shown that only the two 
four-neutrino schemes A and B are compatible with the results of all neutrino oscillation 
experiments. If scheme A is realized in nature, the effect of non-zero neutrino masses 
could be observed in the near future in tritium /3-decay experiments and neutrinoless 
double-/? decay experiments, whereas these effects are strongly suppressed in scheme B. 
Furthermore, we have shown that the results of SBL v e disappearance experiments and 
solar neutrino experiments imply that in both schemes A and B the electron neutrino 
has a small mixing with the neutrinos whose mass-squared difference is responsible for 
the oscillations of atmospheric neutrinos (i.e., V\, v 2 in scheme A and u 3 , z/ 4 in scheme 
B). As a consequence, the transition probability of electron neutrinos and antineutrinos 
into other states in atmospheric and LBL experiments is strongly suppressed. 

Furthermore, it has been shown in Refs. |[417| , f424|| that in the schemes A and B, if the 



standard Big-Bang Nucleosynthesis constraint on the number of light neutrinos (see ||425|| 
and references therein) is less than 4, then one has a stringent limit on the mixing of the 
sterile neutrino with the two massive neutrinos that are responsible for the oscillations of 
atmospheric neutrinos. In this case, v s is mainly mixed with the two massive neutrinos 
that contribute to solar neutrino oscillations (z/3 and z/4 in scheme A and v\ and z/2 in 
scheme B) and the unitarity of the mixing matrix implies that v T is mainly mixed with 
the two massive neutrinos that contribute to the oscillations of atmospheric neutrinos. 
Hence, the two schemes have the form shown in Fig. |6.6| and have the following testable 
implications for solar, atmospheric, long-baseline and short-baseline neutrino oscillation 
experiments: 

• The solar neutrino problem is due to v e — > v s oscillations. This prediction will be 
checked by future solar neutrino experiments that can measure the ratio of neutral- 



current and charged-current events ||356|| 
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The atmospheric neutrino anomaly is due to Vn 



oscillations. This prediction 



will be investigated by LBL experiments. Furthermore, the absence of — > v s atmo- 
spheric neutrino oscillations may be checked in the future in the Super-Kamiokande 
atmospheric neutrino experiment [|2 7 1| , |2 72| , [2 75| , |2 76|| . 



and v e — > v s oscillations are strongly suppressed in SBL experiments [[424 



If these prediction will be falsified by future experiments it could mean that some of 
the indications given by the results of neutrino oscillations experiments are wrong and 
neither of the two four neutrino schemes A and B is realized in nature, or that Big-Bang 
Nucleosynthesis occurs with a non-standard mechanism |426|1 . On the other hand, if the 
standard Big-Bang Nucleosynthesis constraint on the number of light neutrinos is less 
than 4 and one of the two four neutrino schemes depicted in Fig. |6.6| is realized in nature, 
at the leading order the 4x4 neutrino mixing matrix has an extremely simple structure 



in which the v e 



/a 



and Up, v T 



sectors are decoupled: in scheme A 



U ~ 



and in scheme B 
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(6.112) 



(6.113) 



where ^sun and t? atm are, respectively, the two-generation mixing angles relevant in solar 
and atmospheric neutrino oscillations. Therefore, the oscillations of solar and atmospheric 
neutrinos are independent and the two-generation analyses of solar and atmospheric neu- 
trino oscillations yield correct information on the mixing of four neutrinos. The ^v^-^Ve 
transitions observed in the LSND experiment represent only a slight distortion of this 
simple picture. 



7 Conclusions 

In this review, we have discussed the phenomenological aspects of neutrino oscillations. 
The phenomenon of neutrino oscillations is possible if neutrinos are massive and mixed 
particles. The problem of neutrino masses and mixing is the central problem of today's 
neutrino physics and astrophysics. More than 40 different experiments all over the world 
are dedicated to the investigation of this problem and several new experiments are in 
preparation. 

The investigation of the properties of neutrinos is considered as one of the most 
important direction for the search of new physics (see Ref. |31|). Massive neutrinos are 
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among the plausible candidates for dark matter particles and the number of massive light 
neutrinos plays crucial role in Big-Bang Nucleosynthesis |[425| . 

The Super-Kamiokande collaboration obtained recently a rather convincing evidence 
in favour of transitions of atmospheric into other neutrino states iv T or z/ ste riie)- Other 
indications in favour of neutrino mixing were obtained in all solar neutrino experiments 
and in the accelerator LSND experiment. 

If all these indications will be confirmed by future experiments, it will mean that the 
number of massive light neutrinos is larger than the number of flavour neutrinos, i.e. 
that sterile neutrinos exist. This would be an evidence that neutrino masses and mixing 
are phenomena due to physics beyond the Standard Model with right-handed neutrino 
singlets. 

If the indications in favour of a "large" Am 2 ~ 1 eV 2 obtained in the LSND experi- 
ment will not be confirmed by future experiments, the mixing of three massive neutrinos 
with a mass hierarchy is a possible and plausible scenario. In this case, the Dirac or 
Majorana nature of massive neutrinos can be determined only increasing the sensitivity 
of the experiments on the search for neutrinoless double-/3 decay by almost two orders of 
magnitude. It has been shown by an analysis of the existing data |[427| , |428|1 that the effec- 
tive Majorana mass in the case of three massive Majorana neutrinos and a mass hierarchy 
is smaller than 10~ 2 eV, whereas the present upper bound is about 0.5 eV ||429| , |430|| and 
the sensitivity of the next generation of experiments will be of the order of 10 _1 eV |431 



Future neutrino oscillation experiments will make possible to investigate in detail the 
regions of Am 2 and of the elements of the neutrino mixing matrix allowed by the existing 
results of neutrino oscillation experiments. Future solar neutrino experiments (SNO, 
ICARUS, Borexino, GNO and others ||353|| ) will provide model- independent information 
on the transition of high energy 8 B neutrinos into other states and investigate in detail 
the fluxes of medium energy 7 Be, pep and CNO neutrinos. Future accelerator LBL 
experiments (K2K |71, MINOS |93|, ICARUS |29| and others Iggg - gg^gH ) will 
explore the region of Am 2 indicated by the atmospheric neutrino anomaly. Besides 
v„ — 



survival, also za, 



and Un — > u e transitions will be studied. A very important 



task for the neutrino oscillation experiments of the next generation is the check of the 
indication in favour of neutrino oscillations obtained in the LSND experiment. The 
KARMEN experiment will reach the LSND sensitivity in one or two years ||379|| . The 
proposed — > v e SBL experiments BooNE ||387|| at Fermilab, 1-216 ||388|| at CERN, 



ORLaND [|389|1 at Oak Ridge and NESS at the European Spallation Source [[39011 plan 
to achieve a sensitivity about two orders of magnitude better than the one of the LSND 
experiment. 

It is now widely believed that the recent measurement of an up-down asymmetry 
of atmospheric /z-like events in the Super-Kamiokande detector represents a strong evi- 
dence in favour of neutrino oscillations generated by new physics beyond the Standard 
Model (see, for example, Ref. [432]). We think that in order to understand the origin 



of neutrino masses and mixing it is necessary to investigate in much more detail all the 
phenomena generated by these properties: neutrino oscillations, neutrinoless double-/? 
decay, distortion of the high-energy part of the electron spectrum of tritium (3 decay, 
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neutrino magnetic moments, etc. There is no doubt that the next generation of neutrino 
experiments will lead to a great progress in the understanding the properties of neutrinos, 
the most puzzling among the known particles. 



Note added 

After this review was finished and presented to the editor (December 15, 1998), new 
Super-Kamiokande data were reported at WIN99 (Capetown, January 24 - 30, 1999) 
by Y. Suzuki and at the VHIth International Workshop on Neutrino Telescopes (Venice, 
February 23 - 26, 1999) by K. Scholberg and K. Inoue. Also Soudan-2 published an 
updated result. Below we will shortly present the new development. 



Atmospheric neutrinos 



For the double ratio R (|5.8j ) the new values based on 736 days for FC events and 685 
days for PC events are 



R 



0.67 ±0.02 ±0.05 sub-GeV, 
0.66 ± 0.04 ± 0.08 multi-GeV. 



For the up-down asymmetry ( 5.10 ) of /z-like events in the multi-GeV region the new value 



= -0.311 ±0.043 ±0.01 

is more than 7 standard deviations away from zero. Thus, the new data confirm the Super- 
Kamiokande evidence in favour of oscillations of atmospheric neutrinos. The updated 
dependence of the up-down asymmetry A as a function of the lepton momentum for 
e-like and /i-like events is given in Fig. A. The new allowed region in the plane of the 
oscillation parameters under the assumption of — > v T oscillations is shown in Fig. B. 
The region allowed by Kamiokande and the allowed region found from the analysis of 
previous Super-Kamiokande data (see Fig. |5.1|) are also depicted. As it is seen from 
Fig. B, the new Super-Kamiokande data favour higher values of Am 2 than the previous 
ones. The best-fit values of the parameters are 

sin 2 2tf = 1 , Am 2 = 3.5 x 10~ 3 eV 2 . 

This fit corresponds to x m m = 62.1 for 67 DOF (the previous best-fit value of Am 2 was 
2.2 x 10" 3 eV 2 ). The new allowed range for Am 2 at 90% CL is 10~ 3 eV 2 < Am 2 < 
8 x 10~ 3 eV 2 . 

The Super-Kamiokande Collaboration also measured the zenith angle dependence of 
the number of up-going muon events (for a discussion see Section 5.1.2). Also these data 
are compatible with the hypothesis of —>■ v T oscillations. The best-fit values of the 
oscillation parameters obtained from the data on the angular dependence of through- 
going muon events are sin 2 2$ = 1 and Am 2 = 3.2 x 10~ 3 eV 2 (corresponding to x min = ^ 
for 8 DOF). 
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Figure A. The up-down asymmetry A of e-like 
and /i-like events as a function of the lepton mo- 
mentum. The hatched region shows the theo- 
retical expectation without neutrino oscillations. 
The dashed line for /j,-like events represents the fit 
of the data in the case of two-generation — > v r 
oscillations with Am 2 = 3.5 x 10 _3 eV 2 and 
sin 2 2d = 1.0. 
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Figure B. The region allowed at 90% CL in the 
sin 2 2-d-Am 2 plane for Vn — * v T oscillations ob- 
tained by Super-Kamiokande (December 98). For 
comparison, the allowed regions found by the 
Kamiokande Collaboration and by the Super- 
Kamiokande Collaboration (June 98) are also 
shown. 



The atmospheric neutrino FC and PC data can also be described if one assumes that 
v s transitions take place. In this case the best-fit values of the oscillation parameters 
are given by sin 2 2$ = 1 and Am 2 = 4.5 x 1CT 3 eV 2 with Xmin = 64.3 for 67 DOF. 

The Super-Kamiokande Collaboration reported also the observation of 231.8 NC 7r° 
events (see Eq.( p.llj )), after subtraction of background, with two e-like tracks with an 
invariant mass in the range from 90 to 180 MeV. As a preliminary result it was found 
that 



7T°\ / /7T° 



e ' data ' \ ° / MC 



1.11 ± 0.06 (data stat) ± 0.02 (MC stat) ± 0.26 (syst) 



In the case of —>■ v T transitions this ratio should be equal to one, whereas in the case 
of Up — > v s transitions this ratio should be about 0.75. The large systematic error in the 
present data does not allow to reach a definite conclusion on the transition channel. 

The updated value for the double ratio R of the Soudan-2 experiment is R = 0.64 ± 
0.11 ±0.06 (W.W.M. Allison et al, |fiep-ex/9901024l ). 



Solar neutrinos 

In the region of recoil electron energies between 6.5 and 20 MeV, during 708 days, the 
Super-Kamiokande Collaboration has observed 9531 events (about 13.5 events per 
day). For the value of the 8 B neutrino flux it was found $ff = (2.42 ± 0.04 ± 0.07) x 10 6 
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cm 2 s 1 . The ratio of this flux to the flux predicted by the SSM ||304j| is equal to 0.470 ± 
0.008 ±0.013. 

Now Super-Kamiokande has also data with an energy of the recoil electrons in the 
range from 5.5 to 6.5 MeV. This lowering of the electron energy threshold is the most 
important difference of the new data with respect to the previous ones. The best fit to 
the recoil energy spectrum is given by vacuum oscillations with the best-fit parameters 
sin 2 2$ = 0.80 and Am 2 = 4.3 x 10~ 10 eV 2 . This is mainly because of the observed 
excess of events with respect to the SSM prediction in the high energy region of the recoil 
electron spectrum. 

The contribution to the recoil spectrum in the region E e > 14 MeV comes from hep 
neutrinos, which according to the SSM is small (see Table 5.2). If the flux of hep neutrinos 
is considered as a free parameter, then from a fit of the Super-Kamiokande data under 
the assumption of absence of neutrino oscillations it was found = (M.O+q^) 
The contribution of hep neutrinos in the Super-Kamiokande recoil energy spectrum was 
analysed in detail in the papers of J. Bahcall and P. Krastev, Phys. Lett. B 436, 243 
(1998) and G. Fiorentini et al, Phys. Lett. B 444, 387 (1998). 

Super-Kamiokande also presented a new measurement of the day-night asymmetry 
with the result N/D - 1 = 0.060 ± 0.036 ±° f 2 8 7 . 
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A Properties of Majorana neutrinos and fields 

In this Appendix we summarize the main properties of Majorana particles associated 
with Majorana fields. These particles have spin 1/2 and all charges equal to zero. A free 
Majorana field v(x) satisfies the Dirac equation 

(i'fd lt -m)j/(x) = (A.l) 

together with the Majorana condition 

v{x) = v c {x) = CV T {x), (A.2) 

where v c is the charge-conjugate field and C is the charge-conjugation matrix which fulfills 
the relations 

Ct^C" 1 = - 7# , , Ct = C-\ C T = -C. (A.3) 
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From Eqs. flA.l]) and ( |A.2| ) it follows that the Fourier expansion of the Majorana field is 

v{x) = [ J^— « r (P) ^ + 4(P) « P (p) e**] , (A.4) 

J v /(27r) 3 2po r t^i 

where the spinors u r {p) and v r {p) are related by 

v r (p) =C{u r {p)) T . (A.5) 

The operators a r (p) and aj(p) annihilate and create Majorana particles with momentum 
p and helicity r (particle = antiparticle for a Majorana field), respectively. 
Let us consider the quantities 

KOjV,,, (A.6) 
where v a and i/& are two Majorana fields and 

0j = 1, 7 M , <V, 7^75, 75, for J = ^> T ' ^> p - ( A - 7 ) 

Taking into account the minus sign that appears under interchange of two fermion oper- 
ators and the Majorana condition ( |A.2j ), we obtain 

W a OjV h =V- b COjC- 1 v a . (A.8) 

From Eq. ( |A.3| ) we derive 

CaJ/T 1 = -<v and CjJC' 1 = 75 • (A.9) 
With the help of this equation and Eq.( A.3|) it follows that 

Oj v h = TjjVb Oj v a (A. 10) 

with 

tfr = -l for j = V,T and ^ = 1 for j = S,A,P. (All) 
Thus, for a single Majorana field ^ = z/ a = Uj, we obtain 

^7^ = V/i and va\ p v = v<j\ p ^v = VA, p . (A. 12) 

This means that a Majorana particle has neither electric charge nor magnetic moments 
or electric dipole moments. 

For any fermion field there is a relation between its left-handed component and the 
right-handed component of its charge-conjugate field (see Section 2.2): 

^r = W)l) c . (A.13) 

In the case of a Dirac field, ip c is different from ip and, therefore, the left and right-handed 
components are independent, whereas from Eq. ([A.13|) it follows for a Majorana field v 
that 

Pr = {vlY (A.14) 
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and thus the left and right-handed components of a Majorana field are not independent 
(in other words, a Majorana field is a two- component field). 

Though Dirac and Majorana neutrinos differ in the number of degrees of freedom, 
it is nevertheless difficult to distinguish between the two in practice. The reason is 
that in the standard weak interaction Lagrangians ( [2.1| ) and ( |2.2| ) only the left-handed 
components un (£ = e, fi, r) of the flavour neutrino fields are present. These are connected 
with the left-handed components v^l of the neutrino mass eigenfields by the relation 
vil = J2k^k^kL (see Section 2). The mass eigenfields could be of Dirac or Majorana 
type. In the first case we have = v^l + ^kR, where the v^r are independent degrees 
of freedom, whereas in the second case we have vjf = v^l + {ykh) c - The only difference 
between the Dirac and Majorana case with respect to the standard weak interactions 
resides in the mass term. In the limit of vanishing neutrino masses the right-handed 
Dirac fields decouple from the standard weak interaction Lagrangians and the distinction 
between Dirac and Majorana gets lost. Then the physical fields are simply given by the 
In summary, we have arrived at the following conclusion ||433|| : With the standard 



weak interactions one cannot distinguish between the Dirac and Majorana nature of the 
neutrinos in the limit of vanishing neutrino masses. In other words, for the standard 
weak interactions, Dirac neutrinos with negative helicity and Dirac antineutrinos with 
positive helicity are the same as Majorana neutrinos with positive and negative helicities, 
respectively, in this limit. 

In neutrinoless double-/? decay, which can occur only with Majorana neutrinos (see 
Section 6), one needs to calculate the propagator 

(0\T(u eL {x 1 )i^ L (x 2 ))\0). (A.15) 

With the Majorana mass eigenfields Vk satisfying the relation z/jT = — VkC this propagator 
is obtained by the following steps: 

(0\T (u eL ( Xl )uJ L (x 2 )) \0) 



1 - 75 tt2 /mrr / / s,T, s\ , n \ 1 ~ 7 5 



k 



2 k m ki I 7^31-3 ——^ C - ( A - 16 ) 



d 4 p e -M*i-**) l- 7s 
(2vr) 4 p 2 -m\ 



Finally, we want to prove that the CP parity of a Majorana field is ±i. The CP 
transformation of a Majorana field u(x) is given by 

Ucp v{x) U C p = rjcp 70 v(xp) , (A. 17) 

where the phase factor r] CP is the CP parity of the field u(x), Uqp is the unitary operator 
of CP conjugation, x = (x°, x) and xp = (x°, —x). The relation ( |A.17| ) can be written as 

\J GV Cv T {x) U£ = -r]* Cpl0 CV T {x P ) . (A.18) 
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Using the Majorana condition ( |A.2| ) one obtains 

l7 cp v{x) = -V*cp 7o v(x P ) . (A.19) 



From the comparison of the expressions ( [A. 17 ) and ( |A.19j ) one finds that r/cp = —f]cp^ 
which implies that 

V 2 cp = "I (A.20) 
and rjcp — ±i. As it can be seen from Eq.( |A.lS| ) written as 

U CP v\x) C/cp = -Vcp 7o v c (x P ) , (A.21) 

in general the product of the CP parities of particles and antiparticles is equal to — 1. 
The relation ( |A.20| ) results from the fact that in the Majorana case there is no difference 
between particles and antiparticles. 
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